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ABSTRACT 

Algebraic  Methods  are  employed  for  frequency  response  studies  of 
physical  systems.   The  transfer  function  of  each  system's  mathematical 
model  is  composed  of  at  least  two  independent  variables  embedded  in  the 
coefficients.   These,  together  with  magnitude,  angular  frequency,  and 
phase  angle  are  the  variable  parameters  to  be  studied. 

Computer  programs  are  developed  for  graphically  representing  these 
variable  parameters.   Interpretation  of  these  results  is  made  to  evaluate 
the  utility  of  such  methods  in  ^K^ff'^ncT  analysis  studies. 
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1.    Introduction. 

The  dynamic  response  of  linear  systems  can  be  found  basically  by  two 
methods:   experimental  and  analytical-   Continuous  advancements  in  tech- 
nology have  precluded  the  feasibility  of  the  former,  relegating  its  utility 
to  hybrid  studies  at  best.   Increased  attention  has  therefore  been  focused 
on  the  analytical  techniques  available,  on  development  of  new  techniques 
where  necessary,  or  possible  improvements  in  known  techniques.   A  funda- 
mental dilemma  to  the  analytic  approach  is  the  synthesis  of  one  or  more 
linearly  independent  differential  equations  which  adequately  describe  the 
system  to  be  studied.   Whether  the  problem  consists  of  multiple  equations 
involving  eigen  values,  or  a  single  mathematical  statement,  the  system 
limitations  and  characteristics  are  revealed  by  a  thorough  study  of  the 
roots  of  the  characteristic  equation.   This  equation  is  generated  im- 
plicitly by  matrix  representation  of  a  set  of  simultaneous  first  order 
equations,  or  explicitly  by  a  rearrangement  of  a  single  higher  order 
polynomial. 

Many  sound  advances  have  been  made  since  the  turn  of  the  century  in 
the  realm  of  root  location  techniques  for  time  invariant  systems.   As 
control  engineering  emerged  from  its  embrionic  stage  and  became  embedded 
more  firmly  in  all  phases  of  engineering,  added  emphasis  was  placed  upon 
the  knowledge  of  systems  behavior  in  a  dynamic  atmosphere  consisting  of 
at  least  one  varying  parameter.   Among  the  more  important  achievements 
in  system  stability  analysis  are  the  Routh,  Mikhailov,  and  Nyquist  cri- 
teria.  Frequency  response  characteristically  has  been  greatly  utilized  as 
a  fundamental  system  specification  ,  and  the  graphical  techniques  of 
Bode  and  later  Nichols  have  become  widely  used  throughout  the  world  for 
both  analysis  and  design.   Many  other  graphical  tools  have  been  formulated, 


including  sophisticated  root  locus  schemes  involving  cascade  compensa- 
tion, and  various  forms  of  feedback.   Inherent  in  all  of  the  above,  how- 
ever, has  been  the  severe  restriction  of  evaluating  systems  behavior  as 
a  function  of  only  one  variable  parameter.   The  number  of  systems  which 
truly  fall  within  this  category  are  in  the  minority,  and  it  therefore 
becomes  necessary  to  employ  single  variable  techniques  on  multivariable 
problems.   The  result  is  a  non-optimum  compromise  between  system  sensi- 
tivity and  computational  requirements,  i.e.  ,  the  calculations  necessary 
to  accurately  predict  system  behavior  in  multivariable  environments  be- 
come unwieldly. 

Several  noteworthy  achievements  have  been  made  to  partially  overcome 
this  difficulty.   Extensive  application  of  the  analog-computer  has  been 
made  to  problems  suitable  for  its  use.   The  advent  of  the  digital  com- 
puter has  opened  a  magnitude  of  possibilities,  one  of  which  has  been  the 
emergence  of  state  space  theory,   A  large  potential  exists  in  a  relative- 
ly new  field  known  as  Optimal  Control  Theory,  which  is  founded  upon  the 
maximum  principal  approach  of  Pontryagen,  and  Bellman's  dynamic  programming 
method.   Stochastic  processing  and  predicting  filtering  theory  have  become 
very  influential  in  the  field  of  controls  due  in  large  part  to  the  accom- 
plishments of  Wiener  and  Kalman.   A  still  different  approach  has  dealt 
with  the  various  possible  algebraic  manipulations  of  the  characteristic 
equation,  hence  the  name  ALGEBRAIC  METHODS.   It  is  this  topic  to  which  the 
following  study  is  devoted.   Basically,  these  methods  employ  techniques 
which  permit  engineering  systems  to  be  designed  and  evaluated  by  simul- 
taneous variation  of  two  parameters.   Extension  to  one  or  more  additional 
parameter  variations  is  possible,  but  again  the  computational  requirements 
and  intelligent  interpretation  of  results  become  dominant  factors,  just  as 


LO 


in  the  expansion  of  single  parameter  techniques.   Little  imagination  is 
needed  to  impress  in  one's  mind  the  impact  that  this  tool  will  have  in 
optimizing  two  parameter  design  problems  as  its  capabilities  become  better 
defined,  and  engineers  achieve  proficiency  in  its  use. 

The  following  section  presents  the  rudimentary  mathematical  theory 
as  applied  to  electrical  control  systems  employing  two  variable  para- 
meters.  Primary  emphasis  is  placed  upon  formulating  the  necessary  rela- 
tions for  frequency  response  graphical  representations.   Subsequent  sec- 
tions are  devoted  to  the  generation  and  evaluation  of  these  displays  from 
the  standpoint  of  analysis  and  design.   Frequency  response  and  damping 
criteria  are  investigated  by  means  of  various  forms  of  derivative  feed- 
back.  Noise  limitations  in  physically  realizable  systems  preclude  the 
usage  of  derivative  feedback  in  excess  of  second  order  in  most  cases,  and 
for  this  reason  the  study  of  such  stabilizing  schemes  is  not  considered 
practical.   Typical  examples  are  presented  and  a  comprehensive  evalua- 
tion of  each  is  endeavored.   Graphical  representations  are  utilized  and 
analyzed  according  to  their  respective  value  for  information  content. 
The  concluding  section  recapitulates  the  most  important  limitations  and 
advantages  revealed  by  each  graphical  analysis  example,  and  postulates 
further  extensions  of  specific  developments  where  applicable.   A  general- 
ized version  of  computer  programs  employed,  with  an  explanatory  flow- 
graph  of  each  is  included  in  the  Appendix. 
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2.    Mathematical  Theory  of  Algebraic  Methods. 

Conventional  network  analysis  techniques  rely  heavily  upon  pole- 
zero  representations  in  the  complex  frequency  plane,  as  well  as  real  time 
transient  and  frequency  response  curves.   These  methods  are  not  conducive 
to  the  synthesis,  or  design  problems,  however.   System  parameters  are 
seldom  specified  precisely,  and  a  slight  variation  in  one  might  correspond 
to  an  excessive  alteration  in  another  while  still  satisfying  yet  a  third 
design  criterion.   The  single  parameter  variation  is  seldom  a  rewarding 
approach  to  this  type  of  problem,  and  it  therefore  becomes  necessary  to 
transform  response  curves  into  another  reference  frame  where  system  para- 
meters are  the  variables.   Fundamental  complex  variable  transformation 
theory  provides  the  necessary  conversion  media,  based  upon  Cauchy's 
Principle  of  Argument. 

2.1  Transformation  Techniques  Employing  Chebyshev  Polynomials. 

The  ensuing  development  is  based  upon  three  fundamental  steps,  and 
permits  transforming  the  complex  frequency  plant  into  a  coefficient 
plane  and  a  parameter  plane,  where  the  corresponding  variables  are  arbi- 
trarily designated  as  a  ,  a  and  alpha  (  OC     ),  beta  (  P  )  respectively, 
x   y 

The  first  step  involves  choosing  a  contour  to  be  transformed.   The  complex 
frequency  damping  constant  zeta  (  j  )  has  a  radial  line  contour  in  the 
s-plane,  and  is  selected  as  the  transformation  contour.   A  mapping  func- 
tion must  be  considered  next,  and  for  this  the  characteristic  equation  is 
appropriate.   The  third  stipulation  requires  a  transformation  variable 
which  allows  the  mapping  function  to  be  represented  by  two  or  more  inde- 
pendent equations.   In  this  case,  it  was  decided  to  separate  the  complex 
variable  into  its  real  and  imaginary  parts.   Listed  below  are  the  neces- 
sary steps  required  to  achieve  this  transformation.   Considerable  attention 
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is  devoted  to  a  complete  and  detailed  development  in  order  to  provide 
continuity  in  the  fundamental  relationships  which  might  be  readily  refer- 
enced as  major  divergence  points  for  certain  ensuing  pursuits.   The  basic 
relations  presented  in  the  following  diagram  and  adjacent  equations  serve 
merely  to  establish  a  consistent  notation  for  all  further  developments. 
The  approach,  although  exceedingly  elementary,  is  considered  necessary  to 
avoid  perpetuating  existent  ambiguities  where  similar  ground  rules  have 


not  been  employed. 
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Substituting   the   relations    for    o~  and   u)    into   the   top  equation: 

S  =  -    :  ::r>i  '  ~    I l        =    ul^COsQ  +  j  u^  5 in  6 

S  =    (  c  o  s  ©    +-   J 

Employing  the  results  of  the  mathematician  Euler,  this  last  equation  may 

JS 

be  written:     5  =  luT)    6 

And  using  DeMoivre's  Theorem  the  relation  for  any  power  of  S   is  found 
to  be: 


Sn   =    LU    "    e 


K     rtJK0 


uuv  K  (cos  k&  +  j  s'ir>  k  9  ) 


(2.1-1) 
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It  is  now  necessary  to  define  two  new  variables: 

TKf-5)    =    cos  k  0    =    coS(k  cos~'(-  I)) 
[j      /_  o    —  ALIL5  $    sc'tt  (kcqs (-  I)) 

Substitution  of  these  equivalences  into  (2.1-1)  yields: 

5K=     u.,[TK(-5)    +    J*fl-  5*    UkC-$)1 

Since  engineering  problems  involving  stable  systems  are  concerned  with 
positive  values  of  zeta,  it  is  advantageous  to  employ  the  following 
relations  which  are  obtained  by  trigonometric  identities: 

rh  H)  =  -iKTk(S)      ;     uK{-s)  =  -j""'  kk(0 

All  characteristic  equations  can  be  reduced  to  a  general  form  which  is 
herein  defined  as: 

Rewriting  this  in  summation  notation,  and  employing  the  identities 
developed  above,  we  can  write: 

Fes)  =  #    <kv  sK  -   2.    v  a-*'  4-  j  1   a.  wk  =  0 
F  ^)  =  £-    ^k    ^     cos  /<©  4-  J  £    aK  Wy,     sin  k&  -  O 

Fee)  =  £VK  ^W*T/f)+'jl  a.^V-O^jJ-O 

A  polynomial  may  be  set  to  zero  when  its  orthogonal  components  are 
independently  zero;  hence: 


l  \  "/h)*  TKa)=o 


■H  *%    l~"     '  H    <  S;  -O  (2.1-2) 

K  =  0 
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k+-/ 


%  "»    I'D  MK  r f  >   =  0 


(2.1-3) 


It  would  now  seem  plausible  to  inspect  several  values  of  these  new  vari- 
ables T^  (  £)   and  li  y,   (  §  )  ,  in  order  that  any  discernable  character- 
istic traits  might  be  utilized. 

T.m  -  (-\Y  T>(-$)  =  T»  (-U  =   cot(oe)-c>t(o)=l 

T,(5)  =  (-0'  f,  (-^)=  -"T  H)  =  eoi/eoj^Z-S)')  =  S 

ii.<s)-H>'u.H>-«={^=-^-  ° 

»-  (-ozm,  r-s)  =  »i:fc;-.f.{j  j  =  -i 

Comparison  of  these  polynomials  with  a  table  of  Chebyshev  functions  will 
reveal  their  identity.   Accordingly,  the  following  recurrence  relations 
may  be  employed  for  complete  tabulation. 

T,«.,($)     -ZJT»(:5)       *■    T*-,(S)      =     O 
U  «*■'($)     -I  s  U  k(H       +     U  k-i  C  f  3      =     o 

It  can  be  shown  that  T^  ( Q  and  U^(H  are  related  by  the  following: 

tk  ci)  =  s  «k  c  ?)  -  u*.,  en 

Substituting  this  into  (2.1-2)  yields: 


£**<U)"teUk(i)]   -ic\<(-l)Wm=0   (2.1-4) 


which  can  be  rewritten  as  follows: 
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Factoring  a  (-1)  out  of  (2.1-3)  and  multiplying  that  result  through  by 
(  £  )  in  a  separate  step  yields  the  following  sequential  equations,  so 
arranged  for  easy  inspection. 

This  last  form  of  (2.1-3)  is  identically  equal  to  the  first  term  of 
(2.1-4).   Therefore  the  two  fundamental  relations  become: 

^    H) 

K  =  0 
Y) 


H)   «KW/(i  „_,  ^()  =  0  (2.1-5) 

fl)K  %  V  UK  (U   =  0  (2.1-6) 


The  utility  of  (2.1-5)  and  (2.1-6)  is  not  obvious,  but  their  judi- 
cious use  results  in  a  powerful  tool.   Employment  of  the  bilateral  trans- 
formation allows  the  natural  angular  frequency  (  UJ^  )  to  range  from 
minus  infinity  to  plus  infinity,  thereby  eliminating  it  as  an  unknown 
variable  per  se.   Choosing  a  specific  value  for  zeta  (  ?  )  allows  the 
values  of   U  ($)   to  be  obtained  by  the  preceding  method  or  by  refer- 
ence to  precomputed  tables  [1],   Thus  any  two  values  of  Q-^    can  be  al- 
lowed to  vary,  thereby  formulating  two  independent  equations  in  terms  of 
two  variable  coefficients  of  the  characteristic  equation.   Mitrovic  [2], 
who  pioneered  the  above  development,  chose  the  two  lowest  order  coef- 
ficients, (X0    and  CI,  ,  for  his  study.   Elliott,  Heseltine,  and  Thaler,  [3] 
and  summarily  Siljak  [4]  extended  Mitrovic's  approach  to  additional  co- 
efficients.  A  brief  summary  of  the  coefficient  plane  derivation  appears 
on  the  next  page. 
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2.2  Coefficient  Plane 

For  the  generalized  form  of  the  characteristic  equation,  let  two 
variable  coefficients  be  defined  as  CX^  and  (X  y  .  The  order  of  the 
equation  (  77  )  is  no  smaller  than  X  ,  and  X  is  greater  than  y  . 
The  value  of  y  is  non-negative.  n  ^  X  >  Y^O)  Transposing  the  <X^ 
and  (Xy    terms  in  (2.1-5)  and  (2.1-6)  yields  the  following  results: 


(-'J'  A*   <V  Ux_,  (i)     +    ('l)\y      U^  Y    Uy.f     (I) 


H  =  c 


r-o  *  aA  u)^ ux  ($)  +  f-/)  ^/  ^  *'  «y  f  s) 


K7^X 

My 


Solving  the  above  equations  simultaneously  for  ax   and  <X  y  yields: 


71  ,  Jf-X   K-X 

k  =  o 
k^x 


(2.2-1) 
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a_    =  1L    (-0   W  a, 


iW$)<VO  -KyWfc-iW 


Ux  fS)M^f5>-«x./5)My^ 


(2.2-2) 


It  will  now  be  shown  that: 


Ux-vti)   =    Wy(J)  U^CfJ  -11^  (0  Wxfi)   (2-2-3) 


Let   X  =  2 

From  the  preceding  Chebyshev  recursion  relations: 

Ms)  =  i 

Ms)  =  z  J  u,($)-u.(s)  =  Z$(0-o-lf 


Therefore: 


U^     H)  =  ^y  fU  --  2  S  u    ;  (S) 


(2.2-4) 


But  it  has  been  shown  that: 


U,(l)  =-  UK  (S) 


(2.2-5) 


Substituting  (2.2-5)  into  (2.2-4),  we  see  the  general  recursion  formula 
as  previously  stated  for  Chebyshev  polynomials: 
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u.y ah-  Z;iay_,  a)   +  uy_2  n)  -  O 

It  can  be  shown  by  inductive  reasoning  that  the  above  relation  exists 
for  all  values  of  K   [5].   Therefore  we  can  also  write: 

Thus  the  general  solution  for  two  variable  coefficients  becomes: 


ax  =  —  Jl   (—1) 


uK-x  (S) 


7) 


Ux-y  CS) 


A  graphical  display  of  the  characteristic  equation  is  often  formed  with 

as  one  axis  and   -'■   as  the  other  axis.   Such  a  plot  is  termed  a 
coefficient  plane.   With  sufficient  data  and  interpolation,  it  is  pos- 
sible to  derive  a  very  satisfactory  design  in  this  manner  for  systems  whose 
characteristic  equations  possess  the  desired  variables  in  coefficient  form. 
An  excellent  treatise  on  the  applied  coefficient  plane  can  be  found  in 
the  NASA  Contractor  Report  by  Thaler,  Siljak,  and  Dorf  [6]. 
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2.3  Parameter  Plane,  Linear  Case. 

The  coefficient  plane  provides  a  means  of  working  with  two  variable 
coefficients  of  a  polynomial,  as  long  as  each  of  these  variables  does 
not  appear  in  more  than  one  coefficient.   It  can  be  recalled  that  each 
coefficient  within  a  polynomial  is  merely  the  sum  of  the  product  of  the 
roots  taken  TO   at  a  time,  where  Tq    equals  the  order  of  the  equation 
minus  the  order  of  the  particular  coefficient's  independent  variable. 
Therefore  severe  restrictions  are  encountered  when  attempting  to  use 
this  method  for  simultaneous,  but  independent  variation  of  any  two  roots 
of  a  polynomial.   Thus  the  coefficient  plane  method  does  not  include 
the  more  general  class  of  problems  which  have  the  desired  variables  em- 
bedded within  several  coefficients.   Siljak  [7]  postulated  that  this 
class  of  problems  could  be  handled  in  a  manner  similar  to  that  of  the 
coefficient  plane,  provided  these  variables  appeared  linearly.   Desig- 
nating the  variable  parameters  as  alpha  (  oL     )  and  beta  (  8      ),  one 
can  rewrite  any  coefficient  of  a  general  polynomial  as  follows: 

a  K  ^  k  v   <*.  ■+■  C.  (3  -f-  d  K 

Using  this  expanded  form,  the  two  independent  equations  (2.1-5)  and 

(2.1-6)  become: 

*  £     (~\)      t     uu        U       (i) 

k=0        k        n  k-i 


+  i    (-1     c      u,         IX       (I) 

k=0  k       ^  *-/ 
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*$..  f-J-])  c«  a«a) 


'    *  fO  '*  "/  M«  =° 


This  rather  unwieldy  notation  will  be  slightly  modified  by  the  follow- 
ing substitutions: 
Let: 

■■'<,*„>  =  1h)*  a„  «,*  uK  (sj 

S,*„)  -=1    (-)*    ck    yJ/  ^  (l) 

fy-0 

Employment  of  these  substitutions  reduces  (2.1-5)  and  (2.1-6)  to  the 
following  form: 
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These  two  equations  are  quite  easily  solved  for  alpha  and  beta  by  Cramer's 
rule,  or  by  employing  the  algebraic  method  of  elimination.   The  resultant 
solutions  are: 


B,  (%>&»)  C^Cl,^)    -  B%  (i,  uJ^)  c,  (l>u)n)   (2'3_1 


The  stability  variables  most  commonly  employed  by  control  engineers 
are  the  damping  ratio,  zeta,  and  the  undamped  natural  frequency,  (m^). 
Most  of  the  current  literature  places  great  importance  upon  these  terms, 
and  interprets  system  performance  by  correlation  with  the  standardized 
second  order  system.   In  some  cases,  however,  the  rectangular  complex 
frequency  plane  configuration  consisting  of  the  imaginary  frequency 
sigma  (  CT  )  and  the  real  frequency  omega  (  U'  )  is  more  useful  than 
the  polar  coordinate  system  with  damping  ratio,  zeta  (   f  )  and  the 
undamped  natural  frequency  (  oU^).      In  second  order  system  parlance, 
sigma  represents  the  rate  of  exponential  decay  of  the  response  signal, 
while  omega  is  the  angular  frequency  at  which  this  signal  oscillates. 
The  necessary  steps  for  representing  (2.3-1)  and  (2.3-2)  in  terms  of 
sigma  and  omega  involve  the  use  of  the  binomial  expansion.   A  brief 
outline  of  this  method  is  now  presented.   Recalling  the  basic  relations 
between  <T   ,  U)     ,   f   and  U) y,     in  2.1,  it  is  possible  to  rewrite  (2.1-1) 
as  follows: 
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S     =   uy  77      e  =      (  <r  +  J  ix) ) "    =    XH   +  j  yk 

where    the   individual   components   of     VK      and     )K      are   seen  to  have   a 
particular   recurrence   pattern  as   shown  by   the   following  tabulation: 

^°  =  i 

Sx  =  a-1  -LUl     +    Zjlu 

S^rrrcr?     -   3  r  ^      +  j    (     3  <r  *  tU    -     UJ3) 

$**—    cr^-^cr1^2--*-^^    +■  j    (  V  c  2  uJ    -  H  <ru>  ~*  ) 

5  s"  =    cT r  -  /  0  cr2u;3  +  5  <r  uoH  +  j  £  5"  0-  *  u)  - 1  O  <r260  3  +  nJ  *"  ) 

From  which: 

X  I 

X,     —    cr 

2  2.  * 

X3     =    cr3  -  3  cr  K.JJ2- 

G 

y,  =  u) 

yz  =  2  <r  u) 

y,  =  3  a-  z  uj    -  ou3 

yH  =  &<r  *  uj    -  h  <ruo  3 

which  can  be  consolidated  into  binomial  expansion  form  for  the  general 
ized  case. 
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+   J    LrO  LzIh)  <r  uj  J  ^  —  0 


where      I   X  /    =    -7- — tV      rrr       for      K  - 


k   1 

(k-x)  !  (x)! 

and    (  £)  =  Q  whenever  ft  £  X 

Recalling  the  generalized  form  of  the  characteristic  equation  r(S)  =  £^S, 

we  can  again  separate   $  s   into  its  perpendicular  components  X,      and  yK 


n 


K 

srms 


such   that      fr_      ®~KXu         represents   the  summation  of   the   real   tei 
and   4>    n     \/  denotes   the   summation  of  all   imaginary   terms.      Thus: 


K-c 


I   *KxK  =  I    I    «ltLM(tl)<r"-ziu,'-iJ 


K=o  k=o    <=o 


«-c       " 


In  order  to  effectively  utilize  these  equations  for  computer  applica- 
tions, it  is  necessary  that  the  recursion  relation  be  explicitly  formu- 
lated.    This  can  be  achieved  by  observing  the  following: 
S  =   (T  +  J  W   =^>    S  -  cr  -  j  to    =    Q 
S  =   <r  -  j  LU      =^>     5  -o-^j^  =  Q 
multiplying  these  two  equations; 
(S-<r-}uj)(<;  -  (T  +j  UJ)  =  0 
Expanding  this,  we  have: 

Sz   -    Z   cr-s    -h    cr*   +    luz  =  O 

Pre-multiplying  this  by    3  ^     we  may  write 
5K    [   SZ  -ZrS     +  <r2    4-  uu  *■  ]  =  O 
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which,  upon  expansion,  becomes 

SK  +  Z  -  Zr  s"  +  l    +-  <r*  s*  +  u,2-   5K  =  O 

This  may  be  rewritten  in  terms  of  orthogonal  components  where 

s    "      =    xK^z    +  j  yK^2 
s  *       =    '*  k        +  i  y. 


+.  r  <r  *>  w  z)  l xK  +■  J  yK )  =  o 

Separating  this  into  its  real  and  imaginary  parts  results  in 

Noting  that  (T  is  X,   ,  and  vJ    is  y,   ,  we  can  make  these  substitutions, 
resulting  in  the  following  explicit  recurrence  relations: 


Thus  one  can  solve  a  given  coefficient  equation  with  coefficients  of  the 

form  aK  =  Lo(,  ^^kp^^*    for  the  Parameters  alpha  and  beta  in  terms 

of  the  variables  sigma  and  omega.   The  method  is  the  same  as  for  the  case 

involving  the  variables  zeta  and  uO  ~     .   The  two  independent  equations 

become : 

■n  y>  n 

,'-         CK     *k         +       : 

k 


a    < '    ■  tK   yK    <-    /s  I    £*   /«      +•    £     <^k  /«    =  o 


25 


which  may  be  written  in  abbreviated  form  as 
where : 


B, 

1 

b* 

*k 

c, 

= 

I 

«K 

*K 

D, 

== 

e(K 

*K 

&2.= 

£    S  * 

cx= 

1  c*  * 

Dj,= 

jjL  ^  ^ 

As  in  the  preceding  development,  the  solutions  for  alpha  and  beta  take 
the  form: 


ot  =  dC^jiu)    b^CcTjUj)    -    C^(<r,u))  D,(<r;uj)  (2  23) 


8  xt      Pi  (^^>  &x(<n  LV)  .  -.  P*.L<rJw)Bit<nw)  (2 . 3  _4) 

These  equations  are  seen  to  be  of  the  same  form  as  (2.3-1)  and  (2.3-2). 
Thus  the  method  of  solution  is  identical  to  the  one  using  Chebyshev 
polynomials,  with  only  the  recurrence  relation  and  the  solution  variables 
being  different. 

An  analogous  approach  to  this  problem  does  not  involve  the  Chebyshev 
polynomials  method  at  all,  as  can  be  seen  by  the  following  procedure 
which  starts  by  employing  (2.1-1)  directly: 

F  (  S)  =     ji    ak  S  *    =     2l  aK  uUyi  *  cos    k  G 


+  J  £ .    «.K  ^h  K    sin   k  G        --  O 
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K     -h 


Separating  this  polynomial  into  its  orthogonal  components  we  may  write: 

iu{fc$)}  =   i    Q_K     uoTl*  ccs    ke  =  o 

Now  substituting  the  integral  components  of  a    into  each  of  the  above 
yields  the  following  general  terms: 

H 


y>  r, 


o(.  £.    bk  ucr  '  ca  k&     +  8  JL  ck  ijj     cos  ke 


+-        £      d  k    uu'y.     ccs  k  &      —     O 


'n  ft 


k=c     r*      "  > 


n 


k 


Solving   these   equations    simultaneously   for  alpha  and  beta,   we  have: 


,£     cv  ^r;      cos  >f^-£  ^^rj   s<n/r6  -(£cifa;r|,st*y7/f^£.c/fcp«;w  cask© 
^   _     *-<» *Z1 *12. *f* (2.3-5) 

H  =  c  K'-O  "=o     '  fro 


k  v  k  k 

£_AV^   coske£  b^j»  sink0-^.dtiuJnisinKO£_  bKu)n  caskB 
fi_     h«q    K *=o    ' £=£ *££ (2.3-6) 

iL  h    u/'n  cc<k©£cKu;n  sinktf-Z-Pfrtc,,  ''"*#  ^S^  ^  ^0 
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These  two  alternative  parametric  solutions  are  seen  to  be  no  less  compli- 
cated in  form  than  their  Chebyshev  polynomial  or  binomial  expansion  counter- 
parts.  In  fact  their  application  is  more  cumbersome  since  there  is  no 
analogous  recursion  equation  for  the  trigonometric  functions.   Therefore 
the  first  approach  resulting  in  equations  (2.3-1)  and  (2.3-2)  is  utilized 
in  all  but  the  most  specialized  cases  since  it  was  the  first  method  to  be 
successfully  developed,  and  its  application  is  more  familiar  than  the  sigma- 
omega  form. 

EXAMPLE  1  (FEEDBACK  COMPENSATION  FOR  LINEAR  THIRD  ORDER  SYSTEM) 

At  this  point  it  will  be  advantageous  to  select  a  basic  system  for 
initial  studies.   The  linear  case  is  sufficiently  unencumbered  by  mathe- 
matical manipulations  to  properly  illustrate  certain  fundamental  character- 
istics, and  therefore  serves  as  an  appropriate  introduction  to  the  forth- 
coming discussions  of  section  4.   Consider  an  uncompensated  system  with  a 
forward  transfer  function  of  third  order,  whose  gain  is  sufficiently  large 
to  make  it  unstable.   Such  a  system  is  depicted  by  the  unity  feedback 
block  diagram  and  characteristic  equation  shown  in  figure  2-1.   The  corres- 
ponding root  locus  diagram  is  included  in  figure  2-2  and  it  is  seen  by 
Routh's  Criterion  that  the  limit  of  stability  is  at  K"  =  30,000.   For 
larger  gains  than  this  critical  value,  it  will  be  necessary  to  introduce 
either  cascade  or  feedback  compensation  to  regain  system  stability.   Selec- 
tion of  tachometer  and  acceleration  feedback  results  in  the  compensated 
block  diagram  and  characteristic  equation  shown  in  figure  2-3.   Choosing 
a  forward  gain  value  of  100,000  results  in  a  highly  unstable  uncompensated 
system.   The  stabilizing  effect  of  either  one  of  these  forms  of  deriva- 
tive feedback  without  the  other  is  shown  by  the  cross  hatched  root  reloca- 
tion zones  of  figure  2.4.   The  effects  of  simultaneous  variation  of  both 
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forms  cannot  be  as  easily  predicted  on  the  root  locus  plane.   At  best,  one 
can  generate  several  root  locus  families,  each  as  a  function  of  only  one 
variable  parameter,  and  gain  from  these  patterns  an  approximate  idea  of 
the  required  parameter  values  needed  to  place  the  dominant  roots  of  the 
system  in  a  prescribed  position.   The  parameter  plane  technique  becomes 
a  very  useful  tool  in  this  type  of  situation,  however,  since  the  pre- 
scribed root  location  automatically  specifies  the  required  value  of  zeta 
and  also  of  wJn  .   A  cursory  inspection  of  equations  (2.3-1)  through  (2.3-4) 
will  reveal  that  the  following  curves  can  be  drawn  in  the  parameter  plane  ? 
which  by  convention  has  alpha  on  the  abcissa  and  beta  on  the  ordinate: 

1.  Constant  zeta  curves,  in  which  zeta  is  held  fixed  at  a 
specified  value,  and  cvp     is  allowed  to  vary  through  a  pre- 
determined range. 

2.  Constant  undamped  frequency  or  "natural  frequency"  (  tu^   ) 
curves  in  which  uu^     is  held  at  a  specified  fixed  value, 
and  zeta  is  allowed  to  vary  from  0  to  1. 

3.  Constant  C  —  uo-n  curves,  in  which  zeta  and  uJ^  are  both 
varied  such  that  their  product  is  always  constant  at  pre- 
selected curve  values. 

4.  Constant  sigma  (  C  )  curves,  which  represent  lines  paral- 
lel to  the  imaginary  axis  in  the  complex  frequency  plane. 

5.  Constant  "damped  frequency"  (  uJ    )  which  represent  lines 
parallel  to  the  real  axis  in  the  complex  frequency  plane, 
and  are  constructed  by  holding  u)     fixed,  and  varying  sigma 
from  zero  through  a  pre-selected  negative  value  range. 

6.  Constant  real  root  (  o~   +J°)  curves,  which  represent 
specified  fixed  points  on  the  negative  real  axis  of  the 


29 


complex  frequency  plane.   (Thus  it  is  seen  that  certain 
points  on  the  s-plane  are  transformed  into  lines  within 
the  parameter  plane.   This  factor  plays  a  predominant  role 
in  singular  line  theory). 
Not  all  of  these  curves  need  be  used  by  any  means.   The  combination 
of  any  two  of  the  first  five  curves  positively  identifies  a  complex  root 
in  the  left  half  of  the  complex  frequency  plane,  while  the  sixth  type  of 
curve  locates  the  real  roots  on  the  negative  sigma  axis  of  the  s-plane. 
Figure  2-5  is  one  type  of  parameter  plane  representation,  in  which  Ka    ■ 
alpha  and  K±   =  beta.   It  shows  the  constant  zeta,  constant  U)v    and  real 
root  curves.   One  can  see  that  for  dominant  roots  at  5  =  -16.5  +  j33  the 
required  values  of  alpha  and  beta  will  be  at  the  intersection  of  the  £  = 
0.5  and  UJ^    =  33  curves.   In  this  example,  they  intersect  at  alpha  = 
0.0005  and  beta  =  0.036.   The  third  root  is  real  and  can  be  determined  by 
the  real  root  curve  which  coincides  with  these  values  of  alpha  and  beta. 
It  is  seen  that  the  real  root  is  approximately  at  S  =  -80.   These  values 
are  indicated  by  point  A  on  figure  2-5. 

The  curves  for  this  example  are  deceivingly  simple,  and  therefore 
easy  to  interpret.   Additional  care  must  be  employed  in  higher  order  equa- 
tions, where  all  the  root  locations  are  not  so  easily  defined.   In  order 
to  ensure  the  desired  dominant  characteristics,  one  must  certify  the  loca- 
tion of  all  existent  roots.   In  exceedingly  complex  cases,  this  might  neces- 
sitate a  second,  enlarged  parameter  plane  plot  centered  at  the  point  in 
question.   This,  however,  is  a  greatly  simplified  task  compared  to  single 
parameter  techniques. 

2.4  Paramater  Plane,  Non-Linear  Case. 

Equations  (2.2-6)  and  (2.2-7)  have  been  developed  to  handle  two  variables 
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when  they  are  individual  coefficients  of  a  polynomial.   Equations  (2.3-1) 
and  (2.3-2),  or  similar  forms,  apply  to  two  variables  embedded  linearly 
in  the  polynomial  equation.   The  following  development  is  concerned  with 
non-linear  arrangements  of  two  variables  within  the  coefficients  of  a 
given  polynomial.   Consider  first  the  case  where  the  variables  appear  in 
product  form.   This  arrangement  is  termed  the  alpha  -  beta  product  case, 
because  the  coefficients  of  the  polynomial  can  be  represented  as 

aK  =  K.  oc  4-  cK  &    -*-  d  K  oc(3  +  ek 

Substituting  this  into  (2.1-5)  and  (2.1-6),  one  can  write: 


ol  £    (-/)KbK'Vu^(-  (I)    +•/?£   (-ifcxVJ^U^Cs) 


n 


4-  £  f-/)*ek  u;/  UK1,  ($    =  0 


4-  oC(3  £  Hi*  dK  ^r>*  U  K  f  J) 


4-£f-/)*eK  */«H  a)  -0 


Employing  notational  substitutions  similar  to  those  in  the  preceding  sec- 
tion, we  may  write: 

K  =  0 

K= 
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p,  =  %  \K  tf;vft) 


Dz  =. .^   ^  *  ^  .  liK  (  S ) 


E.  =  ^   e.  uJnr  'tf  K_,  CO 

thereby  reducing  the  preceding  equations  to  the  form 

The  functional  subscripts  have  been  included  up  to  this  point  to  emphasis 
the  applicable  dependent  variables.   It  will  now  be  necessary  to  elim- 
inate this  descriptive  subscript  notation  in  order  to  achieve  better 
comprehension  of  mathematical  manipulations.   Henceforth  the  dependent 
variables  will  be  regarded  as  zeta  and  LO^     unless  otherwise  specified, 
and  their  inherent  inclusion  will  be  assumed  in  the  remainder  of  this 
section.   Accordingly,  the  above  equations  may  be  rewritten  as 

ol   B,  +■  (3  C,   +■  oLp    p,  +•  £*,  =  0 

The  solution  of  these  simultaneous  equations  for  either  alpha  or  beta  is 
possible  provided  their  Jacobian  is  not  zero,  where  the  Jacobian  is  de- 
noted: 

j      (b,  +  p  oHCt+oLO j    -fc,  ■*' oca) (Sjl  +  p c J 
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Presuming  this  to  be  the  case,  the  solution  for  beta  is  found  as  follows 

[-B,  -pojC«  b2  +  Pcz  +  h2  ;j  =o 

which  may  be  written  in  determinant  notation  as 


Acc  P2    -   (Ace    *-AEB)(S     +    A 


E  B 


0 


where 


^CD     = 


c,       cx 
o,      a, 


A 


CB 


B,  Bz 


and  so  forth.   Employing  the  quadratic  formula,  beta  is  found  to  be 


2-Acz  (2.4.1) 

Proceeding  in  a  similar  fashion,  it  is  possible  to  solve  for  alpha  by 
the  following  process: 

from  which: 
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or  in  determinant  notation: 

The  desired  solution  for  alpha  then  becomes 


^  =  -  (Asc    +AED)  ±  -j  (Aec   +  A^p)*  -*(AepAec)  (2.4.2) 

While  the  above  solutions  can  be  used  to  obtain  alpha  and  beta,  it  is 
advantageous  to  simplify  them  for  adaptation  to  computerized  programming. 
Analyzing  the  terms  under  the  radical  signs,  it  is  possible  to  show  their 
equivalence  as  follows: 


Noting  that  ACg  '==:^   —  A  or         >  and  squaring  both  sides  yields 
which,  upon  expansion,  becomes 

(Aec)Z  +   l(ABC  A£0)    +    (&£0)Z     -^AeDAec     = 
(Aec;1  -Z(ABCAeo)    +(Aeo)1-HACDAea 

Transposing  and  collecting  terms 

4ABCA£0     -^AeDAec     +^ACDAt-e    =0 

Dividing  through  by  4,  and  expanding  determinants, 

(BlCz-BlCl)(tlDl-E2.Pl)  -  (B,Dj.-BiPl)(£iC2.-E2.Cl) 

4-  (C>DZ  -CZD,)(  E,BZ-£ZBI)  =  0 

B.CjlDj.E-,  -BzCt  OzEl-BlCzDlB1-hBzClP,  Ez  -B,CzPzEi  +  BZ_CZP,E, 
+  B, C,  D2  tz  - S^C,  /),  F2  *  fi^  C,  ^ E,  -  /?;,  £z 0, E,  -B,  C, 0XEZ+B, CzOtEz  =  0 
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Thus,  we  may  define  the  terms  under  both  radicals  as  K    ,  where 


and  the  solutions  for  alpha  and  beta  become 


Oi    = 


-    ~fcw  +    A£D  )  ±  R 


A 


a    =  —  r-Aflr,  +  AED)  ±  H 

of  the  four  possible  combinations  for  alpha  and  beta,  it  is  necessary  to 
determine  which  of  these  satisfy  the  original  simultaneous  equations. 
The  following  results  are  obtained: 

x  +  C,  (3  +■  0,   ol  (3    4-  Ft.  j  =  0 
ll*i-  Cx  p    +    Dz  ot  (3    +■  ETZI/  =  0 
(/?,  p?   -BiD,)*  +  (c,  Pz  '-CzDt ')  (3  -f-  EiPx-B7_D.   =  <9 

In  determinant  form  this  may  be  written: 

*  +  ACD  0     +  AtrD     =0 

Substituting  the  above  relations  for  alpha  and  beta 

2  ~Z  fcL 

Therefore,  the  radical  coefficients  for  one  term,  either  alpha  or  beta, 
must  be  the  negative  of  the  radical  coefficient  for  the  other  term.   This 
permits  the  following  solutions  for  alpha  and  beta: 


0Ck=-   i *2-     r-^ l  '      n  ;  U=\,l  (2.4-3) 


(4 

ED 

+- 

A 

EC)  + 

!- 

o* 

R 

Z 

^BD 

r* 

Et 

- 

A 

ec)  + 

f- 

•l)" 

R 

Pk==      t    -u       ~ ;   k==|jZ     (2.4.4) 


Z  A 
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SINGULAR  POINTS 

Inspection  of  (2.4-3  and  (2.4-4)  reveals  that  limitations  are  en- 
countered whenever  ^  bd   or  ^cd     equal  zero.   For  these  cases,  it  is 
necessary  to  resort  to  the  original  set  of  simultaneous  equations  for 
the  complementary  solution  as  follows: 

Should  Agp  be  zero,  the  two  solutions  for  beta  are  obtained  from 
(2.4-4)  and  then  substituted  into  the  original  simultaneous  equations 
where  the  solution  for  alpha  is  obtained  by  the  following  relationship: 

*y=P»C'  +  E'       =   4"C^Et-   ;  U  =  \,l  (2.4-5) 


When  ACD  equals  zero,  the  analogous  approach  yields  two  solutions 
for  alpha  from  (2.4-3)  which  are  likewise  substituted  into  the  original 
equations,  yielding  the  following  values  of  beta: 

p        **gi  +e»        =  j^e^f^   .   k  =  (Z    (2.4-6) 


SINGULAR  LINES 

The  solution  of  a  set  of  TTJ   simultaneous  equations  in  77  unknowns 
is  not  always  obtainable.   When  a  solution  does  indeed  exist,  it  some- 
times possesses  unusual  characteristics  in  that  it  is  not  a  unique  solu- 
tion.  These  cases  have  been  catagorized  as  singular  lines,  and  their  ap- 
plications to  parameter  plane  theory  are  presented  here  for  completeness. 
The  following  development  is  included  for  the  purpose  of  showing  the  neces- 
sary and  sufficient  conditions  for  singular  lines. 

Any  set  of  simultaneous  equations  may  be  represented  in  the  following 
matrix  notation: 


[  r\l  Lxi  =lb] 
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where  \~f\  J  is  the  coefficient  matrix  of  order  7n  K  Y) 

L  X  _i  is  the  variable  matrix  of  order  y\   x  I 

EB   is  the  constant  matrix  of  order   yd  X  I 
and  L^  Bj  is  the  augmented  matrix. 
Any  set  of  simultaneous  equations  is  said  to  be  consistent,   (i.e.,  a 
solution  is  guaranteed  to  exist)  when  the  rank  of  the  coefficient  matrix 
and  the  rank  of  the  augmented  matrix  are  equal.   A  unique  solution  may  be 
found  only  when  the  rank  of  the  matrix  equals  the  number  of  unknowns,  Y) 
Otherwise  the  matrix  is  said  to  be  singular.    These  three  statements  are 
very  subtle,  because  their  simplicity  belies  their  importance.   As  an 
illustration  of  their  application,  consider  the  linear  second  order 
equations  which  were  derived  in  the  preceding  section  for  coefficients  of 
the  form: 

-    {  F  (  S J }  =  B.  d.    +■  C,  (3  +  P,  =  0  (2-4-7) 

Im  {  F  (S)}  —  E^oL    -f-  Cx   P  ■+-  D±   —    0  (2.4-8) 

from  which: 

Thus  for  any  characteristic  equation  of  the  form 

^     aK   s*     _Q 

where  the  coefficients  are  linear  combinations  of  alpha  and  beta,  it  is 
possible  to  construct  a  planar  representation  as  a  function  of  the  three 
terms,  X,  V  and  ]?  ,  with  alpha  and  beta  as  the  variable  parameters. 
The  inconsistency  arises  when  the  denominator  determinant   5:  is  zero, 
and  when  at  least  one  of  the  numerator  determinants  is  not  zero.   Since 
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division  of  a  finite  number  by  zero  yields  an  infinite  (or  undefined) 
quotient,  it  is  readily  seen  that  the  two  simultaneous  equations  yield 
two  parallel  lines,  which  by  definition  are  everywhere  equidistant  and 
therefore  never  intersect.   Thus  an  inconsistent  set  of  simultaneous 
equations  has  no  solution.   A  unique  solution  is  obtained  whenever  the  H 
term,  or  denominator  coefficient  is  not  zero.   This  is  the  most  general 
type  of  solution,  denoting  the  intersection  of  two  lines  at  a  finite 
point  in  the  alpha  -  beta  plane.   Should  one  or  both  of  the  numerator 
determinants  X  and  Y   be  zero,  it  merely  places  the  unique  value  on  one 
or  both  axes  respectively.   The  latter  case  of  course  is  just  a  solution 
at  the  origin.   The  third  classification  for  simultaneous  equation  solu- 
tions is  the  singular  case  which  exists  if  and  only  if  all  three  deter- 
minants X,  V  >  anc*  2   are  identically  zero.   Such  a  case  can  exist 
when  the  number  of  unknowns  in  a  set  of  equations  exceeds  its  rank.   Thus 
one  or  more  of  the  equations  is  said  to  be  linearly  dependent.   Geometri- 
cally this  represents  two  collinear  lines  which  by  definition  intersect 
everywhere,  thereby  giving  an  infinite  number  of  solutions.   Extensive 
value  can  be  gained  by  the  use  of  singular  lines  for  design  purposes, 
especially  in  self  adaptive  systems.   It  therefore  becomes  highly  desir- 
able to  be  able  to  recognize  the  existence  of  such  lines  when  depicting  a 
system  characteristic  equation  in  the  parameter  plane.   The  theory  of 
linear  dependence  stipulates  that  two  linearly  dependent  equations  differ 
only  by  a  constant.   Thus  we  may  say  that  for  (2.4-7)  and  (2.4-8)  to  be 
linearly  dependent,  the  following  conditions  must  hold: 

S2  =  k  8, 
Cz  =  fi  C, 
D,.  =  k  D, 
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Because  of  the  exceptionally  involved  calculations  using  the  Chebyshev 
polynomial  method  or  the  binomial  expansion  method,  the  trigometric  solu- 
tions (2.3-5)  and  (2.3-6)  are  amenable  to  singular  line  synthesis  develop- 
ments.  Employing  those  equations,  (2.4-7)  and  (2.4-8)  can  be  broken  down 
into  the  following  equivalences: 


k    £      bi   iV    cos  l&     =     B2=     £     Mai*    sin  16 

k  ^     c{oJp{    tos  IB      =    Cz  =      £-   c-^K     sin  18 

k  D,  =    H<     d;  u/„      cos  10      =      0Z=      £.  ef(-  ^  s^  i  0 


Using  this  criterion,  it  is  possible  to  define  the  required  coefficient 
values  necessary  to  give  singular  line  solutions,  or  to  determine  if 
any  such  combinations  are  indeed  possible  for  a  given  polynomial  order. 
Various  extensions  can  be  made  to  include  all  forms  of  non-linear 
polynomials.   Such  applications  merely  increase  the  number  of  trigono- 
metric equivalence  equations,  and  involve  conic  sections  rather  than 
planar  surfaces,  provided  the  independent  variables  are  limited  to  second 
order.   For  higher  orders,  no  geometrical  interpretation  is  possible;  how- 
ever the  fundamental  theory  remains  unchanged. 

SINGULAR  SOLUTIONS 

While  most  equations  have  explicit  solutions,  certain  differential 
equations  cannot  be  explicitly  defined  by  specifying  the  values  of  the 
arbitrary  constants.   In  Laplace  or  Fourier  transform  theory,  this  corres- 
ponds to  non-parametric  equations.   The  solutions  to  equations  written  in 
this  form  are  called  singular  solutions  [8].   Certain  forthcoming  studies 
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can  be  facilitated  by  expressing  solutions  for  alpha  and  beta  in  non- 
parametric  form.   Up  to  this  point,  the  simultaneous  equations  have  been 
written  as: 

x,  C&\  p';  l)  fatf)    =  ■  o 
X3l(«1,P)  %>  *>«)  =   o 

with  their  associated  solutions  for  alpha  and  beta  written  as 
(X  =  ft   (  I  )  ^n) 

The  following  theory,  therefore,  is  applicable  to  reducing  the  above  equa- 
tions into  the  form: 

Y  («>P>  $)  =-° 

with  associated  solutions 
Ot  =  C  C  (3  ,  S  ) 

and  the  alternate  form 
2  (  <*.,  (3,  LUj  =  O 

with  the  solutions 

CX  =  F  (  f3)  LO„) 
(3  =  F  (<*,^n) 

It  will  now  be  shown  that  for  all  integers   j   and   K  , 
i   MK  aK  ^nK  ^    .  (  §)   =0  (2.4-7) 

whenever  the  relations  (2.1-5)  and  (2.1-6)  can  be  employed.   That  this  is 
the  case  may  be  proved  as  follows:    Given  any  polynomial  of  the  form 

it  is  possible  to  multiply  both  sides  by  U  j-i  C  $)  '. 

i   (V  \  «V   Uj-I  f§)  aK^)  =0  (2.4-8) 

K  =  o 
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Given  also   that 

I      (-l)*    a„    a/     (4  ff)=0 

Multiply  this  by   U  j  (  $ )    to  obtain: 

I     (-1)"  %  «•*"  Uj(jj  U  „.,('$)  =  c 

K-  0 
Subtracting  (2.4-8)  from  this  yields: 

Recalling  (2.2-3)  we  may  rewrite  the  above  in  the  following  form: 

I    (-I}"  *m  "^  E  Ufc-j   (S)J  =  0  (2.4-7) 

K  =  0 

With  the  aid  of  this  relationship  ,  it  is  possible  to  generate  additional 
equations  which  will  enable  (2.1-5)  and  (2.1-6)  to  be  written  in  non- 
parametric  form. 

To  derive  non-parametric  equations  for  constant  zeta  loci,  it  is  neces- 
sary to  eliminate   >  ,   from  the  generalized  simultaneous  equations  which 
then  become 

=  o 

The  necessary  steps  involve  Sylvester's  Dialytic  Method  of  Elimination 
[9]  for  two  polynomial  equations  of  different  order.   An  alternate  approach 
restricted  to  equations  of  equal  order,  can  be  used  with  Bezout's  Method 
of  Elimination  but  is  not  presented  here  [10].   Consider  the  following  two 
polynomial  equations  of  order  rn     and  Ti      where  m    may  not  necessarily 
be  equal  to  71       ,  and  where  the  coefficients  of  omega  are  functions  of 
alpha,  beta,  and  zeta. 
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vn 

K  -  ° 

It  will  be  necessary  to  form  a  square  coefficient  matrix  out  of  the 
combined  equations  to  enable  the  determinant  to  be  taken.   This  deter- 
minant is  then  set  equal  to  zero  in  order  to  obtain  the  non-parametric 
solution  for  all  cases  but  the  trival  one  which  occurs  when  uJ n   simply 
reduces  to  zero.   The  square  matrix  is  of  order  m   X  n   and  is  formed 
by  writing   y  (  uun)      and  ^/ %  ( Wn)  as  follows: 


Y3K)^£  «*»i  UK?)]  =£  £    f^iftfW '«\i  =0    (2.4-8) 


y«K)  =  £   wvL^H)j=i  £   9kKM)u'*  »>»=o   (2.4-9) 


J  ~0 


This  corresponds  to  multiplying  each  term  in  V/(wr.)by  I,  wJr<  )  ^  **• 

Y?-3    ...  Tl-Z    ,  t?-/ 

u>n     LUyi        ^  LUn  respectively. 

In  a  similar  fashion  each  term  in  Y^CuO^)  is  multiplied  by  'j^n; 
^n,"-'  wn  "  respectively.  The  matrix  equivalent  is  shown  for  amplifi- 
cation on  the  next  page. 

Equations  y  (u^)   and  y„(^'Vi  were  not  constrained  to  be  of  equal 
order,  so  that  the  preceding  steps  might  show  the  method  to  be  appli- 
cable for  unequal  cases.   Practical  use  of  this  procedure  for  the  elimina- 
tion of  \a)  y.      in  context  with  the  parameter  plane,  however,  results  in  two 
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o      c      c     £   f,  fz   . 

'   f0   f,   fi   *    '   fn-2    fan  f* 
O   f.   f,   ^    '    '   fm-3  «H  ** 


f.   ft   4 


>'(u,) 


0 


o  9.  5.  1 


0 
0    0 

c 


0  3«  £ 
o  o  o  s.  &  9» 

3,   9,  '  U     c 


b 

C 

. 


a/* 

J  i 


WJ 


=  0 

(2.4-10) 


equations  of  equal  order  when  the  following  substitutions  are  made  for  fK 

and  q  K  : 

The  non-parametric  equation  for  constant  LUT)    is  derived  by  tactics 
analogous  to  those  used  to  generate  the  constant  zeta  equation.   Again, 
a  set  of  equations  is  obtained  for  the  purpose  of  taking  the  determinant 
of  the  coefficient  matrix  which,  by  definition,  must  be  a  square  matrix 
if  the  determinant  is  to  exist.   Since  the  variable  matrix  \A^(^)    is 
only  of  the  first  power,  it  is  possible  to  generate   71  -  i    additional 
equations  by  use  of  (2.4-7)  which  was  previously  derived.   Letting  j 
range  from  zero  to   7)  -  I    in  (2.4-7),  the  following  set  of  simultaneous 
equations  is  obtained: 
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£   8r   *■*   «V   UK(§)  =  0 


(2.4-11) 


£  ('l)K      *k  <V  UK-fr-,)(J>  = 


1   fl)K   **  "n*  ^..^^(0=0 


*-o 


Writing  these  equations  in  expanded  form: 


r-i)'a0  ">„°U -(„.,;  (J) ■  ■+-  f-0  a,  w„  U-cn-Dd)   +   -•• 
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These  equations  may  be  simplified  by  substituting  the  following 
relations: 

u0  (v  =  c 


Thus  we  may  write: 

o 


(-0  a,  tu„  u,(f) 
O 


Which  is  seen  to  be  a  set  of  n-1  equations  with  n-1  unknowns.  In 
matrix  notation  these  become: 


+  (-0  tCLxwn%  +-•  H)W 


-\)   a,  ujn 


l3   I 


0 


(2.4-12) 


JL   j 
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The  above  developments  have  shown  how  a  set  of  two  n   order  linear- 
ly independent  equations  with  four  variables  can  be  reduced  to  one  equa- 
tion of  order  n-1  with  only  three  variables.   This  is  accomplished  by 
solving  the  coefficient  matrices  of  (2.4-10)  and  (2.4-12)  which  elimin- 
ate zeta  and  omega  respectively.   The  application  of  these  methods  will  be 
indicated  in  the  following  article  where  the  mathematical  concepts  of  the 
frequency  response  curves  are  presented. 

2.5  Frequency  Response  Applications  of  Parameter  Plane  Theory 

Frequency  response  curves  can  be  portrayed  in  many  ways.   Some  of  the 
more  useful  ones  will  be  considered  in  sections  3  and  4.   The  various 
representations  of  frequency  response  curves  can  be  derived  from  the 
following  fundamental  theory  relating  magnitude  to  frequency.   Consider 
again  the  transfer  function  of  a  mathematical  model  which  represents  the 
system  to  be  studied.   This  can  be  reduced  to  a  ratio  of  two  polynomials 

by  partitioning,  or  other  suitable  means  to  arrive  at: 

n 
-r,0_  ao+*.s+az*fc  +  fi3*3+"  qnsB  jj0  at  s  tf(s) 

I    O;  — — = (2.5-1) 

bc  +b,  S  +  bzs*  4-fc,53+"«  bmS™  3     ,      J  D(S) 


j=0 


To  find  the  absolute  magnitude  of  this  transfer  function,  the  following 

is  noted: 

Z 


MZ  =      TCJw)  =  [TCJw)J  \TC-juj)\ 


(Although  the  notation  standards  of  section  2  clearly  state  the  distinc- 
tion between  UJ  and  'JJ  ^       ,  it  is  important  to  note  that  natural 
angular  frequency  is  no  longer  being  used.)   The  following  steps  will 
illustrate  this  transition: 
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ae   f   aJuj    t  a^Jco)    -h'a3  (juj)3  +.  . .  .  an  (jt*j)  ' 

TCJu-0  = : 

Rewriting  this  to  separate  the  real  and  imaginary  parts,  we  have: 


(ju,-)  = 


T(jou)  = 


br  +  ba.(jou)1  -+-  bM  (ju;)*  +  m*  •#-  bi  j*tf  +•  h3  (jw)  +.». 

[a,  -a^  u^fa^u'4   ••  j  +  j  £  a,  w;  -a2  u;?  +  af  fd««j 


In  abbreviated  notation,  this  becomes: 

T-fji  I  =       Meven    (*U)) +■  <V  odd  (  J  w) 

Proceding  in  a  similar  manner  we  find  the  following: 

r  c-jw)  3= 

[b0  -  b^v*  +  *>«,*"*-  •  •  ■  ]  -j  f.b|W  -^wU  bro>r-«D 

Employing  abbreviated  notation  again,  we  may  write: 
Neve*    (   $  <»)     -Nc<tcH(]w) 

T  (-jw)  =  ~ — — 

Thus   the   square  of   the   absolute  magnitude   is   found  to  be: 

[tieven  (}u>)    '+■  NoddCJu))'1l  IV  even  Ci"^  -NoJJ  G">)~] 


M*  = 


\_tieve*  (JWJ   +■  DodJ  6:"f[  [  De^e«  CM  ~  Dodd  CiuJ)] 


% 


M     = ±, i.  (2.5-2) 


5"   -,2 


M  = 
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This  can  be  written  in  the  following  summation  notation: 


Z-n       K. 
±        i'\)%     c 
K  even 


£  f-l) 

K  odd 


ILL1 
Z 


Xrn 


±      ir\)   *■  dtvj- 


even 


2-m-i  J_~J  j 

£     (-0  %  d.  w 


71 

H=0 

k 

(-1) 

C>tK 

uj 

ZK 

Jt  =  0 

H)i 

^ 

6t; 

11 

X   - 

2 
LV 

9     - 

-.  c, 

k 

Substituting  these  new  variables  above,  a  more  concise  notation  may 
be  formalized  for  easier  illustration  of  pertinent  relationships. 


M  =  -* 


£.    (-0     eK  * 


i-o  * 


(2.5-3) 


In  general,  the  coefficients  of  (2.5-1)  are  non-linear  combina- 
tions of  at  least  two  variable  parameters,  designated  alpha  and  beta. 
The  linear  combinations  were  discussed  in  article  2.3  and  the  associated 
form  of  any  polynomial  coefficient  was  found  to  be 

For  the  non-linear  case,  the  simplest  combination  was  seen  to  be  the 
alpha-beta  product  combination.   This  subject  was  treated  in  article  2.4 
where  the  resultant  coefficients  of  a  given  polynomial  were  broken  down  in 
the  following  manner: 
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Theoretically,  there  is  no  limit  to  the  degree  of  difficulty  which  exists 
in  the  non-linear  case.   A  few  examples  of  the  complexity  involved  in 
non-linear  parameter  combinations  are  presented  here  for  completeness. 


aK  =    bK  *-    +-    ck   P     +-    <*K  ocf   +-  eK     +  f K  oc     ■+-  ^  fi 


Clearly,  it  can  be  seen  that  such  combinations  are  exceedingly  unwieldy. 
Attempts  to  handle  all  but  the  simplest  arrangements  are  extremely  labori- 
ous.  Computer  programming  is  still  in  the  embrionic  stage  for  parameter 
plane  applications.   Attempts  to  utilize  this  tool  have  been  extended 
only  to  the  alpha-beta  non-linear  case.   This  restriction  is  acute  when 
considering  the  frequency  response  curves.   A  cursory  review  of  the 
mathematical  developments  at  the  beginning  of  this  article  will  show  that 
the  polynomial  coefficients  of  (2.5-3)  are  obtained  by  squaring  the  poly- 
nomial coefficients  of  (2.5-1).   Thus,  even  the  simple  linear  case  of 
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article  2.3  results  in  non-linear  combinations  of  second  order  for  fre- 
quency response  studies.   To  handle  this  problem,  it  becomes  necessary 
to  depart  from  the  procedures  which  were  employed  in  previous  develop- 
ments.  By  assigning  a  value  to  either  alpha  or  beta,  the  task  is  reduced 
to  merely  a  quadratic  formula  solution  for  the  remaining  unknown  variable. 
This  can  be  illustrated  in  the  following  manner:   Given  a  non-linear  co- 
efficient combination  of  the  alpha  -  beta  product  form: 
aK=  bKot-t-CK(3  +    c(Ko((3  +  eh 
Assign  a  value  to  alpha,  and  denote  the  value  as  A: 

This  may  be  rearranged  for  clarity, 


a 


[  M  +  enl     +  '[>*  +  ***]  P 


K  —  L    ^H  ''     '     "K-*       '   L"K  '  "K   "J  r  (2.5-4) 


where  the  first  term  on  the  right  hand  side  denotes  a  constant  value,  and 
the  second  bracketed  term  represents  a  coeficient  of  the  independent  vari- 
able beta.   In  a  similar  manner,  a  fixed  value  may  be  assigned  to  beta. 
For  illustration  purposes,  let  it  be  denoted  as  B.   This  gives: 

ah  =  iH  a  -r-   cK  B    •+-   clk'ctB    +■  eH 

which  is  likewise  rewritten  as: 

aK  =  [  cK  b  +  ej  +  [^  +  cIhbJol  (2.5-5) 

Again,  the  first  bracketed  term  represents  a  constant  while  the  second 
denotes  the  coefficient  of  the  independent  variable  alpha.   The  computer 
program  for  frequency  response  curves,  program  PARAMS,  is  based  upon  this 
solution  procedure.   Consider  first  the  case  where  magnitude  and  omega  are 
specified.   By  assigning  incremental  values  to  either  alpha  or  beta,  the 
corresponding  solution  may  be  obtained  for  the  other.   This  is  the  method 
employed  in  the  graphical  representations  which  have  both  alpha  and  beta 
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on  the  axes.   For  the  graphs  where  magnitude  is  an  axis  variable,  it  is 
incremented  while  a  corresponding  value  of  beta  is  obtained  as  a  function 
of  alpha  and  omega.   The  same  philosophy  is  utilized  when  omega  is  an 
axis  variable.   For  the  case  where  both  magnitude  and  omega  are  on  the 
axes,  the  values  of  alpha  and  beta  are  pre-selected  and  the  familiar  Bode 
plot  results.   All  of  these  graphical  arrangements  are  treated  with  con- 
siderable detail  in  section  3.   Their  introduction  at  this  point  was 
simply  to  point  out  that  the  solutions  are  based  on  (2.5-3)  and  either 
(2.5-4)  or  (2.5-5).   Since  (2.5-3)  is  the  fundamental  gain  versus  fre- 
quency equation,  it  can  be  anticipated  that  all  frequency  response  para- 
meter plots  will  contain  the  same  basic  information  as  a  series  of  con- 
ventional Bode  frequency  response  curves.   When  two  variable  parameters 
alpha  and  beta  are  utilized  together  with  magnitude  and  omega,  the  problem 
becomes  one  of  four  dimensional  space.  When  still  another  polynomial  co- 
efficient variable  parameter  is  considered,  an  additional  dimension  is 
added.   In  spite  of  this,  it  is  still  possible  to  derive  one  set  of  para- 
meter curves  from  another,  since  each  is  just  a  different  projection  from 
the  same  multi -dimensional  model.   The  value  of  each  separate  representa- 
tion lies  in  the  ability  to  depict  certain  important  characteristics  which 
would  not  be  obvious  otherwise.   Additionally,  certain  configurations 
facilitate  the  study  of  a  dominant  parameter,  or  the  effect  of  certain 
parameters  on  different  performance  requirements.   Certain  mathematical 
considerations  must  be  pointed  out  at  this  time  concerning  the  inter- 
dependence of  these  graphical  plots.   Clearly,  the  construction  of  PARAM-7, 
the  Bode  plot,  offers  no  problem.   In  this  case,  the  desired  values  of 
alpha,  beta,  and  omega  are  inserted  into  (2.5-3)  and  the  corresponding 
magnitude  squared  value  is  then  computed.   For  the  remaining  programs 
(PARAM-1)  through  (PARAM-6) ,  a  rearrangement  of  (2.5-3)  is  necessary  since 
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magnitude  is  no  longer  the  unknown  variable.   As  an  illustration,  con- 
sider (2.5-4)  where  alpha  has  been  specified.   Beta  then  becomes  the  un- 
known value  to  be  solved.   For  added  clarity  it  is  better  to  revert  back 
to  (2.5-2)  in  order  to  point  out  why  a  second  order  system  is  always  en- 
countered in  the  linear, case,  as  well  as  the  non-linear  alpha  -  beta  pro- 
duct case.   Since  beta  is  the  unknown  variable,  pre-selected  values  of 
magnitude,  omega,  and  alpha  must  be  presumed.   Performing  these  substitu- 
tions, (2.5-2)  may  be  written  as 

where,  using  the  notation  of  (2.5-4) 

f K  =     L  bo  ft  *■  ej  -     [  t,/?  f  ez~]  uj2-   +  jTb^fl  ±eHl^H-  ■  •  • 

%k  =     L  c°  +  ^Pil    -    lei  -h    dzfl]  uj2-    +.    [cti    +JH  nj  u)q-   -  ■> 

>%  —    [  fcfl0  +  e0]-  [bzfl  +  ez  J  u>*   +  ly   ft  +ehr\v,*-"i 


[c0  +    ddfi]     -   [c%+  <jxfll  tuZ    f   [cH   +dnfi~[ 


uo1*. 


Therefore,  one  can  easily  see  that  by  squaring  the  right  hand  bracketed 
values  and  combining  like  terms, 

i_     $NBO      +    c,N  &  I    (Z      +  snBZ    £ 

where  the  terms  SNBO,  SNBI,  SNB2,  SDBO,  SDB1,  and  SDB2  are  so  chosen  be- 
cause they  are  the  variables  used  in  the  computer  program  to  be  dis- 
cussed later.   Since  the  numerical  value  of  M   has  been  specified,  this 
equation  can  be  cross-multiplied.   Transposing  the  right  hand  side  of  the 
resultant  equation,  one  can  write: 
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SNBC   ~M%SDBO    +  I  $  N  B  I   -M2$Di3ll£    +  l$N6Z  '  M*SD82]fi>Z=  0 

This  is  clearly  a  quadratic  in  beta,  and  employment  of  the  quadratic  equa- 
tion yields  the  desired  results.   Imaginary  numerical  values  have  no 
physical  significance  for  graphical  plotting  and  are  therefore  discarded. 
However  both  the  plus  and  minus  signs  in  front  of  the  square  root  term 
of  the  quadratic  formula  must  be  considered.   Thus  it  is  possible  to 
have  two  values  of  beta  for  a  specified  combination  of  magnitude,  omega, 
and  alpha.   This  is  an  important  consideration  in  graphical  plotting 
techniques  which  can  lead  to  spurious  lines  if  not  properly  accounted 
for.   Effectively  the  problem  requires  calculating  all  the  prescribed 
graph  points  using  first  one  sign  of  the  square  root  term,  and  plotting 
them.   Then  the  computational  process  is  repeated  with  the  alternate 
sign  of  the  square  root  term.   This  illustration  has  shown  how  a  solu- 
tion would  be  obtained  for  beta  when  alpha,  magnitude,  and  omega  are 
known.   Precisely  the  same  procedure  would  be  followed  to  obtain  a  solu- 
tion for  alpha  when  the  values  of  beta,  magnitude,  and  omega  are  given; 
only  (2.5-5)  would  be  used  in  lieu  of  (2.5-4) 

RESONANT  PEAK  MAGNITUDE 

Performance  specifications  of  certain  systems  include  stringent  re- 
quirements on  the  presence  of  a  stipulated  resonant  peak  magnitude.   Con- 
versely, other  design  studies  are  concerned  with  the  guaranteed  absence 
of  a  resonant  peak  magnitude  above  a  pre-determined  level.   Either  of 
these  requirements  results  in  the  study  of  frequency  response  curves  with 
particular  emphasis  on  the  peak  magnitude  and  the  frequency  at  which  it 
occurs.   To  facilitate  such  studies,  it  would  be  advantageous  to  construct 
resonant  peak  magnitude  curves  for  the  system  as  a  function  of  the 
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polynomial  variable  coefficients,  alpha  and  beta.   This  involves  the 
utilization  of  the  non-parametric  theory  presented  at  the  end  of  article 
2.4. 

Normally,  one  thinks  of  only  one  maximum  point  in  the  entire  fre- 
quency spectrum  because  of  dominant  root  characteristics  of  low  order 
polynomials.   Certain  high  order  systems,  however,  can  be  represented  by 
characteristic  equations  which  have  more  than  one  pair  of  dominant  com- 
plex roots.   Such  arrangements  will  yield  more  than  one  maximum  point  on 
the  frequency  response  curve.   Only  one  of  these  points  can  be  consider- 
ed as  the  global  maximum,  while  all  other  peaks  must  be  relegated  to  the 
role  of  local  maxima.   To  avoid  ambiguities,  the  notation  Mp(V)will  be 
restricted  to  the  global  maximum  point  only.    The  corresponding  value 
of  omega  will  be  denoted  by  (J„  .   Since  it  is  desirable  to  refer  to 
local  maximum  points  as  well,  the  term   M  ^  (io)   shall  be  construed  to  be 
applicable  to  all  resonance  points,  both  local  and  global.   The  respec- 
tive angular  frequencies  at  which  these  occur  will  be  labelled  as  id      . 
Thus  it  is  possible  to  have  any  number  of  MY(io)   and  LJ  r  values  for  a 
given  transfer  function.   Only  one  of  these  can  be  the  /V|rtMand  UJp   , 
however. 

Recalling  (2.5-3)  where  e         and  fR   are  both  functions  of  the  vari 
able  parameters  alpha  and  beta,  it  is  possible  to  state  this  equation  in 
the  following  functional  notation: 

M*M  -     '        ■   ' »• (2.5-6) 


fj.  (*•>    P;  W) 
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Assuming  at  least  one  pair  of  dominant  complex  conjugate  roots,  it  is 
possible  to  find  the  following  relation  from  (2.5-6): 

Mr%)   -   ft(dL,P,M)         ^_     f,  U,(V<.)  (2*5*7) 


-    rf.,>,  *) 


i/J-'XV 


K  /^    =  Al^A^I  _i£M5?  (2.5-8) 

flKftw)    I  fi  Cot,  j»/4) 

Fundamental  theory  of  equations  dictates  the  following  correlation  be- 
tween (2.5-6),  (2.5-7)  and  (2.5-8): 

1.   A  relation  in  the  form  of  (2.5-7)  exists  whenever  the 
slope  of  (2.5-6)  is  equal  to  zero  and  its  radius  of  curvature  is  nega- 
tive.  Expressing  this  statement  in  mathematical  terms: 

±hll^l         =    C  (2.5-9) 

/r  m  2  (wf\  i 

=   0  (2.5-10) 


2.      Only   that  value  of     Ujt      in   (2.5-7)   which  globally  maxi- 
mizes   (2.5-6)    satisfies    (2.5-8).      Thus   it  is   possible   to  state: 


M     Z(i>j)      =        TOa.xj.nn4m  5    /vf*"^) 


(2.5-11) 

J 


The  preceding  theory  is  adequate  in  that  it  gives  the  necessary  and  suf- 
ficient conditions  for  the  existence  of   H\ «  (uv)»      Satisfying  the  above 
criteria  is  fairly  straight  forward  whenever  the  values  of  alpha  and  beta 
have  been  selected  or  prescribed.   The  dilemma  arises  when  seeking  the 
algebraic  relationship  for   fA  p  \UJ)     as  a  function  of  alpha  or  beta. 
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Although  three  equations  are  available,  the  inequality  in  the  second  de- 
rivative criterion  of  (2.5-10)  imposes  a  severe  limitation  on  the  solu- 
tion technique.   The  most  straight  forward  approach  would  be  to  utilize 
(2.5-7).   From  these  it  is  possible  to  obtain  (2.5-11).   Such  a  method 
of  attack  is  only  useful  when  both  alpha  and  beta  are  given.   This  would 
require  a  two  dimensional  scan  for  all  possible  values  of  alpha  and  beta. 
Various  computerized  programs  exist  in  optimal  control  theory  for  the 
determination  of  (2.5-7).   These,  however,  are  based  on  apriori  knowledge 
or  random  search  techniques  using  Lagrange  multipliers.   Not  only  are 
these  methods  time  consuming,  they  are  not  conducive  to  highly  sensitive 
frequency  response  curves  where  resonant  peaks  are  in  close  proximity  to 
one  another.   An  alternative  solution  employs  the  non-parametric  theory 
which  results  in  (2.4-10).   This  procedure  is  applicable  only  when  the 
two  equations  involved  can  be  set  equal  to  zero.   Therefore  it  is  neces- 
sary to  investigate  the  locus  of  inflection  points  of  the  magnitude  versus 
frequency  curve.   From  this  information,  it  is  possible  to  determine  how 
many  local  maximum  points  exist.   It  can  be  shown  that  for  Tl    inflection 
points,  there  can  be  no  more  than  T)   maxima.   In  practice,  however,  only 
Yl  -  I  will  be  found  for  most  systems.   This  procedure  has  the  important 
restriction  of  eliminating  omega.   Thus  it  is  impossible  to  employ  (2.5-6) 
in  order  to  find  the  value  of   M  [u))   at  each  inflection  point.   At  best, 
this  result  makes  it  possible  to  construct  the  locus  of  inflection  points 
in  the  parameter  plane  when  alpha  and  beta  are  the  axis  variables.   It 
then  becomes  necessary  to  employ  some  other  means  of  determining  which  side 
of  this  locus  represents  the  minima,  or  maxima.   Full  utilization  of  the 
non-parametric  theory  can  produce  this  desired  result.   Inspection  of  the 
developments  leading  to  the  coefficient  matrix  of  (2.4-12)  will  show  that 
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the  combined  application  of  (2.4-10)  and  (2.4-12)  provides  the  necessary 
information  for  (2.5-6)  as  follows: 

1.  The  results  of  (2.4-10)  yield  a  solution  of  the  form: 

X  =   f7  (d,   (3  )  (2.5-12) 

2.  The  results  of  (2.4-12)  provide  the  following  alternate 
solutions: 

V  =  f    (   *,  oj)  (2.5-13) 


Vz  =   f  ^  (  f,  w)  (2.5-14) 

depending  on  which  coefficient  variable  is  eliminated.   From  these  re- 
lations, it  is  possible  to  employ  (2.5-6)  directly  by  use  of  (2.5-12) 
and  (2.5-13),  or  (2.5-12)  and  (2.5-14).   In  this  manner,  it  would  be 
possible  to  employ  the  optimal  control  theory  previously  mentioned. 
However  the  more  direct  approach  would  be  the  utilization  of  the  gradi- 
ent technique  which  is  presently  employed  in  many  optimal  control  pro- 
blems.  Alternately,  either  of  these  pairs  may  be  solved  simultaneously 
to  provide  a  locus  of  inflection  points  as  a  function  of  the  two  re- 
maining variables,  alpha  and  omega,  or  beta  and  omega.   This  latter  ap- 
proach would  be  most  desirable  for  frequency  response  parameter  plots 
employing  these  variables  as  the  graphical  axes. 

The  preceding  theory  has  shown  several  ways  of  handling  the  peak 
resonant  magnitude  problem.  Inherent  in  each  of  these,  however,  is  a 
vast  array  of  mathematical  manipulations  which  is  limited  by  solution 
techniques  presently  available.   The  following  development  will  indicate 


57 


this  more  clearly.   Consider  again  the  fundamental  relationship  of  (2.5-3) 
upon  which  these  derivations  are  based.   The  following  result  is  obtained 
for  (2.5-9): 
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Jt  =  0 


(2.5-15) 


Let  q  =  numerator  coefficient  of  (2.5-15)  with  the  sign  embedded,  and 
let  Y\  n    =  denominator  coefficient.   Thus  (2.5-15)  can  be  rewritten  as: 


d    M     (x.) 
J  * 


x--o  l         xi 


The  equivalent  form  of  (2.5-10)  then  becomes: 


c^M  2(*) 


d'y 


3L    (-4)      h. 


I  'LH'tu'"! 


/=o  , 


(2.5-16) 


w     / 


2/-/ 


L  ho    h*m&')'l'»M*~~] 


Thus  from  (2.5-16)  it  is  noted  that  the  order  of  the  denominator  poly- 
nomial is  raised  to  a  power  of  four  times  the  order  of  the  denominator  in 
the  original  transfer  function.   In  like  manner,  it  is  evident  that  the 
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numerator  polynomial  of  (2.5-16)  is  also  of  very  high  order.   Therefore, 
all  but  the  simplest  of  transfer  functions  would  require  computer  assist- 
ance to  solve  the  determinants  of  (2.5-15)  and  (2.5-16).  Although  many 
programs  exist  for  explicit  determinant  solutions,  these  cannot  be  used 
when  one  of  the  terms  is  an  unknown  variable.   Until  a  generalized  com- 
puter program  is  developed  to  handle  these  cases,  the  construction  of  peak 
resonant  curves  can  only  ba  achieved  by  a  laborious  calculation  of  second- 
ary curves.   Such  procedures  do  exist,  whereby  tangent  lines  are  drawn 
from  a  series  of  curves.   To  be  highly  accurate,  however,  an  excessive 
number  of  these  are  required  for  sensitive  systems. 

PHASE  CALCULATIONS 

The  preceding  discussions  in  this  article  have  been  limited  to  magni- 
tude studies  in  the  frequency  domain.   An  equally  large  field  of  interest 
is  that  of  phase  angle  characteristics.   The  utility  of  the  magnitude 
versus  frequency  plot  lies  not  only  in  analyzing  or  predicting  the  fre- 
quency response  of  a  system,  but  also  in  the  ability  to  check  its  stabil- 
ity as  well.   The  correlation  between  magnitude  and  phase  curves  provides 
this  verification  by  means  of  phase  margin  and  gain  margin  calculations. 
It  therefore  becomes  desirable  to  derive  the  phase  calculations  as  a  func- 
tion of  the  three  variables:   alpha,  beta  and  omega. 

Initial  developments  in  this  article  showed  how  (2.5-1)  could  be  alter- 
ed to  produce 


A/  even  (j iJj)       +      N  odd    (  Jw) 

ja  iv)   = — — 

|     even   (J^)     +     D  odd    (juj) 
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Considering  the  numerator  and  denominator  separately,  either  can  be  re- 
garded as  the  real  and  imaginary  components  of  a  vector  with  a  given  magni 
tude  and  angle.   Magnitude  calculations  have  already  been  developed  and 
are  of  no  interest  here.   Phase  angle  calculations,  however,  are  known 
to  involve  the  arc-tangent  of  the  vector  components  as  follows: 


P '  li  as  e  ft  nql  e    = 


fan 


[TCjuj)} 


Re  {  T(jw)] 


(2.5-17) 


Two  solution  procedures  are  therefore  possible.   One  is  to  solve  independ- 
ently for  the  angle  of  the  numerator  vector  and  subtract  the  denominator 
vector  angle  from  it.   Thus: 


fa.  n 
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(2.5-18) 


The  alternate  approach  involves  only  one  trigonometric  conversion.   Clear- 
ing the  imaginary  term  from  the  denominator  of  T~  L]u>)   , 
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Therefore,  substituting  this  into  (2.5-17), 


6    =    taw 


I   Nodd  (jw)  Oe^en(j^)  "  Never)  ('}">)  Oodi  (^ 
|__Neven(Juj)    Pei/en(jw)-  Mode/  (  ju>)  Oodd0^_ 


(2.5-19) 


Either  (2.5-18)  or  (2.5-19)  may  be  used  for  the  prescribed  phase  calcu- 
lation.  For  the  Bode  plot  (PARAM-7) ,  both  have  the  same  degree  of  diffi- 
culty in  computer  applications,  although  (2.5-18)  is  preferred  on  long 
hand  calculations.   Extensions  to  phase  curves  in  the  remaining  para- 
meter plane  frequency  response  graphs  (PRAAM- 1  through  PARAM-6)  have  not 
been  attempted.   These  would  be  possible  without  the  use  of  the  quadratic 
formula  if  (2.5-18)  were  used.   Selection  of  (2.5-19)  requires  the  usage 
of  the  quadratic  formula  with  both  its  square  root  coefficients  as  il- 
lustrated in  the  magnitude  portion  of  this  article. 
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FIGURE  2-2 
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FIGURE  2-4 
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FIGURE  2-5 
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3.   Computerized  Graphical  Construction  Techniques  For  Frequency 
Response  Studies. 

Graphical  analysis  has  always  held  a  virtuous  distinction  in  all 
facets  of  industry  and  business.   Since  the  Second  World  War,  however, 
the  advances  of  technology  promoted  the  engineering  profession  beyond 
the  existent  graphical  methods  which  were  derived  for  the  single  purpose 
of  varying  only  one  parameter  at  a  time.   In  spite  of  this,  the  ancient 
adage  "one  picture  is  worth  a  thousand  words"  has  withstood  the  test  of 
time,  as  witnessed  by  many  practicing  engineers  today.  Many  pictorial 
representations  of  single  variable  systems  are  familiar;  however  with  the 
advent  of  multiple  parameter  variations,  graphical  representation  in  two 
dimensions  becomes  more  complex.   Considerable  experience  and  moderate 
imagination  are  necessary  to  visualize  three  dimensional  configurations 
when  presented  in  two  dimensional  planes,  in  order  to  adequately  make  an 
intelligent  intrepretation  for  meaningful  results. 

The  recent  developments  employing  algebraic  methods  for  variation 
of  two  parameters  have  resulted  in  the  successful  use  of  digital  com- 
puter facilities  to  perform  the  required  manipulations  for  graphical 
representation  of  mathematical  equations.   The  preceding  section  was 
devoted  to  a  fundamental  development  of  the  underlying  mathematical  theory 
employed  in  two  parameter  design  studies.   This  section  will  include  a 
broad  presentation  of  certain  graphical  methods  which  are  most  likely  to 
be  employed  in  two  parameter  electrical  control  problems,  with  primary 
application  to  frequency  response  studies.   This  section  is  also  con- 
cerned with  a  detailed  inspection  and  description  of  the  computer  pro- 
grams utilized  in  formulating  those  graphs.   The  discussion  is  not  direct- 
ed toward  programming,  but  rather  an  effort  is  made  to  clarify  certain 
computer  oriented  problems  so  that  engineers  who  desire  to  utilize  any 
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program  contained  herein  may  do  so  with  minimum  difficulty  and  yet  gain 
maximum  benefit  from  its  comprehensive  construction.   Illustrative  ex- 
amples of  each  of  these  programs  appear  in  sections  4  and  5  where  several 
case  studies  are  investigated  in  depth. 

3.1   Selection  of  Graphical  Representations. 

For  any  given  set  of  six  variables,  the  study  of  their  respective 
inter-relationship  by  means  of  two  dimensional  projections  results  in  a 
myriad  of  graphical  portrayals.   The  magnitude  of  the  problem  can  best 
be  seen  in  the  following  development:   Assume  that  an  arbitrary  number 
of  variable  parameters  are  to  be  considered.   If   T)    is  the  number 
design  variables,  it  is  possible  to  show  that  the  number  of  possible 
graphical  representations  is  71  -factorial  (  7?  .  ).   All  graphs  reveal 
essentially  the  same  information  content  when  one  axis  variable  is  switch 
ed  from  the  ordinate  to  the  abcissa,  and  vica  versa.   Therefore,  redun- 
dancy is  eliminated  by  allowing  the  abcissa  and  ordinate  variables  to  be 
interchangeable.   This  reduces  the  number  of  combinations  by  a  factor  of 
two,  thereby  making  the  revised  number  of  graphical  representations  equal 
to  %  (  7)  [   ).   Although  a  considerable  reduction  in  the  size  of  the  ar- 
ray has  been  effected  by  interchanging  the  axis  variables,  it  is  highly 
desirable  to  reduce  this  size  further,  since  a  six  variable  array  means 
360  graphs,  a  five  variable  array  requires  60  representations,  and  so 
forth.   To  further  reduce  this  compilation,  it  is  often  possible  to  inter- 
change any  two  parameters  alpha  (  °^  )  and  beta  (  f3   ).   Such  would  be 
the  case  in  the  parameter  plane  technique  where  they  are,  by  definition, 
arbitrary  variables  of  the  characteristic  equation.   Even  this  reduction 
leaves  a  formidable  task  of  evaluating  over  180  graphs  for  a  system  in- 
corporating six  variables,  and  any  system  with  five  variables  requires 
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an  array  of  33  different  graphical  representations.   For  this  reason 
alone,  the  present  state  of  the  art  is  limited  to  at  most  four  variables 
for  practical  design  techniques  employing  graphical  displays.   Any  system 
with  four  variables  requires  basically  seven  different  graphical  presenta- 
tions for  comprehensive  evaluation.   Seldom  are  all  of  these  required  for 
a  satisfactory  design,  although  barring  the  elements  of  time  and  expense 
they  would  most  likely  be  necessary  for  optimal  results.   Effectively  the 
problem  at  hand  involves  the  inspection  of  four  independent  variables. 
This  can  be  regarded  as  a  four  dimensional  space  which  is  not  physically 
realizable.   By  assigning  a  constant  value  to  one  of  these  four  variables, 
the  problem  is  reduced  to  a  three  dimensional  sub-space  representation 
which  does  have  a  physical  shape.   It's  meaning,  however  is  restricted  to 
that  particular  value  of  the  fourth  variable.   To  present  this  three  di- 
mensional sub-space  on  a  two  dimensional  plane  requires  that  a  large 
number  of  curves  be  projected  onto  this  plane,  where  each  curve  represents 
a  different  value  of  the  third  variable.   The  following  discussion  is 
based  on  this  concept,  and  the  resultant  graphs  should  be  interpreted  ac- 
cordingly. 

Most  specifications  for  the  dynamic  performance  of  linear  systems  are 
written  in  terms  of  the  time  domain  and  usually  include  such  terms  as  peak 
overshoot,  transient  duration,  and  steady  state  error.   It  therefore  be- 
comes necessary  for  the  control  engineer  to  bridge  the  time  domain  with 
the  more  utilitarian  frequency  domain.   Correlation  is  accomplished  by 
means  of  the  Fourier  integral  when  rigorous  methods  are  called  for.   This 
treatment  is  avoided  in  practice  however,  because  of  the  inherent  arduous 
task  involved  in  all  but  the  simplest  of  integrations.   It  can  be  shown 
that  by  suitable  approximation,  many  linear  problems  fall  into  the  cate- 
gory of  second  order  systems.   For  engineers  not  experienced  in  partitioning 
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methods  and  root  locus  techniques,  this  statement  would  seem  to  be  easily 
contested;  however,  such  an  analogy  actually  exists  in  most  cases.   For 
systems  which  cannot  be  resolved  into  an  equivalent  second  order  approxi- 
mation, the  absence  of  dominant  complex  roots  is  evidenced  by  the  shape  of 
the  closed  loop  frequency  response  curve.   Since  considerable  information 
exists  in  present  day  literature  on  the  generalized  second  order  system, 
much  time  and  labor  is  saved  by  manipulating  a  given  problem  into  this 
form  whenever  possible,  and  extracting  the  equivalent  damping  and  band- 
width specifications  from  pre-computed  curves  and  tables. 

A  cursory  review  of  the  terms  used  in  control  engineering  applica- 
tions will  reveal  that  some  of  the  more  common  design  criteria  are  steady 
state  magnitude  (  M   ),  phase  angle  designated  theta  (  S      ),  natural 
angular  frequency  designated  omega  (  uU ^       ),  and  damping  ratio  designa- 
ted zeta  (   ?   ).   Some  problems  are  more  amenable  to  the  rectangular 
coordinates  of  the  complex  frequency  plane.   These  are  designated  as  the 
damping  constant  sigma  (  cr  )  and  the  damped  angular  frequency  omega 
(  uJ     ) .   It  is  again  emphasized  that  two  forms  of  angular  frequency  are 
designated  by  the  word  omega.   Much  of  the  existing  literature  uses  this 
term  very  loosely,  and  extreme  care  must  be  exercised  not  to  mistake  un- 
damped, or  natural  angular  frequency  (  kJTi  )  for  damped  angular  frequency 
(  U)     ).   The  relation  between  these  two  variables  is  found  in  article  2.1, 
and  it  should  be  noted  that  these  terms  are  synonomous  only  when  the  damp- 
ing ratio  (  f   )  equals  zero. 

Section  2  developed  the  notational  concept  of  alpha  (  oi    )  and  beta 
(   (3  )  as  representing  two  arbitrary  variables  embedded  within  the  co- 
efficients of  the  characteristic  equation.   The  preceding  paragraph  has 
shown  that  eight  variables  are  frequently  encountered  in  the  resultant 
electrical  engineering  control  problems.   As  previously  stressed,  the 
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simultaneous  usage  of  all  of  these  is  not  feasible.   Fortunately,  however, 
one  may  break  the  control  problem  down  into  two  basic  categories:   stabil- 
ity and  frequency  response.   These  two  facets  are  interrelated,  however, 
and  for  that  reason  some  attention  must  be  given  to  system  stability 
when  performing  frequency  response  studies. 

The  parameter  plane  stability  problem  may  effectively  be  considered 
as  the  study  of  the  relation  between  the  characteristic  equation  poly- 
nomial and  the  stability  performance  variables  £   and  iOn    ,  or  <T      and 

U)      .   When  the  two  parameters  alpha  and  beta  are  embedded  within  the 
characteristic  equation,  the  scope  of  the  problem  increases  to  a  study  of 
at  least  four  variables.   Thus  one  might  be  inclined  to  consider  seven 
possible  graphical  arrangements.   Unlike  the  frequency  response  studies, 
however,  the  stability  problem  has  been  attacked  from  a  slightly  differ- 
ent viewpoint.   Only  alpha  and  beta  have  been  plotted  along  the  axes,  and 
no  attempt  has  been  made  to  develop  the  necessary  equations  and  programs 
to  portray  a  different  set  of  graph  coordinates.   Thus  only  one  basic 
graphical  presentation  is  used  for  parameter  plane  stability  studies.   The 
example  in  section  2  illustrated  this  method  of  stability  study  by  employ- 
ing the  following  alpha  -  beta  parameter  plane  curves:   constant  zeta  (  £  ) 
curves;  constant  omega  (  ^  n   )  curves;  constant  zeta  -  omega  (  £-6^n) 
curves;  constant  sigma  (  CT'   )  curves;  constant  omega  (  lP      )  curves;  and 
real  root  (  CT  +  J  C  )  curves.   Although  the  constant  zeta  -  omega  ( $-<^r>) 
curves  are  identical  with  the  constant  sigma  (   CT   )  curves,  the  alge 
braic  methods  required  for  each  are  different.   The  names  assigned  to  each 
have  therefore  been  kept  separate  to  avoid  confusion  in  the  respective 
graphical  construction  techniques.   Additionally,  the  real  root  (  (X  +  J  °    ) 
curves  represent  constant  sigma  points  in  the  complex  frequency  plane  which 
are  transformed  into  lines  on  the  parameter  plane.   As  such,  they  are  not 
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to  be  confused  with  the  constant  sigma  (  CT'  )  or  constant  zeta  -  omega 
(   c  -  uJ  n    )  curves.   Although  partial  programs  exist  for  these  curves,  a 
program  entitled  PARAM-D  has  been  written  for  the  primary  purpose  of 
consolidating  these  programs  into  one  package  for  added  flexability  and 
versatility.   This  permits  any  combination  of  the  above  curves  to  be  pro- 
cessed on  one  graph. 

Complex  engineering  problems  often  involve  a  stability  study  of  systems 
with  more  than  two  arbitrary  characteristic  equation  coefficients.   Pro- 
gram PARAM-D  has  been  extended  to  include  a  third  variable  which  is  here- 
in defined  as  gamma  (  Y     ) .   This  scheme  will  allow  stability  analysis 
studies  to  be  visualized  in  three  dimensions,  rather  than  two.  When  con- 
fronted with  this  type  of  problem,  it  then  becomes  essentially  a  four 
variable  problem.   Therefore  it  is  necessary  to  hold  one  of  these  four 
independent  variables  constant,  and  generate  a  series  of  curves  using 
the  remaining  three.   This  is  accomplished  by  assigning  a  constant  value 
to  one  of  the  stability  performance  variables:  f  u)^     ^~6<J  <r  to    or  cr+'jO. 
The  remaining  three  variables  are  those  in  the  characteristic  equation  co- 
efficients, namely  alpha,  beta  and  gamma.   With  these,  the  program  at- 
temps  to  develop  a  three  dimensional  surface  for  the  fixed  value  of  the 
specified  stability  performance  variable.   The  x,  y  and  z  axis  variables 
are  alpha,  beta  and  gamma  respectively.   Graphical  presentations  are  re- 
stricted to  two  dimensions,  however.   The  dilemma  therefore  is  that  this 
three  dimensional  surface  must  be  projected  onto  a  two  dimensional  plane, 
which  is  herein  defined  by  the  x  and  y  axes.   The  only  way  that  this  pro- 
jection can  be  accomplished  is  to  construct  a  large  number  of  third  vari- 
able (gamma)  curves  on  the  graph.   Each  curve  then  represents  a  cross 
section  of  the  three  dimensional  surface.   The  physical  shape  of  the  sur- 
face can  then  be  visualized  providing  a  sufficiently  large  number  of 
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these  curves  are  thus  constructed. 

The  frequency  response  study  has  been  viewed  in  a  manner  which  is 
slightly  different  than  the  above  stability  program  approach.   It  was 
initially  decided  to  make  a  thorough  analysis  of  the  four  variables:  OC   , 
8   »  M  .  and   KJ  .   As  such,  alpha  and  beta  can  be  regarded  as  two 
variables  embedded  in  the  characteristic  equation  coefficients,  while 
magnitude  and  omega  are  the  two  frequency  response  performance  criteria 
of  interest.   Treating  these  as  four  independent  variables  results  in 
seven  basic  graphical  combinations.   These  have  been  titled  PARAM-1 
through  PARAM-7.   Additionally,  phase  is  investigated  in  PARAM-7,  where 
its  correlation  with  magnitude  is  most  easily  made.   All  possible  graphi- 
cal presentations  of  interest  have  been  utilized,  and  a  summary  of  these 
is  presented  in  figure  3-1.   They  have  been  implemented  in  one  consoli- 
dated computer  program  PARAMS  which  is  included  in  appendix  I.   A  brief 
summary  of  these  graph  characteristics  is  presented  below,  while  a  more 
thorough  discussion  is  reserved  for  article  3.2. 

PARAM-1  is  a  graph  of  alpha  versus  beta,  with  magnitude  held  at  a 
fixed  value  while  drawing  constant  omega  curves.   Interchanging  alpha  and 
beta  produces  absolutely  no  alteration  of  curves,  other  than  a  ninety  de- 
gree re-orientation.   PARAM-2  is  a  graph  of  alpha  versus  beta  again,  with 
omega  held  at  a  fixed  value  this  time,  while  constant  magnitude  curves  are 
drawn.   As  in  PARAM-1,  no  difference  other  than  a  ninety  degree  re-orienta- 
tion is  obtained  when  the  axis  variables  are  interchanged.   PARAM-3  places 
alpha  on  the  ordinate,  and  magnitude  on  the  abcissa.   Omega  is  held  at  a 
fixed  value  throughout  the  computational  process,  while  constant  beta 
curves  are  generated.   Now  interchanging  alpha  and  beta  in  PARAM-3  will 
result  in  an  entirely  different  set  of  curves,  since  the  sensitivity  of 
these  two  variables  is  normally  different.   PARAM-4  is  a  graph  of  alpha 
on  the  ordinate  versus  omega  on  the  abcissa.   Beta  is  held  constant  while 
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magnitude  curves  are  generated.   Again,  different  curves  will  be  obtained 
when  alpha  and  beta  are  interchanged,  since  holding  one  characteristic 
equation  parameter  constant  will  produce  an  entirely  different  effect  if 
another  parameter  is  held  at  a  fixed  value  instead.   PARAM-5  generates 
constant  omega  curves  with  beta  held  at  a  fixed  value,  and  graphs  these 
curves  against  magnitude  on  the  abcissa,  and  alpha  on  the  ordinate.   As 
with  PARAM-4,  the  PARAM-5  curves  are  quite  different  when  the  roles  of 
alpha  and  beta  are  interchanged.   PARAM-6  is  a  graph  of  alpha  on  the  ord- 
inate and  omega  on  the  abcissa.   Magnitude  is  held  at  a  fixed  value  through- 
out the  calculations  while  constant  beta  curves  are  generated  on  the  alpha- 
omega  plane.   This  is  the  same  arrangement  as  PARAM-3  with  the  roles  of 
magnitude  and  omega  being  reversed.   Interchanging  the  alpha  and  the  beta 
in  PARAM-6  will  also  cause  a  sharp  alteration  in  the  curve  configurations. 
The  familiar  Bode  plot  is  herein  reproduced  as  PARAM-7.   Needing  no  re- 
view, it  is  simply  a  logarithmic  plot  of  magnitude  on  the  ordinate,  and 
omega  along  the  abcissa.  With  beta  held  fixed  at  some  constant  value, 
constant  alpha  curves  are  generated.   Conversely,  when  alpha  is  held 
fixed  at  some  constant  value,  the  beta  curves  are  generated  and  plotted 
on  the  magnitude  -  omega  plane.   One  problem  exists  with  the  logarithmic 
plot,  however.   In  order  to  draw  a  useful  logarithmic  grid,  it  was  neces- 
sary to  shift  the  origin  to  the  lower  right  hand  corner  of  the  graph  paper. 
In  this  manner,  the  decade  span  could  cover  a  total  of  15  inches  rather 
than  a  maximum  of  nine.   The  net  effect  has  been  to  require  the  abcissa 
to  be  plotted  along  the  right  hand  edge  of  the  graph  while  the  ordinate 
is  plotted  along  the  bottom  where  the  abcissa  otherwise  would  be,  thereby 
coinciding  with  the  X-axis.   Thus  the  ordinate  originates  in  the  bottom 
right  hand  corner  and  increases  toward  the  lower  left  hand  corner,  or 
original  origin.   The  same  process  is  required  for  PARAM-4  and  PARAM  6 
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whenever  specifying  the  logarithmic  plots.   Detailed  coverage  of  the 
available  options  will  be  found  in  article  3.2  which  is  devoted  to 
program  PARAMS. 

The  programs  PARAM-D  and  PARAMS  utilize  only  four  variables,  al- 
though their  adaptation  to  others  has  been  made.   The  method  of  solution, 
however,  has  been  with  two  variables  along  the  axes,  while  the  remaining 
two  have  been  employed  as  constant  curves,  and  in  the  case  of  program 
PARAMS  one  variable  has  been  ehld  fixed  at  a  constant  value.   In  program 
PARAM-D,  the  axis  variables  have  been  selected  as  alpha  and  beta  only; 
no  attempt  has  been  made  to  develop  any  other  projection  whereby  other 
variables  might  replace  alpha  and/or  beta  as  axis  coordinates.   Program 
PARAMS  has  attempted  to  focus  attention  on  several  two  dimensional  pro- 
jections in  the  hopes  that  peculiarities  not  evident  in  one  planar 
presentation  might  become  conspicuous  in  another  mode.   Due  to  the  ex- 
tensive additional  projections  with  the  incorporation  of  one  more  vari- 
able, it  was  decided  to  refrain  from  including  phase  (  0  )  as  a  para- 
meter in  all  but  the  Bode  diagram,  PARAM-7,  where  its  utility  and  value 
have  already  been  demonstrated. 

The  following  article  will  develop  the  basic  capabilities  of  program 
PARAMS  with  regards  to  its  application  in  frequency  response  studies. 

3.2  Digital  Computer  Program  PARAMS. 

The  preceding  article  and  earlier  example  in  section  2  presented  in 
some  detail  the  capabilities  of  program  PARAM-D  for  stability  studies. 
Program  PARAMS  has  been  developed  to  study  some  pertinent  factors  in  fre- 
quency response  which  might  prove  useful  in  systems  and  controls  design. 

For  each  of  the  graphical  representations  entitled  PARAM-1  through 
PARAM-7,  a  separate  computer  program  was  developed  in  FORTRAN  II  language 
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for  use  on  the  Computer  Data  Corporation  (C.D.C.)  Model  1604  digital 
computer  [11].   Each  of  these  programs  was  compiled  in  a  similar  fashion, 
and  employed  the  algebraic  methods  derived  by  Doctors  G.  J.  Thaler  and 
A.  G.  Thompson  [12].   As  previously  discussed  in  this  section,  each  graphi- 
cal configuration  has  been  developed  in  order  that  the  controls  design 
engineer  might  gain  additional  tools  to  facilitate  studies  in  analysis 
and  synthesis.   For  the  study  of  one  particular  problem,  it  is  often 
necessary  to  analyze  the  systems  behavior  from  uiany  different  angles,  and 
then  select  a  smaller  number  of  characteristics  for  further  intensive 
study.   Unfortunately,  every  real  life  problem  involves  a  different  set 
of  constraints,  and  for  that  reason  there  can  be  no  set  pattern  for  suc- 
cessful design  application  in  all  cases.   Many  situations  reduce  to  trial 
and  error  techniques  on  preliminary  studies  in  order  to  obtain  the  best 
compromise  when  a  number  of  possible  solutions  exist.   The  ability  to  ac- 
complish this  objective  using  many  separate  computer  programs  is  indeed 
possible  but  cumbersome.   Instead  of  this  method,  one  consolidated  com- 
puter program  titled  PARAMS  has  been  written  in  FORTRAN  IV  language  for 
use  on  the  International  Business  Machine  (IBM)  model  360  computer.   The 
entire  computer  program,  together  with  explanatory  flowgraphs  and  illustra- 
tive printouts,  has  been  included  as  Appendix  I.   Details  of  its  use,  to- 
gether with  the  advantages  and  disadvantages,  will  now  be  presented. 

Versatility  is  indeed  a  virtue,  provided  it  does  not  gain  domination 
over  simplicity.   While  considerable  apprehension  might  be  the  first  re- 
action upon  reviewing  Appendix  I,  a  systemized  pattern  will  become  evident 
upon  further  inspection.   With  this  program,  complexity  does  not  enter  into 
the  computational  procedures,  but  rather  a  considerable  amount  of  detail 
has  been  included  in  the  printout  format  for  the  sake  of  clarity  when  re- 
viewing the  material  contained  within  the  graphs.   At  most,  the  computational 
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process  is  intricate  in  the  initializing  stage  wherein  the  various  axis 
maximum  and  minimum  limits,  together  with  the  curve  and  fixed  constant 
values  are  redesignated  for  the  ensuing  common  computational  and  graphi- 
cal construction  sections.   This  can  be  verified  by  noting  that  the  numera- 
tor and  denominator  coefficient  calculation  sections  are  common  to  PARAM-1 
through  PARAM-6.   It  will  further  be  seen  that  the  gain  computation  for 
PARAM-7  employs  different  formulae  but  precisely  the  same  logic,  and  for 
that  reason  is  incorporated  in  the  same  numerator  and  denominator  coef- 
ficient calculation  sections  as  used  for  PARAM-1  through  PARAM-6.   Only  the 
phase  calculations  for  PARAM-7  require  a  different  process  for  generation 
of  the  coefficients.   For  this  reason,  this  computational  section  is  by- 
passed for  all  other  processes. 

The  logical  commencement  point  for  design  studies  is  to  define  a 
range  of  parameter  values  which  guarantees  all  roots  to  lie  in  the  left 
half  plane,  for  clearly  the  question  of  stability  must  remain  foremost  in 
the  designer's  mind.   Program  PARAM-D  has  been  devised  to  meet  this  demand, 
and  its  capabilities  to  provide  a  stability  perspective  have  already  been 
delineated.   In  determining  the  various  possible  frequency  response 
characteristics  for  a  given  system,  several  graphical  representations 
will  be  invariably  studied.   Most  likely,  however,  program  PARAM-7  will 
be  employed  in  the  final  analysis  as  an  overall  check  on  system  specifica- 
tion criteria.   The  remaining  programs,  PARAM-1  through  PARAM-6  are  most 
informative  in  the  initial  and  intermediate  stages  of  design  study.   It 
would  therefore  seem  that  the  inclusion  of  all  frequency  response  para- 
meter plane  programs  into  one  package  would  be  most  useful.   Different 
design  studies  will  normally  emphasize  certain  performance  criteria  more 
than  others.   Thus  seldom  would  the  same  sequence  of  graphical  presenta- 
tions be  called  for  in  two  consecutive  studies.   In  fact,  PARAM-7  might  be 
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omitted  for  the  initial  coarse  adjustments,  but  yet  included  in  the  final 
stages  of  refined  improvements.   In  this  manner  each  system  study  can  be 
tailored  to  meet  certain  specifications  and  suit  specific  needs.   Also, 
certain  systems  will  be  more  complex  than  others.   Therefore,  more  than 
one  mathematical  model  might  be  considered  for  simultaneous  inspection. 
The  program  is  constructed  to  handle  any  number  of  such  studies  at  one 
time.   Each  of  these  normally  requires  a  complete  set  of  data  cards,  and 
a  brief  explanation  of  their  usage  will  follow.   Additionally,  an  abbrevi- 
ated description  section  appears  in  the  program  preface,  found  in  Appendix 
I.   The  total  number  of  these  data  card  sets  must  be  indicated  on  the 
first  card  following  the  main  program.   This  desired  number  of  runs  is 
designated  by  the  variable  NRUN.   The  program  simply  processes  the  re- 
maining data  cards  by  commencing  a  "DO  LOOP"  with  the  variable  IRUN  going 
from  1  to  NRUN.   The  first  card  of  each  data  set  indicates  which  graphi- 
cal presentations  (PARAM-1  through  PARAM-7)  are  desired.   Again,  there  is 
no  limit  to  the  number  of  these.   The  program  simply  commences  another 
"DO  LOOP"  at  statement  number  26,  with  IPARAM  going  from  1  to  KPARAM. 
This  cycles  each  of  the  graphical  constructions  in  the  order  called  for. 
In  this  manner,  a  great  multiplicity  of  variations  exist  for  the  simultan- 
eous study  of  different  preoblems.   At  this  point  it  is  pertinent  to  in- 
dicate a  technique  by  which  additional  finesse  may  be  utilized.   Although 
most  studies  involve  extensive  rearrangement  of  a  data  set  from  one  system 
to  another,  some  require  only  a  slight  modification.   In  these  latter  cases, 
it  is  not  necessary  to  submit  an  entire  additional  set  of  cards  for  each 
study.   By  first  specifying  the  desired  graphical  plots  of  one  system,  it 
is  possible  to  consecutively  list  the  desired  sequence  for  another  system. 
Then  a  simple  insertion  of  the  desired  data  modification  immediately  follow- 
ing statement  number  26  will  alter  the  original  data  in  the  specified  manner 
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As  an  example,  suppose  one  desired  to  first  construct  PARAM-5  with  alpha 
on  the  ordinate,  and  then  with  beta  on  the  ordinate.   Simply  list  KPARAM 
as  2,  and  then  indicate  5  twice  for  two  consecutive  plots  of  PARAM-5. 
Insert  the  following  statement  on  an  additional  card  behind  statement  26: 
'IF  (IPARAM.  GE.  2)  MODORD  =  2".   Thus  the  required  change  in  the  input  data 
is  effected  at  the  desired  time  by  one  card  rather  than  a  complete  set  of 
data  cards. 

Although  no  standard  names  have  been  assigned  to  these  different  pro- 
grams the  title  of  each  graphical  representation  (PARAM-1  through  PARAM-7) 
has  been  retained  to  eliminate  confusion  in  terminology  for  those  persons 
previously  exposed  to  the  partial  programs  which  are  already  in  existence. 
The  result  has  been  a  rather  complex  programming  sequence  for  printout  of 
graphical  information  and  preliminary  initializing  computations,  but  as 
has  already  been  stated,  the  basic  manipulations  are  common  to  all  pro- 
grams, and  are  centrally  located  in  the  program  from  statement  number  350 
through  statement  426.   A  cursory  inspection  of  these  76  statement  numbers 
as  well  as  the  remainder  of  the  program  will  reveal  that  its  construction 
is  modular  in  form,  and  each  section  is  a  separate  entity  unto  itself  with 
adequate  illustrative  sub-titles.   The  overall  program  has  been  so  ar- 
ranged that  many  sections  can  be  removed  if  it  is  desired  to  employ  their 
use  elsewhere.   Conversely,  additional  extensions  are  easily  facilitated 
by  this  manner  of  program  fabrication. 

Careful  consideration  has  been  used  in  this  program  design  to  retain 
flexibility  for  use  on  multiple  problem  types.   Although  the  program  is 
slightly  larger  in  size  than  any  of  its  component  programs  which  were 
previously  mentioned,  it  is  smaller  than  the  combination  of  any  two  of 
them.   While  intensive  effort  to  streamline  programs  is  not  of  particular 
significance,  the  repetition  of  computer  processing  is  needless  and 
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inefficient.   Also,  the  preparation  of  data  cards  for  a  large  number  of 
programs  serves  only  to  increase  the  probability  of  error  with  consequent- 
ial waste  in  time,  expense  and  labor.   More  important  is  the  fact  that 
computer  programming  is  still  a  dynamic  field  which  often  necessitates 
program  modifications  to  meet  the  requirements  of  more  sophisticated  equip- 
ment.  Thus  one  might  be  motivated  to  use  a  single  computer  program  in 
lieu  of  eight  separate  ones,  and  to  this  extent  consolidation  is  warranted. 

To  ensure  adequate  clarity  of  results,  an  elaborate  printout  of  in- 
formation has  been  established.   This  information  can  be  basically  cata- 
gorized  into  three  groups,  each  of  which  is  controlled  by  a  single  input 
data  variable  IWRITE.   The  first  reproduces  all  input  data  for  verifica- 
tion, and  is  effected  whenever  the  value  of  IWRITE  is  any  non-negative 
number.   The  second  category  gives  detailed  information  about  how  each 
graph  will  be  constructed.   Included  in  this  section  are  the  graph  titles, 
the  range  of  values  for  the  ordinate  and  abcissa  variables,  the  curves  to 
be  generated,  and  the  parameter  that  is  to  be  held  constant  throughout  the 
computational  process  of  that  particular  program.   This  explanatory  sec- 
tion will  be  printed  for  all  positive  values  of  IWRITE.   Finally,  the  com- 
puted results  for  each  curve  are  printed  out  at  pre-selected  intervals. 
The  value  assigned  to  IWRITE  determines  at  what  intervals  the  computation- 
al results  will  be  printed,  and  the  range  of  permissible  values  is  from 
one  to  900.   Thus  an  assigned  value  of  one  would  result  in  the  reproduction 
of  every  calculated  curve  value,  whereas  a  value  of  ten  would  give  a  print- 
ed value  for  every  tenth  calculated  point.   Approximately  one  third  of  the 
program  is  devoted  to  printing  out  the  above  information  for  feedback  pur- 
poses.  If  desired,  these  sections  can  be  partially  or  entirely  bypassed 
by  the  appropriate  selection  of  the  input  data  variable  IWRITE.   Any  nega- 
tive value  assigned  to  this  variable  will  provide  complete  calculation, 
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but  without  any  printout  whatever. 

It  has  already  been  pointed  out  that  if  the  roles  of  alpha  and  beta 
are  allowed  to  be  interchanged,  a  much  smaller  number  of  graphs  are  re- 
quired.  In  order  to  accomplish  this  objective,  some  criterion  must  exist 
for  the  choice  of  which  variable  (alpha  or  beta)  will  be  plotted  on  the 
ordinate  in  the  case  of  PARAM-1  through  PARAM-6,  and  which  variable  will 
be  held  at  a  fixed  constant  value  while  the  other  generates  the  desired 
curves  in  the  case  of  PARAM-7.   The  input  data  variable  titled  MODORD  has 
been  designated  for  this  purpose,  and  so  by  this  selection  the  type  of 
curves  is  controlled.   With  a  value  of  one  for  MODORD,  alpha  curves  will 
be  generated  in  the  case  of  PARAM-7,  and  alpha  will  be  plotted  on  the 
ordinate  for  the  rest  of  che  programs,   Conversely  beta  will  be  used  in 
lieu  of  alpha  for  those  functions  when  the  value  of  MODORD  is  set  equal 
to  two.   This  flexibility  for  PARAM-1  and  PARAM-2  might  be  unnecessary, 
because  as  indicated  earlier,  a  ninety  degree  reorientation  of  the  graph 
is  the  only  effect  of  interchanging  alpha  and  beta  when  these  variables 
are  on  the  axes.   However,  the  present  day  graphical  plotting  facilities 
are  generally  constructed  in  such  a  manner  that  the  abcissa,  or  X-axis, 
is  more  severely  restricted  in  length  than  the  ordinate,  or  Y-axis.   The 
facilities  available  at  the  inception  of  program  PARAMS  were  in  fact  such 
that  the  ordinate  length  could  be  approximately  twice  as  great  as  the 
length  of  the  abcissa.   It  was  envisioned  that  in  some  cases  an  expansion 
of  one  variable  might  be  more  desirable  than  another.   Extending  the  ap- 
plicability of  MODORD  to  PARAM-1  and  PARAM-2  makes  this  provision  possible, 
which  otherwise  would  require  modification  of  the  main  program  to  achieve 
such  results.   Furthermore,  it  is  advantageous  to  retain  a  common  orienta- 
tion of  graphical  construction  when  viewing  several  projections  simultan- 
eously. 
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The  value  of  magnitude  often  occurs  in  several  different  forms,  and 
thus  it  provides  a  source  of  confusion  if  the  terminology  is  not  clearly 
specified.   The  input  data  variable  IMAG  serves  to  eliminate  this  poten- 
tial error  by  distinctly  categorizing  magnitude  into  one  of  three  forms. 
The  program  is  constructed  in  such  a  manner  that  the  form  of  magnitude  as 
specified  by  IMAG  is  maintained  consistently  for  all  graph  plots.   Al- 
though this  is  strictly  a  matter  of  convenience,  it  was  found  that  various 
approaches  to  design  require  magnitude  to  be  handled  in  different  forms; 
hence  this  inclusion.   Thus  for  magnitude  to  be  interpreted  as  decibels, 
a  value  of  zero  would  be  assigned  to  IMAG,  and  the  curve  values  of  PARAM-2 
and  PARAM-4  would  be  listed  in  decibels,  while  the  abcissa  values  of  PARAM- 
3  and  PARAM-5  would  likewise  be  interpreted  accordingly.   If  actual  magni- 
tude values  were  to  be  used  instead,  IMAG  would  necessarily  have  to  be 
changed  to  one.   Assigning  a  value  of  two  to  IMAG  would  dictate  all  values 
to  be  handled  as  magnitude  squared.   In  all  facets  of  life  it  seems  that 
most  conventions  are  fraught  with  exceptions,  and  this  program  has  not 
escaped  that  dilemma.   Since  all  Bode  plots  associate  magnitude  with  deci- 
bels, it  was  decided  to  maintain  this  practice;  consequently  the  selection 
of  IMAG  has  no  effect  on  the  form  of  PARAM-7  graphs.   They  automatically 
plot  magnitude  along  the  ordinate  in  the  form  of  decibels.   Conventional 
practice  has  also  dictated  that  omega  be  logarithmically  plotted  along  the 
abcissa  for  all  Bode  diagrams.   Other  frequency  response  graphs  which  also 
plot  omega  on  one  axis  have  not  acquired  any  standardized  form,  and  for  that 
reason  provisions  have  been  made  to  select  a  square  linear  grid  or  a  logari- 
thmic grid  as  deemed  appropriate.   The  input  data  variable  IGRAPH  provides 
this  option  for  the  various  possible  combinations  of  PARAM-4,  PARAM-6 
and  PARAM-7.   Additionally,  this  variable  is  used  to  prevent  any  graphs  from 
being  generated  by  assigning  to  it  a  value  of  minus  two.   This  provision  is 
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useful  only  for  efficient  use  of  available  hardware  when  desired  graph 
scales  are  initially  unknown. 

In  order  that  maximum  versatility  might  be  gained,  the  memory  stor- 
age requirements  of  this  program  have  been  expanded  to  present  day  capa- 
bilities.  A  review  of  the  program  DIMENSION  statement  just  preceding  the 
input  data  read  in  section  will  show  that  the  bulk  of  the  storage  capa- 
city is  consumed  by  the  two  dimensional  variables  S   through   W. 
These  are  employed  for  the  various  optional  arrangements  of  PARAM-7  and 
may  be  combined  into  one  common  storage  array  for  all  other  calculations 
when  employing  this  program  in  smaller  computer  facilities.   In  this 
manner,  the  memory  core  may  be  reduced  from  500K  bytes  to  less  than  200K 
bytes.   This  is  accomplished  by  inserting  the  following  statement  after 
the  dimension  statement:   EQUIVAIANCE  (S,T,U,V,W). 

The  additional  variable,  phase,  has  been  incorporated  in  PARAM-7 
with  several  variations  possible.   The  input  data  variable  IPHASE  is 
used  to  govern  the  relation  between  the  magnitude  and  phase  plots. 
No  phase  curves  will  be  drawn  when  IPHASE  is  set  equal  to  a  negative  num- 
ber.  Setting  IPHASE  equal  to  zero  will  enable  phase  curves  to  be  super- 
imposed upon  the  magnitude  curves  for  ease  in  investigating  stability 
by  means  of  phase  and  gain  margins.   Assigning  any  positive  number  to 
IPHASE  will  provide  separate  phase  curves.   To  allow  for  various  super- 
positions, and  various  separate  phase  curve  configurations  as  well,  it 
was  decided  to  employ  an  additional  variable  IQUAD  which  allowed  the  phase 
graphs  to  be  plotted  in  the  following  two  manners:   When  IQUAD  has  an  as- 
signed value  of  zero,  phase  curves  will  be  automatically  scaled  at  45 
degrees  per  inch,  with  180  degrees  on  the  magnitude  axis.   When  the  vari- 
able IQUAD  is  positive,  however,  the  phase  curves  will  be  plotted  with  a 
scale  equal  to  the  value  of  IQUAD,  starting  with  0  degrees  at  the  top  of 
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the  graph.   In  this  manner,  it  is  possible  to  cover  any  desired  range  of 
phase  angles. 

Basic  conformal  mapping  theory  indicates  that  a  polar  plot  can  be 
generated  from  the  information  contained  in  the  magnitude  and  phase  plots. 
Particular  importance  might  be  attached  to  the  polar  plot  representation 
of  a  system  when  the  Nyquist  criterion  is  preferred  over  the  Bode  stabil- 
ity analysis  method.   Although  it  is  a  duplicative  process  and  therefore 
not  strictly  required  for  information  content,  its  inclusion  has  been 
assumed  appropriate  for  the  sake  of  completeness.   The  basic  construc- 
tion of  the  polar  plane  and  polar  curves  are  contained  within  the  program 
statement  numbers  537  and  454  under  the  heading  "PARAM-7  Polar  Plot  Curves". 
Its  usage  is  governed  by  the  rather  appropriate  input  data  variable  IPOLAR, 
and  a  polar  plot  will  be  drawn  when  its  value  is  one.   Conversely,  no 
polar  plot  will  be  drawn  for  any  other  value  of  IPOLAR. 

An  additional  provision  has  been  incorporated  in  program  PARAMS 
for  extended  usage  of  the  Bode  plot.   Certain  advanced  feedback  design 
techniques  call  for  complex  systems  to  be  studied  from  a  cumulative  stand- 
point.  Such  design  measures  require  vector  addition  of  component  Bode 
diagrams  and  as  such,  the  linear  addition  of  each  magnitude  plot  and  each 
phase  plot  will  yield  erroneous  results.   Instead,  the  vector  quantities 
of  each  design  segment  are  summed  in  complex  notation,  and  the  final  re- 
sult is  then  plotted  as  a  single  magnitude  curve,  and  also  a  single  phase 
curve  when  the  latter  is  requested.   The  program  subsection  designated  to 
accomplish  this  process  is  titled  "PARAM  -  7  Superposition  Section",  and 
is  contained  by  the  program  statement  numbers  520  and  529.   It  is  con- 
trolled by  the  input  data  variable  NPAR7 .  where  the  magnitude  of  this 
variable  dictates  how  many  component  terms  are  desired  to  be  added  toget- 
her.  It  is  necessary  to  emphasize  at  this  point  that  this  section  performs 
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a  vector  superposition  of  feedback  terms,  and  not  merely  a  scalar  super- 
position.  Thus  whenever  NPAR7  is  greater  than  one,  this  section  will 
add  NPAR7  minus  one  terms  vectorially  to  the  first  term  before  graphing 
the  results.   In  order  that  these  components  may  be  properly  calculated, 
provisions  have  been  made  to  insert  the  component  equation  terms  in  the 
same  coefficient  format  as  normally  submitted  for  the  conventional  cal- 
culations.  These  additional  terms,  however,  are  appended  to  the  initial 
group  of  input  data  cards.   Therefore  NPAR7  minus  one  groups  of  coeffici- 
ents must  be  submitted  directly  behind  the  title  cards  of  the  rest  of  the 
program. 

A  brief  summary  of  the  basic  operational  procedures  which  are  common 
to  each  of  the  programs  will  now  be  undertaken.   The  schematic  corollary 
in  appendix  I  will  parallel  the  following  outline  and  should  be  consulted 
together  with  the  main  program  for  more  detail. 

The  first  six  statements  in  the  "Input  Data  Read  in  Section"  pro- 
vide in  proper  FORTRAN  language,  the  necessary  format  for  most  of  the 
read  and  write  statements  of  the  program.   Their  use  is  specified  in  the 
"Data  Card  Information  Section"  which  appears  at  the  beginning  of  the 
program.   The  remainder  of  the  input  data  read  in  section  is  devoted  to 
reading  all  the  19,  21,  or  23  data  cards  as  outlined  in  the  preceding  in- 
formation section.   The  next  step  is  to  print  out  the  input  data  when  so 
requested  for  verification.   This  procedure  is  accomplished  by  the  section 
bounded  by  the  statement  numbers  seven  through  23,  inclusive.   Following 
the  input  data  verification  section,  the  program  initializes  the  various 
parameters  and  computes  the  necessary  controlling  variables  from  such  in- 
put data  as  axes  maximum  and  minimum  values.   Within  this  section,  con- 
tained by  statements  24  to  54,  the  desired  value  of  NPARAM  is  selected  and 
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this  value  guides  the  rest  of  the  program  sequence  by  assigning  the  proper 
axis  variables  and  their  limits,  the  correct  number  and  type  of  curves  to 
be  generated  etc.   When  the  logarithmic  plot  is  requested  for  PARAM  -  4 
or  PARAM  -  6,  and  whenever  the  Bode  plot  is  to  be  generated,  the  section 
titled  "Log  Grid  Graphical  Construction"  is  executed.   Extending  from 
statement  number  55  to  58,  this  sub-section  utilizes  the  information 
generated  in  the  initializing  section  to  determine  the  ordinate  scale 
and  origin  shift  if  applicable.   Then  it  calculates  the  required  number 
of  decades  by  rounding  up  the  maximum  limit  of  omega  to  the  nearest  power 
of  ten,  and  rounding  down  the  minimum  omega  value  to  its  nearest  power  of 
ten.   Utilizing  these  factors  and  the  axes  lengths,  a  logarithmic  grid  is 
constructed  in  such  a  manner  that  the  omega  values  are  properly  labelled 
along  the  right  hand  edge  of  the  graph  paper,  while  the  ordinate  values 
are  indicated  across  the  bottom  in  such  a  manner  that  their  magnitudes  are 
correct  but  the  signs  reversed.   Thus  the  ordinate  scale  must  be  correct- 
ed to  read  increasing  from  right  to  left.   Although  this  arrangement  is 
less  than  optimum,  it  has  proved  to  be  much  better  than  some  alternatives 
which  were  previously  attempted. 

The  most  prolific  section  of  the  program  ranges  from  statement  number 
59  to  200  under  the  sub  heading  "Graphical  Information  Printout  Section". 
Its  utility  lies  in  the  desire  to  establish  positive  correlation  between 
printed  data  and  graphical  presentations  when  several  similar  graphical 
analyses  are  computed  simultaneously.   Following  the  extensive  information 
printout,  the  main  computational  process  commences  by  initiating  calcula- 
tions for  the  first  curve  with  statement  number  202.   Next,  the  coefficient 
preceding  the  quadratic  formula  square  root  is  assigned  by  the  variable 
ZSIGN.   Following  this  designation,  the  abcissa  is  incremented  900  times 
from  the  minimum  specified  value  to  the  maximum  value.   Additionally  the 
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curve  value  is  selected,  based  upon  the  value  of  the  variable  in  statement 
number  202.   Then  the  appropriate  parameter  which  is  to  be  held  at  a 
fixed  constant  value  is  assigned  its  pre-selected  value.   The  ensuing 
calculations  commencing  with  statement  number  350  are  in  conformance  with 
the  theory  presented  in  section  2.   The  magnitude  protion  of  the  Bode 
program  requires  the  coefficients  designated  SUMN  for  the  numerator  and 
SUMD  for  the  denominator,  while  the  remaining  programs  utilize  the 
numerator  variables  SNB2,  SNB1,  and  SNBO  in  addition  to  the  denominator 
variables  SDB2,  SDB1,  and  SDBO.   Such  requirements  are  generated  under 
the  sub  headings  "Numerator  Coefficients  Calculations"  and  "Denominator 
Coefficient  Calculations",   Following  these  sections,  the  necessary  phase 
coefficients  are  computed  from  statement  number  386  through  400,  and  the 
additional  Bode  curve  ordinate  value  calculations  are  continued  through 
statement  number  410.   The  section  titled  "Quadratic  Formula  Solution  of 
Dependent  Variable,  PARAM-1  through  PARAM-6"  is  self  explanatory  and 
simply  computes  the  ordinate  value  commensurate  with  the  assigned  fixed 
value  of  the  parameter  which  is  kept  constant,  together  with  the  specified 
value  of  the  curve  which  is  being  generated,  and  the  corresponding  point 
on  the  abcissa.   The  computed  results  are  then  printed  in  accordance  with 
the  provisions  contained  between  statement  number  426  and  448.   Having 
computed  this  ordinate  value,  the  abcissa  which  is  the  independent  vari- 
able is  again  incremented  and  a  corresponding  ordinate,  or  dependent  vari- 
able point  is  calculated.   After  reaching  the  maximum  permissible  abcissa 
value,  the  process  is  again  repeated  for  the  alternate  solution  of  the 
quadratic  formula.   When  both  possible  ranges  of  values  have  been  scanned, 
the  next  curve  is  specified,  and  the  entire  procedure  repeated.   Upon  ter- 
minating the  curve  requirements,  the  program  then  reverts  to  statement 
number  26  for  the  next  cycle  in  the  run.  When  all  requested  graphical 
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construction  cycles  have  been  accomplished,  the  next  run  is  processed 
in  like  manner.   The  overall  program  is  terminated  when  all  desired  runs 
have  been  computed  and  graphed. 

The  graphical  plotting  procedures  will  be  discussed  separately  be- 
cause their  intricate  operations  tend  to  destroy  the  continuity  of  the 
preceding  discussion  on  the  overall  program  operation.   Basically,  the 
graphical  plotting  subroutine  is  executed  whenever  the  following  condi- 
tions are  engaged:   1)  upon  initially  setting  up  the  logarithmic  plot 
when  specified  in  accordance  with  the  above  dictates;  2)  when  two  or  more 
computed  points  fall  within  the  prescribed  graph  range  and  the  next  point 
lies  outside  the  ordinate  graph  limits;  3)  when  the  end  of  each  scan  of 
the  abcissa  has  been  reached  -  once  for  each  value  of  the  quadratic  equa- 
tion square  root  coefficient;  4)  at  the  end  of  each  set  of  curves  for  one 
graph,  as  accomplished  by  the  graphical  termination  section  found  at  the 
end  of  the  program  beginning  with  statement  number  999. 

The  versatile  arrangement  necessitates  several  provisions  for  the 
many  possible  combinations  of  the  Bode  plot.   Since  the  quadratic  formula 
is  not  employed  in  these  calculations,  only  one  incrementing  scan  of  the 
abcissa  is  performed.   Additionally,  the  computed  magnitude  values  have 
been  programmed  in  such  a  manner  that  when  the  ordinate  graph  limits  are 
exceeded,  the  external  points  are  fictitiously  redefined  to  lie  on  the 
appropriate  ordinate  graph  limits.   This  procedure  was  necessitated  by  the 
additional  phase  calculations  which  remained  within  the  graph  limits 
throughout  the  entire  range  of  the  prescribed  omega  values.   Finally,  a 
separate  provision  had  to  be  made  for  the  construction  of  the  phase  curves, 
since  an  additional  logarithmic  graph  was  required  when  phase  curves  were 
plotted  separately.   This  procedure  has  been  accordingly  programmed  in  the 
section  titled  "PARAM  -  7  Phase  Curves"  found  between  statement  numbers  529 
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and  537.   The  section  titled  "PARAM  -  7  Polar  Plot  Curves"  has  already 
been  discussed  and  is  a  separate  entity  unto  itself,  utilizing  the  pre- 
computed  values  of  magnitude  and  phase,  and  plotting  same  in  polar  coordi- 
nates.  This  process  is  accomplished  between  statement  numbers  537  and  999 
in  the  program. 

The  preceding  development  has  summarily  shown  that  this  program  con- 
tains broad  applications,  and  is  well  suited  for  additional  modifica- 
tions to  meet  particular  needs.   Blanket  provisions  for  another  graphical 
configuration  has  been  purposely  included  in  the  program  sequence  labelled 
NPARAM-  8,  and  should  future  expansion  be  deemed  appropriate,  the  required 
modification  will  merely  entail  filling  the  appropriate  spaces  with  the 
desired  initializing  variables  and  applicable  computational  equations. 
To  facilitate  such  extensions  and  also  to  enable  maximum  comprehension 
of  the  program  operations,  two  logic  flow  diagrams  have  been  included  in 
Appendix  I.    The  first  is  a  generalized  sequentail  listing  of  logical 
operations  which  are  applicable  tc  all  graphical  configurations.   The 
second  is  a  detailed  flow  chart  for  PARAM  -  6,  so  chosen  because  it  in- 
cludes most  of  the  alternatives  presented  in  this  discussion. 

A  final  comment  is  appropriate  concerning  the  application  of  computer 
program  PARAMS.   As  indicated  earlier  in  this  article,  considerable  atten- 
tion has  been  given  to  versatility.   The  associated  price  of  computational 
time  and  memory  storage  has  consequently  been  paid.   The  procedures  to  be 
used  in  reducing  memory  storage  for  PARAM-1  through  PARAM-6  has  already 
been  indicated.   Compilation  time  is  roughly  1.5  minutes.   An  additional 
30  seconds  is  consumed  for  each  printout  of  preliminary  graphical  informa- 
tion.  The  actual  calculation  of  16  curves  for  each  graph  requires  approxi- 
mately 3  minutes.   Thus  a  considerable  amount  of  computer  time  is  consumed 
when  a  complete  set  of  curves  is  requested  for  all  possible  graphical 
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representations.   Various  methods  exist  to  reduce  this  time  requirement. 
By  submitted  an  object  deck,  compilation  time  is  deleted.   This  imposes 
the  restriction  of  submitting  a  complete  data  card  set  for  each  run, 
however.   Eliminating  the  input  data  printout,  and  the  preliminary  graphi- 
cal information  printouts  will  also  reduce  the  time  requirements  some- 
what.  This,  however,  is  not  recommended  because  the  elaborate  informa- 
tion presentation  provides  an  exceedingly  valuable  check  on  the  validity 
of  the  graphical  results.   Finally,  the  most  time  is  consumed  by  a  highly 
sensitive  incrementation  of  the  abcissa  variable.   It  will  be  recalled 
that  the  range  of  the  abcissa  variable  is  divided  into  900  evenly  spaced 
steps.   For  each  of  these  abcissa  values,  a  corresponding  ordinate  value 
is  computed.  Although  many  situations  lend  themselves  to  a  larger  in- 
crementation, it  was  decided  to  maintain  the  fine  step  sizes  for  highly 
sensitive  system  curves.   The  maximum  sensitivity  of  900  points  per  graph 
line  is  the  restriction  placed  upon  the  program  by  the  CALCOMP  plotter  in 
use  at  the  Naval  Postgraduate  School  computer  facility.   Reduction  of  this 
incrementation  sensitivity  is  therefore  another  means  by  which  computation- 
al time  can  be  further  reduced.   The  associated  loss  of  information  con- 
tent must  be  considered  in  such  cases.   As  will  be  pointed  out  in  the 
next  section,  not  all  PARAMS  frequency  response  curves  would  be  required 
or  desired  in  every  design  study.   The  judicious  choice  of  the  most  valu- 
able plots  for  a  given  set  of  constraints  will  normally  eliminate  approxi- 
mately half  of  the  program  possibilities.   Only  by  experience,  insight  and 
intelligent  design  procedures  can  a  good  trade-off  be  achieved  between  the 
tools  available  and  the  data  required. 
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4.    Tachometer  and  Acceleration  Feedback  Compensation  Studies. 

The  mathematical  developments  of  section  2  have  provided  the  basis 
for  constructing  the  computer  programs  which  were  discussed  in  Section  3. 
This  section  will  apply  these  methods  to  a  system  study  employing  tacho- 
meter and  acceleration  feedback  as  compensation  techniques.   The  frequency 
response  parameter  curves  will  be  analyzed  according  to  the  respective 
utility  of  each  as  a  design  tool.   Meticulous  detail  is  employed  in  the 
initial  discussion  of  these  curves.   This  is  to  provide  a  clear  under- 
standing of  the  limitations  and  advantages  of  each  respective  graphical 
representation. 

4. 1  Third  order 

Although  an  entire  library  of  control  engineering  studies  has  been 
created,  there  has  been  no  attempt  to  distinguish  a  set  of  example  prob- 
lems as  superlative  for  generalized  case  studies.   Therefore  a  simple 
third  order  system  has  been  considered  for  an  initial  application  of  the 
parameter  plane  frequency  response  curves.   The  illustrative  example  which 
was  introduced  in  section  2  will  now  be  explored  in  depth.   It  will  be  re- 
called that  Figure  2-3  presented  the  schematic  diagram  of  the  third  order 
system  with  compensation  included.   The  stability  curves  for  the  system 
where  mapped  into  the  parameter  plane  shown  in  Figure  2-5.   In  that  pre- 
sentation alpha,  the  abcissa  variable,  represented  acceleration  feedback 
(  Y\  ct    ).   Similarly  beta,  the  ordinate  variable,  indicated  tachometer 
feedback  (  Kt  ).   Reference  to  the  characteristic  equation  of  the  sys- 
tem shown  in  Figure  2-3  shows  that  the  steady  state  velocity  error  is 
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Inspection  of   Kv   shows  that  the  smallest  allowable  value  of  beta 
will  minimize  the  steady  state  velocity  error.   It  is  presumed  that  this 
is  one  of  the  predominant  performance  criteria.   Furthermore,  one  can 
arbitrarily  stipulate  the  dominant  complex  roots  to  have  a  damping  ratio 
(   ?   )  of  0.5.  With  these  system  stability  requirements,  one  can  see 
that  optimum  performance  will  occur  at  point  A  in  Figure  2-5.   Thus  the 
optimal  value  of  beta  is  0.036,  while  the  optimal  value  of  alpha  is 
0.00047. 

A  tacit  assumption  up  to  this  point  has  been  that  frequency  perform- 
ance constraints  are  of  minor  importance.   Most  design  studies  do  not 
conform  to  this  type  of  ommission;  however  stability  is  of  interest  in 
all  systems.   It  has  therefore  been  assumed  that  the  two  stability  per- 
formance criteria  above  are  desirable  but  not  mandatory.   This  permits 
frequency  response  specifications  to  be  considered  with  some  initial 
reference  point.   Figure  4-1  shows  an  approximate  correlation  between 
frequency  response,  transient  performance,  and  root  locations  in  the  com- 
plex frequency  plane.   Additionally  the  separate  root  loci  with  alpha  and 
beta  varied  in  Figures  4-2  and  4-3  respectively  show  the  approximate  vari- 
ation of  root  locations  as  a  function  of  one  parameter  only.   The  corres- 
ponding effect  in  the  frequency  domain  is  represented  by  the  Bode  plots 
in  Figures  4-4  through  4-7.   The  graphical  relationships  presented  in 
Figures  4-1  through  4-7  are  certainly  not  new  to  the  design  engineer. 
Rather,  they  are  the  fundamental  tools  used  in  most  system  studies  today. 
As  such,  they  are  included  primarily  to  establish  a  firm  foundation  upon 
which  parameter  plane  frequency  response  curves  can  be  built.   In  this 
manner,  a  better  understanding  can  be  gained  concerning  the  utility  of 
the  latter  curves  which  are  unfamiliar  to  most  engineers. 
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PARAM-1 

Complete  analysis  of  the  information  contained  in  PARAM-1  plots 
will  now  be  endeavored.   Figure  4-8  shows  the  relationships  between  alpha, 
beta  and  omega  which  exist  for  a  constant  magnitude  of  12  decibels.   It 
is  clear  from  the  legend  at  the  bottom  of  the  graph  that  alpha  is  the 
abcissa  variable,  while  beta  values  are  indicated  on  the  ordinate.   Each 
curve  in  the  alpha  -  beta  plane  represents  a  fixed  magnitude  of  12  deci- 
bels, but  a  different  value  of  omega.   Thus  the  curve  of  UJ    =10  rad/sec 
barely  appears  in  the  lower  right  hand  corner.   As  omega  increases,  the 
curves  can  be  seen  to  shift  toward  higher  values  of  beta  and  lower 
values  of  alpha.   It  can  also  be  seen  that  the  area  enclosed  by  each 
curve  decreases  as  omega  increases  to  U)   =  70  rad/sec.   The  meaning  of 
each  of  these  14  curves  can  be  explained  as  follows:  For  the  third  order 
system  in  question,  any  point  on  a  specified  omega  curve  depicts  the 
locus  of  alpha  and  beta  values  which  will  guarantee  a  magnitude  of  12 
decibels  at  that  specific  frequency.   This  by  itself  does  not  portray 
much  information,  because  it  says  nothing  about  the  area  enclosed  by  the 
curves,  nor  the  area  outside  the  locus  of  these  curves.   In  most  situations, 
it  is  possible  to  glean  considerably  more  information  out  of  these  curves. 
Consider  for  a  moment  that  value  of  alpha  and  beta  at  which  two  omega 
curves  intersect.   Point  B  of  Figure  4-8  shows  one  intersection  of  the  ^=   25 
rad/sec  and  UJ  =  30  rad/sec  curves.   For  this  value  of  alpha  and  beta,  a 
magnitude  of  12  decibels  is  guaranteed  to  exist  at  both  values  of  omega 
indicated.   This  observation  is  important,  for  it  means  that  the  magni- 
tude is  either  greater  than  12  decibels,  or  less  than  that  value  for  the 
intermediate  range  of  omega.   Furthermore,  it  is  possible  to  assume  that 
the  magnitude  cannot  be  both  greater  than  and  less  than  12  decibels  for 
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these  intervening  values  of  omega.   The  reason  for  this  is  simply  that 
the  magnitude  curve  is  continuous,  and  must  pass  through  12  decibels  in 
order  to  change  from  a  high  value  to  a  lower  value,  or  vica  versa.  If  a 
sufficiently  large  number  of  omega  curves  are  drawn  in  the  region  of 
interest,  it  can  be  seen  that,  at  Point  B,  only  two  omega  curves  do  actual- 
ly intersect.   Therefore,  this  range  of  omega  values  represents  the  12 
decibel  bandwidth  of  a  system  which  satisfies  these  values  of  alpha  and 
beta.   With  this  insight,  it  can  be  concluded  that  a  larger  bandwidth  is 
obtained  for  a  point  slightly  below  point  B.   Thus  at  a  value  of  alpha 
and  beta  indicated  by  point  C,  interpolation  of  the  omega  curves  shows 
an  approximate  bandwidth  from  LO  =  24  to  iO     =31.   Conversely,  proced- 
ing  in  the  opposite  direction,  the  bandwidth  is  reduced  to  the  point  D 
at  which  only  one  frequency  curve  satisfies  the  specified  value  of  alpha 
and  beta.   No  other  frequency  curves  will  intersect  this  point,  for  it 
lies  on  the  tangent  locus  of  all  frequency  curves  evaluated  at  12  deci- 
bels magnitude.   Thus  this  tangent  locus  defines  all  the  values  of  alpha 
and  beta  which  guarantee  a  peak  resonant  magnitude  of  12  decibels.   The 
value  of  the  resonant  peak  frequency  is  not  specified  by  the  tangent  locus. 
The  reason  for  this  is  because  omega  is  an  implicit  function  of  the  tan- 
gent locus  equation.   Article  2.5  discussed  various  methods  to  generate 
the  required  non-parametric  form  of  this  l^\  o    Luj)      equation,  and  the 
limitations  encountered  in  each  of  those  approaches.   For  those  reasons, 
the  tangent  locus  curve  was  not  constructed.   Instead,  a  sufficiently 
large  number  of  omega  curves  were  constructed  to  show  the  approximate  loca- 
tion of  the  tangent  locus.   Although  the  resonant  peak  frequency  is  not 
directly  obtainable  from  the  tangent  locus,  it  is  clearly  seen  that  such 
information  is  readily  found  by  the  intersection  of  a  particular  omega 
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curve  with  this  tangent  line. 

The  discussion  thus  far  has  presumed  that  all  values  inside  the 
tangent  locus  represent  peak  resonant  points  greater  than  12  decibels. 
This  is  intuitively  obvious,  but  no  proof  has  been  made  of  this  fact. 
Figures  4-9  through  4-12  have  been  included  to  verify  this  presumption, 
and  also  to  illustrate  additional  phenomena.   Figure  4-9  clearly  shows 
that  when  the  magnitude  is  reduced  to  6  decibels,  a  greater  area  is  en- 
closed by  each  of  the  omega  curves.   This  results  in  a  tangent  line  rep- 
resenting the  locus  of  peak  resonant  magnitude  equal  to  6  decibels. 
Comparison  of  this  tangent  locus  with  the  one  in  Figure  4-8  will  reveal 
that  the  12  decibel  tangent  locus  falls  within  the  area  bounded  by  the 
omega  curves  of  Figure  4-9.   The  points  B  and  D  of  Figure  4-8  have  been 
reproduced  on  the  constant  6  decibel  graph  for  comparison.   It  can  be 
noted  that  point  B  represents  a  12  decibel  bandwidth  of  5  rad/sec  and  a 
6  decibel  bandwidth  of  approximately  13  rad/sec.   Additionally,  point  D 
depicts  a  resonant  peak  magnitude  of  12  decibels  at  28  rad/sec  with  a  6 
decibel  bandwidth  of  13  rad/sec.   Similar  information  can  be  obtained 
for  other  values  of  alpha  and  beta  throughout  the  parameter  plane.   By 
similar  analogy,  it  is  possible  to  determine  what  range  of  alpha  -  beta 
values  will  guarantee  a  peak  resonance  not  to  exceed  a  pre-selected  value. 
Thus  all  values  between  the   M  f  (w)  =  12  and  the   M  p  (w)   ■  6  curves 
will  guarantee  all  magnitudes  to  be  less  than  12  decibels,  with  a  6 
decibel  bandwidth  of  variable  size  depending  on  the  choice  of  alpha  and 
beta.   A  closer  inspection  of  Figure  4-9  will  reveal  that  the  maximum  6 
decibel  bandwidth  attainable  is  approximately  17  rad/sec  for  the  situation 
just  described.   This  maximum  bandwidth  occurs  approximately  alpha  equal 
to  0.0001  and  beta  equal  to  0.018.   At  this  point  the  resonant  peak  magni- 
tude is  slightly  less  than  12  decibels,  with  a  resonant  angular  frequency 
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of  approximately  40  rad/sec.   The  cut-off  frequencies  occur  at  UJ  ■  32 
rad/sec  and  CO     =  49  rad/sec.   All  of  the  above  information  has  been 
derived  from  the  evaluation  of  two  PARAM-1  graphical  plots.   It  will  be 
recalled  that  the  PARAM-7  curves  of  Figures  4-4  through  4-7  merely  por- 
trayed Bode  curves  as  a  function  of  one  parameter  only.   These  PARAM-1 
curves  denote  a  portion  of  many  Bode  curves  as  a  function  of  two  system 
parameters.   The  value  of  these  two  is  restricted  to  only  two  magnitudes, 
but  certainly  the  analysis  of  any  two  such  graphs  provides  a  certain  in- 
sight as  to  the  system  behavior  for  other  magnitudes.   Interpolation  can 
be  used  for  the  intermediate  magnitude  range.   Such  a  procedure  is  not 
recommended  without  some  caution  for  extensions  of  magnitude  studies  out- 
side the  graphed  limits.   The  consideration  of  Figures  4-10  through  4-12 
will  illustrate  this  point  more  clearly.   For  a  constant  magnitude  of  3 
decibels,  Figure  4-10  shows  a  dramatic  increase  in  the  applicable  range 
of  alpha  -  beta  values.   The  graph  scales  of  all  PARAM-1  plots  were  kept 
the  same  simply  to  illustrate  the  large  increase  in  area  encompassed  by 
each  constant  frequency  curve  as  the  value  of  magnitude  is  decreased. 
As  might  be  anticipated,  a  much  larger  3  decibel  bandwidth  exists  for 
the  points,  B,  C  and  D  shown  on  Figure  4-8.   The  bandwidth  of  these  three 
points  do  not  vary  greatly,  however.   Each  is  approximately  20  rad/sec. 
This  illustrates  the  sensitivity  of  beta,  the  tachometer  feedback  in  the 
area  under  discussion.   For  a  given  value  of  alpha  equal  to  0.008,  a 
slight  change  in  beta  greatly  affects  the  size  of  the  peak  reasonant 
magnitude.   Conversely,  an  equal  change  in  alpha  reflects  a  vast  difference 
in  omega,  the  location  of  this  peak  resonant  magnitude,  plus  a  slight 
change  in  its  magnitude  value.   Although  the  sensitivity  of  the  alpha  vari- 
able might  not  seem  apparent  by  looking  at  the  curves,  it  must  be  remembered 
that  there  is  a  different  graph  scale  for  each  of  the  axis  variables. 
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Figures  4-11  and  4-12  show  most  clearly  why  the  analysis  of  merely 
two  PARAM-1  frequency  response  representations  might  be  misleading. 
Consider  point  B  of  Figure  4-8  which  is  reproduced  on  the  zero  decibel 
magnitude  graph  of  Figure  4-11.   The  constant  omega  curves  of  this  plot 
show  that  no  two  curves  intersect  at  point  B.   Clearly,  it  is  seen  that 
all  values  of  omega  less  than  40  rad/sec  encompass  this  choice  of  alpha 
and  beta.   Therefore  a  cut-off  frequency  of  approximately  39  rad/sec  is 
obtained  for  the  zero  decibel  bandwidth  of  that  system.   Similarly, 
from  Figure  4-12,  the  -3  decibel  bandwidth  is  about  43  rad/sec.   Using 
only  these  latter  two  PARAM-1  plots,  one  cannot  discern  any  resonant  peak 
for  the  system  frequency  response  associated  with  these  values  of  alpha 
and  beta.   It  is  only  through  the  apriori  information  portrayed  in  any 
one  of  the  Figures  4  8  through  4-10  by  which  we  are  able  to  make  such  a 
statement  that  a  resonant  magnitude  actually  exists.   Only  by  the  combined 
interpolation  of  all  the  Figures  4-8  through  4-10,  can  we  state  with  a 
fair  degree  of  confidence  that  the  value  of  this  resonant  magnitude  is 
approximately  15  decibels.   This  of  course  could  be  checked  by  construct- 
ing another  PARAM-1  curve,  or  perhaps  a  Bode  plot.   Neither  of  these  are 
necessary  for  initial  coarse  approximations.   It  will  soon  be  shown  that 
the  combined  usage  of  PARAM-1  and  PARAM-2  provides  a  much  better  method. 
The  reason  for  the  inclusion  of  so  many  PARAM-1  curves  was  simply  to  il- 
lustrate how  a  valuable  insight  can  be  obtained  for  system  frequency  re- 
sponse studies  when  considering  two  parameter  variations  simultaneously. 
Within  the  limitations  just  indicated,  it  is  possible  to  determine  how 
two  parameter  variations  affect  bandwidth  and  peak  resonance  from  a  care- 
ful inspection  of  several  PARAM-1  graphs. 

Before  proceding  to  PARAM-2  graphical  configurations,  it  will  be 
advantageous  to  review  Figures  4-8  through  4-12  in  a  three  dimensional 
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concept.   Consider  for  a  moment  a  curved  conical  surface  similar  in 
shape  to  a  bull's  horn.   Also,  let  this  curved  cone  be  elliptical  in 
shape,  rather  than  round.   Next  assume  that  the  base  of  this  surface  is 
in  the  lower  right  hand  corner,  nearest  to  the  observer.   The  smaller  end 
of  this  surface  would  therefore  be  in  the  upper  left  hand  corner  of  the 
observation  sector,  at  a  certain  distance  away  from  the  observer.   Now 
superimpose  the  omega  axis  on  the  centerline  of  this  conical  surface  in 
such  a  manner  that  the  low  values  of  omega  are  at  the  base  of  this  sur- 
face, or  nearest  to  the  observation  point.   The  higher  values  of  omega 
would  consequently  lie  further  away  from  this  observation  point,  in- 
creasing in  value  toward  the  curved  conical  tip.   With  this  three  dimen- 
sional surface  concept  in  mind,  look  at  Figure  4-8.   The  curved  omega 
axis  originates  in  the  lower  right  hand  corner  and  extends  through  the 
alpha  -  beta  plane,  curving  upward  and  slightly  to  the  right.   Figure  4-9 
depicts  this  curved  conical  surface  as  larger  in  size,  (i.e.,  closer  to 
the  observer),  with  the  base  of  this  surface  shifted  slightly  to  the  left, 
thereby  bringing  it  more  into  the  line  of  sight.   Bringing  this  curved 
conical  surface  still  closer  to  the  observation  point  makes  it  appear 
larger,  as  seen  in  Figure  4-10.   The  zero  decibel  PARAM-1  graph,  Figure 
4-11,  portrays  an  even  larger  shape  with  the  base  shifted  even  more  to 
the  left.   Finally  Figure  4-12  can  be  regarded  observing  this  curved  sur- 
face with  its  base  very  close,  similar  in  appearance  to  passing  through 
a  curved  tunnel  under  a  river. 

The  preceding  illustrative  example  can  be  considered  as  elementary, 
but  is  fundamental  to  the  understanding  of  the  various  graphical  rep- 
resentations to  follow.   No  attempt  will  be  made  to  correlate  the  il- 
lustration to  the  PARAM-7  curves  of  Figures  4-4  through  4-7,  since  those 
are  essentially  single  parameter  variations  of  the  series  of  two  parameter 
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surfaces  just  described.   Better  use  of  the  illustrative  example  is  made 
by  correlating  PARAM-1  with  PARAM-2. 

PARAM-2. 

Consider  now  a  cross  sectional  cut  of  the  curved  conical  surface 
described  above.   PARAM-2  effectively  maps  constant  magnitude  curves  onto 
the  alpha  -  beta  parameter  plane  for  a  fixed  value  of  angular  frequency. 
In  this  manner,  it  is  possible  to  show  how  a  particular  omega  value, 
such  as  UJ  =  10   rad/sec  will  vary  as  a  function  of  magnitude  in  the 
alpha  -  beta  plane.   Consider  again  the  (xv  =  10  rad/sec  curve  in  Figures 
4-8  through  4-11.   The  associated  PARAM-2  plot  is  shown  in  Figure  4-13. 
Since  the   W  =  10  rad/sec  curve  represented  a  cross  section  of  the  curved 
conical  surface  near  its  base,  it  is  possible  to  observe  how  this  base 
is  rotated  toward  the  left  as  magnitude  is  decreased.   Further  inspec- 
tion of  Figure  4-13  shows  how  the  curve  increases  in  size  with  a  reduc- 
tion of  magnitude.   The  -3  decibel,  uJ    =  10  rad/sec  curve  is  in  fact 
shifted  completely  off  the  graph  of  Figure  4-13.   This  effectively  shows 
that  no  alpha  -  beta  value  within  the  prescribed  range  will  have  a  magni- 
tude of  -3  decibels  at  an  angular  frequency  of  10  rad/sec.   No  information 
is  conveyed  about  any  other  angular  frequency  in  Figure  4-13.   To  obtain 
such  information  would  necessitate  a  different  PARAM-2  plot,  or  a  compari- 
son with  one  or  more  PARAM-1  curves.   This  does  not  preclude  the  correla- 
tion of  other  frequency  response  curves.   They  can  be  consulted  for  pert- 
inent relationships,  and  will  be  discussed  later.   Figure  4-14  is  a  graphi- 
cal representation  of  constant  magnitude  curves  for  a  fixed  angular  fre- 
quency of  20  rad/sec.   Inspection  of  these  curves  shows  that  they  lie  very 
closely  to  the  centroid  of  the  curved  conical  surface  discussed  in  the  il- 
lustrative example  of  this  article.   This  can  be  verified  by  their 
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concentric  shape.   Due  to  the  graphical  construction  technique,  these 
curves  do  not  extend  to  the  bottom  of  the  graph.   This  can  be  explained 
by  the  fact  that  these  curves  are  almost  vertical  at  the  lower  graph 
edge.   Even  though  the  abcissa  variable  is  incremented  900  times  during 
each  curve  computation,  this  is  not  sensitive  enough  to  allow  the  detec- 
tion of  such  a  steep  slope.   If  an  infinitessimally  small  step  size  were 
possible,  these  curve  values  would  be  detected.   This  limitation  of  the 
graphical  plotting  procedure  must  be  kept  in  mind.   The  mathematical 
theory  of  section  2  indicated  that  two  quadratic  formula  solutions  were 
necessary  for  the  complete  solution  of  any  curve  generated  by  PARAM-1 
through  PARAM-6.   The  computer  technique  discussion  of  section  3  dis- 
cussed how  each  graph  curve  was  constructed  first  for  one  sign  of  the 
square  root  coefficient  in  the  quadratic  formula,  and  then  for  the  other. 
Subsequent  graphs  will  show  how  these  two  curve  solutions  do  not  always 
connect.   The  reason  for  this  again  lies  in  the  finite  incrementation  of 
the  abcissa  variable.   Returning  once  again  to  Figure  4-14,  it  is  possible 
to  note  certain  characteristics  regarding  the  angular  frequency  of  20  rad/ 
sec  as  a  function  of  the  two  parameters  alpha  and  beta.   The  intersection 
of  the  peak  values  of  each  magnitude  curve  defines  a  vertical  line  in  the 
parameter  plane.   Therefore  it  is  possible  to  conclude  that  variations  in 
magnitude  for  an  angular  frequency  of  20  rad/sec  are  highly  sensitive  to 
beta  values,  and  quite  insensitive  to  changes  in  alpha.   This  says  nothing 
about  resonance  magnitudes.   Rather,  each  curve  guarantees  a  particular 
magnitude  value  for  a  given  angular  frequency.   Also,  each  curve  defines 
an  alpha  -  beta  region  within  which  the  magnitude  will  always  be  greater 
than  the  boundary  value  for  that  value  of  omega.   As  a  further  analysis  of 
the  system  response  for  point  B  of  Figure  4-8,  it  can  be  seen  that  such  a 
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system  will  have  a  magnitude  of  6  decibels  at  a  frequency  of  20  rad/sec. 
Figures  4-15  and  4-16  show  considerably  more  information  than  the  previous 
two  PARAM-2  graphs,  because  the  center  of  these  eliptical  curves  is  indi- 
cated.  This  center  corresponds  to  the  centerline  of  the  three  dimensional 
curved  conical  surface  which  has  been  repeatedly  referenced  up  to  this 
point.   It  should  be  observed  that  the  center  of  these  ellipses  in  Figures 
4-15  and  4-16  denote  the  maximum  magnitudes  attainable  at  their  respec- 
tive angular  frequencies.   Thus  by  analogy,  the  centerline  of  the  three 
dimensional  surfaces  in  Figures  4-8  through  4-16  serves  a  dual  purpose. 
Not  only  does  it  specify  a  particular  angular  frequency,  but  it  denotes 
the  locus  of  the  highest  peak  resonant  magnitudes  for  any  alpha  -  beta 
combination  of  the  system!   This  is  an  important  conclusion,  because  it 
indicates  whether  a  desired  peak  resonant  magnitude  value  can  ever  be 
reached  for  the  system  under  consideration.   The  comparison  of  PARAM-1 
and  PARAM-2  graphs  further  shows  how  the  conical  surfaces  expand  about 
this  centerline  with  decreasing  magnitude.   From  these,  it  is  possible 
to  visualize  the  locus  of  constant  magnitude  conical  surfaces,  each 
surrounding  the  centerline  just  discussed.   It  is  therefore  seen  that  a 
reduction  of  tachometer  gain  coupled  with  an  increase  in  acceleration 
gain  will  shift  the  resonance  peak  to  a  higher  frequency.   Point  B  of 
Figure  4-8  has  also  been  reproduced  on  Figures  4-15  and  4-16.   These 
show  that  a  system  with  the  indicated  alpha  -  beta  values  will  have  a 
nagnitude  of  12  decibels  at  uJ  -   30  rad/sec,  and  -1  decibels  at  an  angular 
frequency  of  40  rad/sec.   At  this  point,  enough  information  has  been  ob- 
tained to  obtain  a  fairly  accurate  Bode  plot  for  the  parameter  values  given 
at  point  B.   Further  comparison  of  the  15  decibel  magnitude  curve  in 
Figures  4-14  through  4-16  show  that  point  B  does  in  fact  have  a  magnitude 
equal  to  15  decibels.   This  confirms  the  assumption  of  this  value  based  on 
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interpolation  of  the  PARAM-1  curves  alone.   Since  enough  information  has 
been  depicted  in  the  PARAM-1  and  PARAM-2  curves  to  portray  a  Bode  plot 
of  point  B,  it  is  obvious  that  similar  Bode  plots  could  be  obtained  for 
any  other  desired  combination  of  alpha-beta  points.   This  is  exceedingly 
valuable,  because  it  allows  the  design  engineer  to  construct  any  number 
of  coarse  Bode  plots  by  evaluating  the  data  contained  in  Figures  4-8 
through  4-16.   Therefore  it  is  concluded  that  the  combined  utilization 
of  PARAM-1  and  PARAM-2  is  a  fundamental  starting  point  for  all  two  para- 
meter frequency  response  design  studies. 

The  preceding  discussion  concerning  the  use  of  PARAM-1  and  PARAM-2 
was  very  comprehensive  because  of  the  following  three  key  factors:  First, 
a  thorough  understanding  of  new  concepts  can  only  be  gained  by  a  slow  and 
detailed  presentation  of  the  uses  and  limitations  of  such  innovations. 
Secondly,  the  combined  application  of  PARAM-1  and  PARAM-2  plots  provides 
the  broadest  information  spectrum  for  frequency  response  design  studies. 
Finally,  a  thorough  understanding  of  these  two  parameter  plots  is  pre- 
requisite to  visualizing  the  following  frequency  response  projections 
which  are  essentially  a  function  of  one  system  parameter.   The  remainder 
of  this  article  will  cover  the  various  roles  played  by  PARAM-3  through 
PARAM-6. 

PARAM-3 . 

Section  3  developed  the  major  graphical  presentations  which  were 
considered  to  be  of  most  use  to  design  engineer.   In  that  discussion, 
it  was  indicated  that  two  basically  different  configurations  were  avail- 
able in  PARAM-3.   Figure  3-1  showed  in  schematic  form  these  differences. 
Attention  will  first  be  focused  on  the  M0D0RD=1   variety  as  shown  in 
figure  4-17.   Before  discussing  the  information  content  of  this  graph, 
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an  apparent  mistake  in  the  graphical  plotting  procedure  is  noteworthy. 
This  set  of  curves  was  purposely  offset  to  illustrate  a  limitation  en- 
countered by  careless  selection  of  axis  variable  ranges.   Graph  scales 
are  limited  to  one  significant  figure.   As  such  the  automatic  scaling  of 
the  program  PARAMS  does  not  always  produce  the  values  which  might  be  ex- 
pected.  The  requested  range  of  alpha  values  was  from  -0.001  to  0.007. 
Since  the  ordinate  length  was  specified  as  ten  inches,  the  program  comput- 
ed a  grid  scale  of  0.0008  units  per  inch  for  alpha.  Accordingly  this 
scale  was  used;  however,  it  placed  the  origin  one  inch  above  the  minimum 
value  in  accordance  with  program  plotting  limitations.   Thus  the  requested 
range  of  alpha  values  was  plotted  properly,  even  though  the  choice  of 
ordinate  values  was  improperly  chosen.   The  preceding  was  for  information 
purposes  only;  PARAM-3  plots  will  now  be  studied. 

Figure  4-17  portrays  a  set  of  constant  beta  curves  on  the  magnitude 
-  alpha  plane  for  a  given  fixed  angular  frequency  value  of  uO  =  30  rad/ 
sec.   This  graph  can  best  be  understood  by  comparing  it  with  the  PARAM-2 
plot  of  the  same  constant  omega,  namely  Figure  4-15.   Effectively,  the 
PARAM-3,  M0D0RD=1  curves  take  the  information  contained  in  the  horizontal 
cross -sections  of  PARAM-2  and  map  this  data  onto  the  magnitude  -  alpha 
plane.   This  can  be  most  easily  visualized  by  taking  the  base  of  Figure 
4-15  and  comparing  it  with  the  beta  =0.0  curve  in  Figure  4-17.   With 
these  two  viewpoints,  one  readily  sees  another  facet  of  PARAM-2  plots. 
Consider  for  a  moment  a  topographical  map  of  land  terrain.   Such  maps 
represent  constant  altitude  as  a  continuous  curve.   Each  such  curve  is 
always  a  closed  contour,  with  either  peaks  or  valleys  contained  within 
the  curve  boundaries.   PARAM-2  is  such  a  topographical  plot,  with  magni- 
tude equal  to  elevation.   This  concept  is  borne  out  by  the  magnitude- 
alpha  plane  of  Figure  4-17.   The  constant  beta  curves  of  this  latter 
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graph  can  be  regarded  as  vertical  slices  into  the  topographical  hill- 
side. As  such,  one  sees  which  cross  sections  have  actually  been  sliced 
out  of  the  PARAM-2  magnitude  contours.   The  complete  range  of  PARAM-2 
alpha  and  beta  values  has  been  spanned  by  the  ordinate  range  and  constant 
beta  curves  in  Figure  4-17.   Accordingly  one  might  conclude  that  a  good 
indication  of  the  system's  magnitude  behavior  has  been  conveyed  by  the 
PARAM-3,  M0D0RD=1  graph  for  u)   =  30  rad/sec.   This  is  most  emphatically 
not  true!  Analysis  of  Figure  4-17  indicates  that  the  greatest  magnitude 
occurs  at  beta  =0.0  and  that  the  magnitude  contours  decrease  as  beta  ap- 
proaches 0.05.   Careful  review  of  Figure  4-15  will  disclose  a  serious 
fallacy  in  such  a  conclusion.   Figure  4-17  shows  no  indication  of  the 
magnitude  peak  resonance  at  beta  =  0.004.   Without  this  apriori  knowledge, 
a  completely  erroneous  interpretation  would  have  been  made  by  analyzing 
the  curves  in  PARAM-3,  M0D0RD=1.   This  example  stresses  once  again  the 
importance  of  viewing  system  frequency  response  as  a  function  of  the  two 
system  parameters. 

Consideration  will  now  be  given  to  the  alternate  PARAM-3  curve  desig- 
nated by  M0D0RD=2.   This  graph  presents  constant  alpha  curves  in  the  beta 
magnitude  plane  for  a  fixed  omega.   As  such,  these  curves  simply  slice 
vertical  cross-sectional  areas  out  of  the  magnitude  contours  of  PARAM-2. 
Figure  4-18  shows  such  a  set  of  curves  for  a  constant  angular  frequency  of 
30  rad/sec.   Looking  at  Figure  4-15  again,  all  curve  values  lie  within  a 
fairly  small  range  of  alpha  values.   Accordingly,  one  would  only  expect 
graph  curves  for  alpha  less  than  0.0021  to  lie  in  the  beta  -  magnitude 
plane  for  magnitudes  greater  than  or  equal  to  -3  decibels.   Furthermore 
one  would  expect  a  series  of  magnitude  peaks  to  occur  at  beta  ■  .004  with 
the  maximum  amgnitude  occuring  on  the  alpha  ■  0.0005  curve.   These  factors 
are  verified  in  Figure  4-18  where  the  alpha  =  0.0005  curve  has  maximum 
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magnitude  greatly  in  excess  of  21  decibels.   All  other  curves  represent 
lower  magnitude  ranges  as  was  anticipated.   In  this  set  of  alpha  curves, 
the  maximum  magnitude  curve  was  intentionally  constructed  to  indicate 
the  large  magnitude  at  peak  resonance.   This  curve  might  have  been  over- 
looked, were  it  not  for  the  apriori  knowledge  gained  from  Figure  4-15. 

The  preceding  discussion  of  PARAM-3  curves  showed  serious  limita- 
tions in  their  usage  when  the  PARAM-2  curve  of  the  same  constant  omega  is 
not  utilized.   When  used  in  conjunction  with  PARAM-2  graphs,  however,  the 
PARAM-3  curves  do  provide  some  useful  information.   They  are  particularly 
useful  for  the  area  enclosed  by  the  largest  value  of  magnitude  in  the 
PARAM-2  plots.   This  area  tells  nothing  about  the  peak  resonant  magnitude, 
other  than  the  frequency  at  which  it  occurs.   By  constructing  one  PARAM-3 
curve  along  the  bisector  of  these  magnitude  curves,  one  quickly  gains  the 
value  of  peak  magnitude  and  the  magnitude  sensitivity  about  this  resonant 
point.   Should  the  precise  peak  magnitude  location  not  be  so  obvious  on 
other  PARAM-2  plots,  a  series  of  curves  in  the  PARAM-3  graph  will  clearly 
indicate  its  location.   Due  to  the  elliptical  shape  of  the  PARAM-2  magni- 
tude curves  for  this  system,  the  PARAM-3,  M0D0RD=2  configuration  would 
be  better  suited  than  the  PARAM-3,  M0D0RD=1  graph.   The  reason  for  this 
is  simply  that  the  magnitude  is  more  sensitive  to  alpha  than  beta.   There- 
fore a  much  smaller  range  of  alpha  values  need  be  considered  than  would  be 
required  for  beta  curves. 

Up  to  this  point,  only  one  value  of  omega  has  been  discussed.   The 
PARAM-3  curves  for  other  angular  frequencies  can  be  visualized  by  studying 
Figures  4-13,  4-14,  and  4-16.   The  PARAM-2  plot  of  Figure  4-13  shows  a 
very  gentle  magnitude  slope  along  constant  alpha  and  constant  beta  lines. 
Thus  for  constant  alpha  lines,  one  would  expect  the  magnitude  to  decrease 
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slowly  with  an  increase  in  beta.   The  magnitude  would  decrease  more 
rapidly  for  larger  constant  alpha  curve  values  than  they  would  for  smaller 
curve  values.   Conversely,  for  constant  beta  lines,  it  is  seen  that  magni- 
tude will  increase  slowly  as  alpha  is  increased.   Large  constant  beta 
curve  values  would  show  a  slower  increase  in  magnitude  than  small  beta 
curve  values.   Figures  4-19  and  4-20  are  included  for  verification  of 
these  trends. 

Inspection  of  Figure  4-14  shows  that  the  PARAM-3,  M0D0R=1  graph  will 
have  a  maximum  magnitude  at  alpha  =  0.0018  for  all  beta  curves.   The  maxi- 
mum value  of  magnitude  for  any\alue  of  alpha  will  be  obtained  on  the  beta 
=  0  curve.   All  other  beta  curves  will  lie  inside  this  magnitude  boundary. 
The  shpae  of  these  curves  will  be  similar  to  those  shown  on  Figure  4-18. 
Since  the  magnitude  is  less  sensitive  at  20  rad/sec  than  at  30  rad/sec,  the 
curves  will  be  less  pointed.   The  PARAM-3,  M0D0RD-2  curves  will  not  be  of 
the  shape  just  discussed.   This  is  because  the  vertical  cross-sections  of 
the  PARAM-2  plot  show  a  decrease  in  magnitude  as  beta  increases.  All  alpha 
curves  will  lie  inside  the  constant  alpha  curve  equal  to  0.0018,  for  this 
curve  represents  a  ridge  of  maximum  magnitude.   The  constant  CO  =  40  rad/ 
sec  magnitude  curves  of  Figure  4-16  show  that  the  peak  resonance  occurs  on 
the  beta  axis  at  an  alpha  value  of  approximately  0.012.   The  sensitivity 
of  magnitude  is  greatest  in  this  graph.  Accordingly,  the  PARAM-3,  M0D0RD=1 
curves  will  be  much  more  pointed  than  those  shown  in  Figure  4-17.   Also 
the  center  of  these  magnitude  peaks  will  be  shifted  down  to  the  abcissa. 
The  PARAM-3,  M0D0RD=2  curves  will  be  likewise  sharper  than  those  shown  in 
Figure  4-18,  but  this  time  the  centerline  of  the  magnitude  peaks  will  lie 
on  the  beta  =  0.012  line. 


108 


It  has  been  shown  that  for  this  system  the  PARAM-3  curves  are  of 
limited  use.   Detailed  correlation  between  PARAM-2  and  PARAM-3  plots  was 
made  to  show  the  inter-relationships  which  exist.   The  reason  for  discus- 
sing PARAM-3  graphs  was  to  show  how  such  information  can  be  picked  off  of 
the  PARAM-2  plots,  thereby  eliminating  the  requirement  for  another  graphi- 
cal presentation.   A  more  important  reason  for  this  detailed  analysis  was 
to  point  out  the  correlations  applicable  to  this  system  under  study.   Many 
systems  will  undoubtedly  be  more  complex  than  this.   Such  cases  might  re- 
quire the  use  of  these  curves.   Therefore  it  is  essential  to  obtain  a  firm 
understanding  of  these  graphical  advantages  and  limitations. 

PARAM-4. 

The  graphical  concepts  of  PARAM-4  will  now  be  discussed  for  the  third 
order  system  under  consideration.   Again,  two  basically  different  plots  are 
obtained  for  this  program.   The  PARAM-4,  M0D0RD=1  configuration  graphs 
constant  magnitude  curves  onto  the  alpha  -  omega  plane  for  a  fixed  value 
of  beta.   Figure  4-21  shows  one  such  plot  for  a  given  fixed  beta  value  of 
0.036.   At  first  glance,  these  curves  seem  very  strange.   Once  again,  how- 
ever, recourse  to  the  two  parameter  frequency  response  curves  of  PARAM-1 
will  show  a  striking  relationship.   Recall  for  a  moment  the  three  dimen- 
sional concept  that  was  introduced  in  the  study  of  PARAM-1  graphs.   Each 
three  dimensional  curved  elliptical  cone  shaped  surface  was  a  function  of 
alpha,  beta,  and  omega  for  a  given  value  of  magnitude.   As  magnitude  was 
decreased,  the  conical  surface  base  was  observed  to  rotate  into  the  obser- 
vation line  of  sight,  while  the  overall  size  expanded.   The  simultaneous 
rotation  and  expansion  caused  the  base  portion  to  expand  more  than  the 
smaller  end.   Furthermore,  recall  that  omega  increased  in  value  from  the 
base  toward  the  small  end.   The  PARAM-4,  M0D0RD=1  graph  effectively  slices 
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a  horizontal  cross  section  out  of  this  PARAM-1  three  dimensional  sur- 
face at  the  prescribed  value  of  beta.   To  see  this  more  clearly,  Figures 
4-8  through  4-12  will  be  reviewed  in  order.   Consider  first  the  constant 
decibel  magnitude  PARAM-1  graph  of  Figure  4-8.   By  cutting  a  perpendicu- 
lar cross  section  along  the  beta  =  0.036  line,  it  is  seen  that  only  one 
omega  curve  is  severed  namely  the  UJ   =  65   rad/sec  curve.   Note,  however, 
that  omega  is  a  continuous  elliptical  surface;  only  a  small  number  of 
cross  sections  are  depicted  in  Figure  4-8.   By  interpolation,  it  can  be 
seen  that  a  perpendicular  cross  section  of  PARAM-1  will  intersect  the 
range  of  omega  values  between  approximately  62  through  67.   Furthermore 
note  that  no  other  omega  curves  have  been  severed  in  this  manner.   This 
intersection  occurs  at  a  range  of  alpha  values  between  -0.002  and  -0.004. 
Now  turn  to  Figure  4-21  and  inspect  the  12  decibel  magnitude  curves  re- 
produced there.   By  a  similar  analysis  of  Figures  4-9  through  4-12,  it  is 
possible  to  establish  a  positive  correlation  between  the  PARAM-1  and  PARAM- 
4,  M0D0RD=1  graphs.   Next  consider  the  PARAM-4,  M0D0RD=2  plot  shown  in 
Figure  4-22.   This  configuration  graphs  constant  magnitude  curves  onto 
the  beta  -  omega  plane  for  a  fixed  value  of  alpha.   By  methods  analogous 
to  those  presented  above,  it  is  possible  to  slice  a  vertical  cross  section 
of  the  PARAM-1  graph  at  the  prescribed  value  of  alpha.   The  information 
thus  obtained  is  seen  to  appear  as  PARAM-4,  M0D0RD=2.   The  PARAM-4  curves 
are  effectively  topographical  maps  with  magnitude  equivalent  to  elevation. 
Their  utility  lies  in  the  ability  to  correlate  magnitude,  frequency  and  one 
system  parameter.   They  are  restricted  by  the  second  system  parameter  which 
has  been  held  at  a  constant  value.  As  such  they  would  only  be  suited  for 
supplying  additional  frequency  response  information  after  an  initial  choice 
of  one  variable  has  been  made.   Providing  a  sufficient  number  of  magnitude 
curves  have  been  constructed  on  these  graphs,  they  become  in  essence  a 
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locus  of  Bode  plots.   When  used  in  this  manner,  they  become  valuable 
tools.   The  constant  alpha  and  beta  values  chosen  for  figures  4-21  and 
4-22  were  so  chosen  because  they  coincided  with  the  values  used  in 
figures  4-4  through  4-7,   Thus  a  comparison  of  these  PARAM-4  curves  with 
the  familiar  Bode  magnitude  plots  can  easily  be  made.   The  PARAM-4 
curves  can  therefore  be  regarded  as  most  useful  for  coarse  Bode  plot 
studies,  since  they  portray  magnitude -frequency  information  for  a  much 
larger  number  of  alpha  or  beta  values  than  would  be  feasible  in  a  con- 
ventional Bode  plot  study.   In  this  manner  a  design  engineer  would  most 
likely  construct  enough  PARAM-1  and  PARAM-2  plots  to  make  a  judicuous 
first  approximation  for  each  of  the  system  parameters.   Using  these, 
he  could  then  construct  both  PARAM-4  plots  and  check  his  system  perfor- 
mance with  a  fair  degree  of  accuracy.  After  suitable  corrections,  he 
would  then  use  PARAM-7  for  a  final  refined  system  frequency  response 
performance  check. 

PARAM-5 


The  discussions  of  PARAM-3  and  PARAM-4  have  considered  cross  section- 
al cuts  of  the  PARAM-2  and  PARAM-1  planes  respectively.   A  slightly  dif- 
ferent concept  is  portrayed  in  PARAM-5  and  PARAM-6,   Consider  first  the 
PARAM-5,  M0D0RD=1  graph  shown  in  figure  4-23.   This  shows  a  series  of 
omega  curves  plotted  upon  the  alpha  magnitude  plane  for  a  fixed  value 
of  beta  =  0.025.   Each  constant  magnitude  depicts  a  horizontal  cross 
section  of  the  PARAM-1  graph  corresponding  to  that  constant  magnitude. 
The  location  of  the  horiztonal  cross-section  is  fixed  by  the  constant 
beta  =  0.025  for  the  PARAM-5  plot.   This  is  most  clearly  shown  by  the 
following  comparison.   Arbitrarily  choose  the  constant  magnitude  line  of 
6  decibels.   This  automatically  selectes  the  PARAM-1  plot  shown  in 
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figure  4-9.   Since  beta  is  fixed  at  a  value  of  0.025,  this  is  the 
horizontal  line  of  figure  4-9  to  be  considered.   Positive  correlation 
between  the  omega  curves  crossing  this  beta  ■  0.025  line  in  figure  4.9 
is  thereby  made  with  those  omega  curves  crossing  the  magnitude  line  of 
6  decibels  in  figure  4-23.   Effectively  the  PARAM-5,  M0D0RD  graph  shows 
the  series  of  PARAM-1  plots  for  any  given  constant  beta  value.  Merely 
studying  the  different  magnitude  lines,  one  can  see  how  a  constant  beta 
cross  section  would  look  on  different  PARAM-1  plots.   The  PARAM-5, 
M0D0RD=2  graphs  show  this  same  behavior,  with  the  roles  of  alpha  and 
beta  interchanged. 

A  comparison  can  also  be  made  with  a  series  of  PARAM-2  plots,  since 
they  represent  another  three  dimensional  surface  projected  into  the 
alpha  beta  parameter  plane.   In  that  case,  however,  each  omega  curve  in 
the  PARAM-5  graph  represents  a  different  PARAM-2  plot.   Consider  for  a 
moment  the  U>=20  rad/sec  curve  of  figure  4-23.   This  curve  depicts  the 
topological  relief  of  a  cross  section  of  the  PARAM-2  magnitude  contour 
graph  when  that  cross  section  occurs  at  beta  =  0.025.   Similarly,  dif- 
ferent omega  curves  of  figure  4-23  represent  such  a  cross  section  in 
other  PARAM-2  plots. 

The  correlation  between  PARAM-5  and  PARAM-1  as  well  as  PARAM-2  pro- 
vides a  rather  important  conclusion.   One  normally  wants  to  get  a  good 
indication  of  system  behavior  by  simultaneously  varying  alpha  and  beta. 
It  is  not  necessary  to  construct  a  large  number  of  PARAM-1  plots  to  do 
this.   Neither  is  it  necessary  to  construct  a  large  number  of  PARAM-2 
plots.   By  constructing  only  a  few  PARAM-1  graphs,  as  well  as  two  or 
at  most  three  PARAM-2  graphs,  one  can  gain  valuable  insight  into  the 
system  behavior  in  the  voids  by  studying  the  PARAM-5  plots.   Thus  for  an 
initial  system  design  study,  one  might  be  inclined  to  use  the  following 


112 


six  graphs:   two  PARAM-1  plots,  two  FARAM-2  plots,  and  two  PARAM-5 
plots.   The  constant  magnitude  values  of  the  PARAM-1  plots  would  most 
likely  be  chosen  from  the  bandwidth  and  resonant  peak  criteria.    The 
constant  omega  values  for  the  two  FARAM-2  plots  would  normally  be  dictat- 
ed by  the  cut-off  frequencies.   Alternatly,  the  peak  resonant  frequency 
might  be  used  as  a  governing  choice.   By  choosing  an  intermediate  value 
of  alpha  and  beta,  the  PARAM-5  curves  would  indicate  how  other  PARAM-1 
and  PARAM-2  graphs  would  have  appeared.   For  illustrative  purposes, 
suppose  a  problem  called  for  a  -3  decibel  bandwidth  of  60  rad/sec,  with  a 
peak  resonance  magnitude  of  nc  more  than  6  decibels  around  40  rad/sec. 
Construction  of  figures  4-9,  4-12,  and  4-16,  one  would  be  able  to  con- 
clude the  following  information: 

1)  The  tangent  locus  of  the  omega  curves  in  figure  4-9  define 
the  limiting  values  of  alpha  and  beta  as  dictated  by  the 
resonance  magnitude  constraint. 

2)  The  intersection  of  the  tangent  locus  and  the  constant 
U)  -   40  rad/sec  curve  will  determine  which  values  of 
alpha  and  beta  will  give  this  desired  resonant  frequency. 

3)  Looking  at  figure  4-12,  one  sees  that  the  alpha  and  beta 
values  determined  above  will  not  provide  a  -3  decibel  band- 
width of  60  rad/sec.   Instead,  it  can  be  seen  by  interpola- 
tion that  the  bandwidth  will  be  approximately  57  rad/sec. 

4)  Figure  4-16  shows  the  magnitude  relationship  at  the  de- 
sired resonant  frequency.   Thus  by  studying  this  graph, 
one  sees  that  a  resonant  magnitude  of  approximately  7 
decibels  will  result  from  a  60  rad/sec  bandwidth. 

The  design  example  above  has  shown  a  common  problem  which  exists  for 
most  engineering  studies.   The  desired  performance  measures  cannot  be 

113 


attained;  instead  some  compromise  must  be  made.   This  conclusion  was 
derived  by  simply  observing  three  frequency  response  parameter  plane 
graphs.   The  advantage  of  PARAM-5  graphs  lies  in  the  ability  to  corre- 
late performance  sensitivities  with  parameter  variations.   Even  though 
the  curves  shown  in  figure  4-23  do  not  coincide  with  the  beta  most  likely 
to  be  selected,  they  do  illustrate  the  behavior  of  magnitude  and  omega 
at  a  neighboring  beta  value.   Thus  it  is  possible  to  project  that  in- 
formation onto  one  of  the  PARAM-1,  PARAM-2  graphs  and  extrapolate  the 
data  thus  obtained.   Additionally,  one  can  do  the  same  with  the  PARAM-5 
curves  shown  in  figure  4-24.   This  figure  shows  the  magnitude  and  omega 
variations  as  a  function  of  beta  for  a  constant  alpha  also  near  the 
point  of  interest.   In  effect  the  curves  in  figure  4-24  reflect  the  be- 
havior about  the  peak  resonant  point  of  UJ=35  rad/sec.   It  is  seen  that 
a  variation  in  beta  will  increase  the  resonant  magnitude  for  smaller  beta 
values,  with  a  corresponding  decrease  in  magnitude  for  higher  beta  values. 
Figure  4-24  shows  another  important  factor.   At  beta  =  0.036  one  sees  that 
a  flat  frequency  response  of  0  decibels  will  be  obtained  for  a  bandwidth 
of  approximately  35  rad/sec.  This  can  be  verified  in  figures  4-11  and 
4-25. 

Figures  4-25  and  4-26  have  been  included  to  show  more  vividly  how  the 
interpolation  of  PARAM-5  plots  might  be  made.   Consider  the  PARAM-5  graphs 
as  topological  portrayals  of  omega.   Figure  4-25  shows  how  these  omega 
contours  change  from  those  of  Figure  4-23  as  beta  increases.   Note  that 
all  omega  curves  are  reduced  in  magnitude  when  beta  is  increased.   Note 
further  that  a  reduction  in  the  number  of  omega  curves  will  result  in  a 
loss  of  important  information.   The  eight  curves  shown  in  figure  4-25 
might  lead  the  observer  to  conclude  that  magnitude  decreases  continually 
as  alpha  is  decreased.   Only  by  careful  inspection  of  the  u)  -   40  rad/sec 

114 


curve  can  one  see  that  this  is  not  the  case.   By  inspecting  the  complete 
array  of  curves  in  figure  4-23,  however,  one  sees  a  resonance  trend  for 
negative  alpha  values.   In  this  figure,  resonance  occurs  at  approximately 
55  rad/sec.   Using  this  information,  another  PARAM-2  graph  can  be  visual- 
ized for  constant  U)   ■  55  rad/sec.   The  curves  would  be  similar  in  shape 
to  those  shown  in  figure  4-15  and  4-16,  but  smaller  in  size.   The  center 
of  these  elliptical  curves  would  occur  at  beta  =  0.036  and  alpha  =  -0.003. 
The  analogous  procedure  could  also  be  employed  using  either  figures  4-24 
or  4-26. 

PARAM-6 

The  discussions  of  PARAM-3  through  PARAM-5  have  shown  how  various 
graphical  configurations  result  from  taking  cross  sectional  cuts  in 
the  PARAM-1  and  PARAM-2  frequency  response  parameter  plane  curves.   The 
curves  shown  in  PARAM-6  depict  another  such  cross  sectional  projection. 
Consider  first  the  curves  shown  in  figure  4-27.   They  correspond  to 
constant  beta  curves  mapped  onto  the  alpha-omega  plane  for  a  given  pre- 
selected magnitude  value  of  12  decibels.   Accordingly,  this  graph  corres- 
ponds to  the  PARAM-1  plot  shown  in  figure  4-8.   The  full  utilization  of 
the  entire  set  of  graphs  will  now  become  clear.   Effectively  the  PARAM-6 
graphs  show  the  third  variable  of  the  curved  elliptical  surface  which 
was  projected  onto  the  parameter  plane  in  PARAM-1.   In  this  manner,  the 
third  variable,  omega,  now  becomes  one  of  the  axis  variables  while  cross 
sectional  cuts  of  the  three  dimensional  alpha-beta-omega  surface  are  ob- 
tained by  slicing  down  along  constant  beta  lines  or  constant  alpha  lines. 
The  latter  cross  sections  are  pictured  in  figure  4-28.   Thus  when  the 
corresponding  PARAM-1  and  PARAM-6  graphs  are  viewed  simultaneously,  one 
can  visualize  the  three  dimensional  alpha-beta-omega  surface  in  its  true 
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form.   The  illustrative  example  of  an  elliptically  shape  bull's  horn  is 
applicable  in  all  but  one  facet  which  will  now  be  corrected.   It  was 
originally  stipulated  that  the  omega  axis  passed  down  the  conical  surface, 
and  was  rotated  and  shifted  slightly  as  magnitude  decreased.   This  stipula- 
tion was  made  to  simplify  the  visualization  of  the  surface.   Inspection  of 
the  PARAM-6  curves  will  show  a  slight  modification  of  this  illustration. 
The  omega  axis  does  not  pass  down  the  surface,  but  instead  projects  per- 
pendicularly through  the  origin  of  the  PARAM-1  plots.   In  this  manner  it 
is  seen  that  the  surface  has  shifted  slightly,  and  been  rotated  in  the 
three  dimensional  alpha-beta-omega  space.'  Now  the  illustration  is  cor- 
rect in  its  entirety.   Analyzing  the  corresponding  PARAM-1  and  PARAM-6 
graphs  gives  a  complete  description  of  the  portrayed  surface  in  three 
dimensions.   Inspection  of  any  one  of  the  PARAM-2  plots  reveals  the 
change  in  shape  of  this  three  dimensional  surface  as  the  fourth  dimen- 
sional parameter,  magnitude,  is  varied.   Figures  4-29  through  4-32  are 
included  for  amplification  and  verification  of  the  above  results. 
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5.    CONCLUSIONS 

Section  4  presented  an  exhaustive  study  of  the  application  of  fre- 
quency response  parameter  plane  curves  to  one  particular  system.   The 
selection  of  the  third  order  system  was  based  on  the  ease  of  interpre- 
tation of  each  frequency  response  graph  discussed.   Extensions  to  systems 
of  higher  order  will  undoubtedly  change  the  graphical  configurations  to 
some  extent.   The  degree  of  this  variation  was  smaller  than  expected  for 
a  similar  fourth  order  system.   Figure  5-1  shows  the  particular  fourth 
order  tachometer  and  acceleration  feedback  system  which  was  studied. 
Figure  5-2  presents  the  system  stability  performance  curves  in  the  alpha- 
beta  parameter  plane.   Figures  5-3  through  5-12  depict  the  associated 
frequency  response  curves.  Comparison  with  the  third  order  system  discuss- 
ed in  section  4  shows  strikingly  similar  results. 

Several  concluding  comments  are  noteworthy  concerning  the  use  PARAM-1 
through  PARAM-6.   The  alpha  -  beta  parameter  plane  provides  more  informa- 
tion than  any  other  representation  for  system  studies  involving  the  simul- 
taneous variation  of  two  parameters.   Several  different  PARAM-1  plots 
are  considered  to  be  most  useful  for  initial  system  studies,  when  they 
are  accompanied  by  their  corresponding  PARAM-6  graphs.   In  this  manner, 
a  three  dimensional  space  of  alpha,  beta,  ?"d  omega  can  be  visualized  for 
each  of  the  constant  magnitudes  used.   The  additional  incorporation  of  one 
or  two  PARAM-2  graphs  will  describe  the  effect  of  a  variation  in  the  fourth 
dimension.   Stated  in  other  words,  the  information  in  PARAM-2  shows  how  the 
three  dimensional  PARAM-1,  PARAM-6  surface  varies  in  magnitude.   PARAM-3 
and  PARAM-5  curves  yield  a  more  sensitive  portrayal  of  the  magnitude  varia- 
tions; however  these  are  restricted  to  a  single  parameter  variation.   They 
are  useful  for  accurately  presenting  a  maximum  peak  resonance  magnitude 
when  this  information  is  not  available  from  the  peak  resonant  tangent 
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locus  of  PARAM-1.   PARAM-4  effectively  shows  a  series  of  Bode  curves 
as  a  function  of  one  system  parameter.   As  such  they  would  be  advantag- 
eous for  intermediate  studies,  when  the  final  adjustments  are  being 
contemplated.   Finally,  the  familiar  Bode  plot  representations  of  PARAM-7 
show  the  most  accurate  frequency  response,  and  are  therefore  most  useful 
in  the  terminal  stages  of  design. 
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APPENDIX  J 

C*******  PROGRAM  PARAMS  ********** 

C   THIS  PROGRAM  CONVERTS  A  GIVEN  TRANSFER  FUNCTION  TO  A  SET  OF  FREQUENCY 
CDOMAIN  REPRESENTATIONS  FOR  ANALYSIS  AND  DESIGN  STUDIES.  THIS  TRANSFER 
CFUNCTIUN  IS  REAO  AS  A  RATIO  OF  TWO  COMPLEX  FREQUENCY  POLYNOMIALS,  WITH 
CTWO  VARIABLE  PARAMETERS  ALPHA  AND  BFTA  EMBEDDED  IN  THE  COEFFICIENTS. 
CTHE  PROGRAM  THEN  TRANSFORMS  THIS   TRANSFER  FUNCTION  INTO  THE  ANGULAR 
CFREQUENCY  DOMAIN,  YIELDING  AN  EQUATION  IN  TERMS  OF  FOUR  VARIABLES: 
CALPHA.  BETA.  MAGNITUDE,  AND  0»«EGA.  SEVEN  TWn-DI  MENS  I QNAL  GRAPHICAL 
CREPRESENTATIONS  ARE  GENERATED  BY  FIXING  ONF  VARIABLE  AS  A  CONSTANT: 
CASSIGNING  DESIRED  CURVE  VALUES  TO  ANOTHER;  AND  LETTING  THE  REMAINING 
CTWO  VARIABLES  BE  THE  ORDINATE  AND  ABCISSA  VALUES.  A  THIRD  COEFFICIENT 
CVARIABLE  PARAMETER.  GAMMA,  MAY  ALSO  BE  EMBEDDED  INTO  THE  POLYNOMIAL. 
CIN  THIS  CASE,  THE  ABOVE  PROCESS  IS  REPEATED  FOR  EACH  VALUF  OF  THE  THIRD 

CVARIABLE.  DATA  CARD  INFORMATION  FOLLOWS   

C******************DATA  CARD  INFORMATION   SECTION********************* 
CCARD  FORMAT   SYN " " 


81  10 


NORN 
NORO 
MODORD 


NALPHA 

N3ETA 

NMAG 

NOMEGA 

NGAMMA 


LX 
LY 

CONALP 
CONBET 


USAGE  OF  SYMBOL 
NUMBER  OF  DESIRED  RUNS.  FACH  MUST  HAVE  A  COMPLETE 
SET  OF  DATA  CARDS  LISTED  BELOW.  NOT  EVERY  SYMBOL 
IS  USED  IN  ALL  CASES.  LEAVE  UNUSED  SYMBOLS  BLANK. 
00  NOT  OMIT  ANY  BLANK  CARDSj  UNLESS  SO  SPECIFIED. 
INSERT  ADDITIONAL  CAROS  BEHIND  EACH  FULL  CARD. 
NUMBER  OF  GRAPHICAL  PRESENTATIONS  TO  BE  PROGRAMMED 
NUMERICAL  TITLES  (1  TO  7)  OF  EACH  PARAM  TO  BE  PUN; 
MUST  AGREE  IN  NUMBER  WITH  KPARAM  APOVE. 
ORDER  OF  NUMERATOR  POLYNOMIAL  (LIMITED  TO  49 ) 
ORDER  OF  DENOMINATOR  POLYNOMIAL  (LIMITED  TO  49) 
Y-AXIS  VARIABLE,  AND  CURVE  VARIABLE  CONTROL 
MODORD=1:  ALPHA  IS  ORDINATE  VARIABLE  (PARAM  1-6) 
M0D0R0=2:  BETA   IS  ORDINATE  VARIABLE  (PARAM  1-61 
MCD0RD=1:  ALPHA  CURVES  WILL  BE  PLOTTED  (PARAM  7) 
M0DURQ=2:  BETA  CURVES  WILL  BE  PLOTTED  (PARAM  7) 
NUMBER  OF  ALPHA  CURVES  DESIRED  (LIMITED  TO  16) 
NUMBER  OF  BETA  CURVES  DESIRED  (LIMITED  TO  IS) 
NUMBER  OF  MAGNITUDE  CURVES  DESIRED  (LIMITED  TO  16) 
NUMBER  OF  OMEGA  CURVES  DESIRED  (LIMITED  TO  16) 
POSITIVE  INTEGER  DENOTES  NUMBER  OF  THIRD  PARAMETER 
VALUES  TO  BE  EMBEDDED  INTO  EACH  PARAM  RUN.  INSERT 
COEFFICIENTS  FOR  THESE  GAMMA  VALUES  ON  SEPARATE 
CARDS  DIRECTLY  BEHIND  REST  OF  COEFFICIENTS;  SEE«19 
NEGATIVE  INTEGER  DENOTES  NUMBER  OF  PARAM-7  VECTOR 
ADDITIONS  FOR  GENERALIZED  BLOCK  DIAGRAM  STUDIES. 
THE  (NEGATIVE)  VALUE  OF  NGAMMA  DETERMINES  HOW  MANY 
BODE  CALCULATIONS  WILL  BE  MADE  BEFORE  PLOTTING  THF 
VECTOR  RESULTANT  AS  ONE  MAGNITUDE,  AND  ONE  PHASF. 
INSERT  ADDITIONAL  SETS  OF  COEFFICIENTS  AT  FND  . 
MAGNITUDE  VALUE  FORM  CONTROL 

IMAG=0:  ALL  VALUES  IN  DECIBELS  (PAPAM  7  THIS  FORM) 
IMAG=1  MEANS  ALL  VALUES  IN  ACTUAL  MAGNITUDE 
IMAG=2  MEANS  ALL  VALUES  IN  MAGNITUDE  SQUARED 
PRINTOUT  CONTROL 
IWRITE--1  MEANS  NO  PRINTOUT 

IWRITE=0  MEANS  ONLY  INPUT  DATA  WILL  BE  PRINTED 
IWRITE=1  TO  900  YIELDS  INPUT  DATA,  GRAPHIC4L 
INFORMATION,  AND  COMPUTATIONAL  RESULTS.  THE  NUMBER 
OF  COMPUTATIONS  BETWEEN  EACH  PRINTOUT  IS  DICTATED 
BY  THE  VALUE   OF  IWRITE.  IF  IWRITE  IS  GREATER  THAN 
900,  ONLY  INPUT  DATA  6  GRAPH  INFO.  WILL  BF  PLOTTED 
FREQUENCY  SCALING  CONTROL 
IGRAPH=-2:  NO  GRAPHS  WILL  BE  PLOTTED. 
IGPAPH=-1:0NLY  PARAM  7  WILL  BE  PLOTTED  ON  LOG  GRID 
IGRAPH=0:  PARAMS  6E7  WILL  BE  PLOTTED  ON  LOG  GRID. 
IGRAPH=1 :PAPAMS  4,6, £7  WILL  BE  PLOTTFD  ON  LOG  GRIP 
IGPAPH=2:  PARAMS  4&7  WILL  BE  PLOTTED  ON  LOG  GRID. 
PARAM  7  PHASE  PLOT  CONTROL 
IPHASE=-1  MEANS  NO  PHASE  CURVES 

IPHASE=0:  PHASE  CURVES  SUPERIMPOSED  ON  MAGNITUDE 
IPHASE=l:  PHASE  CURVES  PLOTTED  ON  SEPARATE  GRAPH 
PARAM  7  PHASE  PLOT  SCALING  CONTROL 

IQUAD=0  YIELDS  PHASE  CURVES  FROM  0  TO  -180  DEGREES 
IQUAD=1  YIELDS  PHASE  CURVES  FROM-90  TO  *90  DEGREES 
PARAM  7  POLAR  PLOT  CONTROL 
IPOLAR=0  MEANS  NO  POLAR  PLOT 

IP0LAR=1  GIVES  POLAR  PLOT  ON  SEPARATE  GRAPH 
LENGTH  OF  X-AXIS   (LIMITED  TO  9  INCHES) 
LENGTH  OF  Y-AXIS   (LIMITED  TO  15  INCHES) 
VALUE  OF  ALPHA  WHEN  USED  AS  A  CONSTANT 
VALUE  OF  BETA  WHEN  USED  AS  A  CONSTANT 


C  CUNMAG   MAGNITUDE  VALUE  WHFN  USED  AS  4  CONSTANT:  SEE  IMAG 

C  CONOMt   VAl'JE  OF  OMEGA  WHEN  USFD  AS  A  CONSTANT  (NOT  ZERO) 

C***************r>ATA  CARD  INFORMATION  SECTION  t CONT I N'IFO I  ******  ********* 
CCAKD  FORMAT   SYMROL         USAGE  OF  SYMBOL 

C   6   3E10.3   ALPHMN   MIN.  VALUE  OF  ALPHA  WHEN  USED  AS  AN  AXIS  VARIABLE 
Alphmx   MAX.  VALUE  OF  ALPHA  WHEN  USEn  AS  AN  AXIS  VARIABLE 
BETAMN   MIN.  VALUE  OF  BETA  WHEN  USED  AS  AN  AXIS  VAPIARIF 
RETAMX   MAX.  VALUE  OF  BETA  WHFN  USED  AS  AN  AXIS  VARIABLE 
XMAGMN   min.  MAG.  VALUE  WHEN  USED  AS  AXIS  VAR.;  SFF  I*AG 
XMAGMX   MAX.  MAG.  VALUE  WHFN  USED  AS  AXIS  VAR.;  SEE  IMAG 
OMEC.mn   MIN.  OM£GA  VALUE  WHEN  USED  AS  AXIS  VARIABLE 
CMEGMX   max.  OMEGA  VALUE  WHEN  USED  AS  AXIS  VARIABLE 
C  NOTE: IF  LUC  GRID  IS  REQUESTED, < SEE  IGRAPH),  OMEGmn 

WILL  ^£  AUTO.  CHANGED  TO  NEARFST  POWFP  OF  10,  AND 
^EGMX  CHANGED  LIKEWISE.  THF  NUMBFf  OF  DECADES  IS 
AUTOMATICALLY  COMPUTED  FROM  THESE  NEW  VALUES. 
C   7   3E10.3   ZALPHA   DESIRFD  VALUES  FOR  ALPHA  CURVES  (SEE  NALPHAI 
C   8   Pfc.10.3   ZBETA    DfcSIRFO  VALUES  FOR  BETA  CURVES  (SFF  NRFTA) 
C   9   8E10.3   ZMAG     DFSIREO  VALUES  FOR  «»AGNITUDC  CURVES(SEE  NMAGf.IMAC) 
C  10   6E10.1   ZJMFGA   DESIRED  POSITIVE  OMEGA  CURVE  VALUES  (SEF  MOMFGA ) 
C  11   SE10.3   HA       NllMfHATOP  CONSTANT  COEFFICIENTS  IN  ASCENDING  ORDER 
C  12   8E10.3   HAT      NUMERATOR  ALPHA  COEFFICIENTS  IN  ASCENDING  ORDER 
C  13   6E10.3   HBT      NUMERATOR  BETA  COEFFICIENTS  IN  ASCENDING  ORDER 
C  U   8E10.3   HCT      NUM.  ALPHA-BETA  COEFFICIFNTS  IN  ASCENDING  ORDER 
C  15   8E10.3   HO       DENOMINATOR  CONST.  COEFFICIENTS  IN  ASCFNDING  ORDER 
C  IS   HE10.3   HOT      DENOMINATOR  ALPHA   COEFFICIENTS  IN  ASCENDING  OROFP 
C  17   ^EIO.3   HET      DENOMINATOR  BETA  COEFFICIFNTS  IN  ASCENDING  ORDER 
C  18   3E10.3   HFT      DtNO".  ALPHA-BETA  COEFFICIENTS  IN  ASCENDING  ORDFR 
C«*********  QMJT  CARDS  «  19  THRU  *  27  IF  NGAMMA  LESS  THAN  ONE. 
C  19   8E10.3   ZGAMMA   VALUES  OF  THE  THIRD  PARAMFTFP;  MUST  MATCH  NGAMMA 
C  20   8E10.3   GAMHA    Nil*.  Gamma-CONST.   COEFFICIENTS  IN  ASCENDING  OROFP 
-10.3   GAMhAT   NUM.  uamma-ALPHA    COEFFICIENTS  IN  ASCENDING  ORDEP 
C  22   SE10.3   C.AMmot   NUM.  GAMMA-BETA     COEFFICIENTS  IN  ASCFNDING  ORDFR 
C  25       6E10.3   GAMHCT   NUM.  G AMMA-ALPH A-BE T A  COEFFS.     IN  ASCENDING  0RDEc 
C  Z<*       8E10.3   GAMHD    DfcN.  G  4MMA-C  ON*  T .  COEFFICIENTS  IN  ASCENDING  ORDEP 
C  25   3E10.3   GAMHDT   DEN.  GAMma-ALPHA   COEFFICIENTS  IM  ASCENDING  ORDER 
C  26   SE10.3   GAMHET   Ot * .  GAMma-BETA    COEFFICIENTS  IN  ASCFNDING  ORDFR 
C  27   ?E10.3   GAMHFT   oen.  GAMMA-ALPHA-BETA  COEFFS.   IN  ASCENDING  ORO^R 
C  26   t,A«       ITITLE   FIRST  LINE  GRAPH  TITLE  (COLUMNS  1THPU43) 
C  29   6AS       ITITLE   StCONO  LINE  OF  GRAPH  TITLE  (COLUMNS  1THRII4S) 

THE  NUMBER  OF  TITLES  FOP  °AKAM  7  M|JST  MATCH  THF 
NIimbfr  OF  GRAPHS  REQUESTED.  (SFF  IPHASF  f,  IPOLAP) 
C  30   6A8       JTITLE   1ST  1/2  PHASE  TITLE;OMIT  UNLFSS  NP  AR  4«=  7r.  I  P'(  A* F  =  1 
C  31  e>Uf>  JTITLE   2NC  1/2  PHASE  TITtE:0*lT  UNLESS  NPAR  Am=7£  I PHAS  F=  1 

C  32   6A8      KTITLl   1ST  1/2  POLAR  TITLF;OMIT  UNLFSS  NP  AR  AM  =  7F.  1  POL  AR=  1 
C  '3   >-A3      -TITLE   2ND  1/2  POLAR  TITLE:OM[T  UNLFSS  NP  AR  AM  =  7r.  1  POL  AR=  1 
r ************************************************  *********************** 
'***«**»**«**•*******«******************************«*******«********««« 
FEAL*BITITLF(12),JTITLE(12),KTITLF(12),J0KFR(9)/1H  , SHPAR AM- , 7HF|jp 
ETHER,5HALPHAf4HBETA,5H0MFGA,BHMAGNITUD,6HG4MMA,BHSFPARATF/ 

REAL  DM/1H1/,LA(  661/4H  0.1, 4H  C.2.4H  0.3.4H  0.4.4H  0.C,4H  0.6.4H  C 
f..7,4H  0.8.4H  0.9.4H  l.C»H  1.0, 4H  2  .  0 , 4H  3.0.4H  4.0,4H  5.0,4H  6.0, 
f,<.H  7.0.4H  8.0.4H  C>.0,4H1C.O,4H10.0,4H20.0  ,4H30 . 0 , 4H40 . 0, 4H5Q .0 , 
E  4H60.  0.4H71.0,  :.MfiC.0,4H90.0,4H100.  , 4H100.,4H2  00.,4H300.,4H400., 
t4H500. ,4H6  00. .4H700.  ,4H800. ,4H900. ,  4H1 000 ,4H 1000 , 4H2 000 , 4H 3000, 
e4H4000,4H5  00G,4H60C3,4H700  0,4HfiOOO,4H9  000,lH  , 1H1 , 1H2.1H3, IH4, 1H5, 
C1H6,1H-',H8,1H9,2H10,7H11,2H12,2H13,2H14,?H1S,?H16/ 

DIMENSION  JPAKAM(f.  )  ,  A(  50)  ,B<  50)  ,0(50  )  ,F(  SO)  ,HA(  SO)  ,HAT(50I  ,HBT<  50) 
F.  ,HCT(  50)  ,HD<50)  ,H0T(  50)  ,HET(50)  ,HFT<  50)  ,OLDHA(50(  ,OL0HAT(50)  ,01.  DHR 
&T( 50)  ,OLDHCT( SO)  ,CLDHD( 50) ,OLDHDT(50) , OLDHETISO ) ,01 DHFT( 50),GAmha( 
C50)  ,GAMHAT(  50)  ,GAmh-)T(  50  )  ,GAMHCT(50)  ,GAMHD(50),GAMHDT(5  0),GAMHFT(S 
&0)  ,GAMHFT( 50)  ,IC( 16)  ,ZALPHA(16),ZBETA(16)  , ? 0 AMM A  I  1 M  , 7 MAGSQ ( I  6 ) , 70 
{.M£GA(16),S(16,90C),T(  16,900  1  ,  U  (  1  6  ,  900  I  ,  V  (  1*  ,  900  )  ,  W  (  1  6  ,  900  )  ,  X  (  900  )  , 
IVI900) 
CALL  CANCcL(2) 
EQUIVALENCE! S,T,U,  V,W) 

1  FORMAT(PIIO) 

2  FORMAT! 1CX.8I 10) 

3  F0RMAT(6A8) 

4  FORMAT( 10X.SA8) 

5  FORMAT!   PE10.3) 

fc  FORMAT! 10X.1P10E12.3) 
C *************************  I NPUT  DATA  READ  IN  *************************** 
READ! 5, 1 ) MRUN 
DO  1000  MUN  =  1,NRUN 

READI 5,1 )  KPARAM,(JPARAM(I),I=1,KPARAM) 

REAO! 5, 1 )NORN,NORO,MODORD,NALPHA,NBFTA,NMAG,NOMCGA,NGAMMA 
REAO!5,l)  IMAG  .IWPITE  , I GR A PH , I PH ASE , I QUAn, I  POL AR , LX , L Y 
READ! 5,5 (COMAL P . CORBET ,CONMAG,CONOMF 


RE AD ( 5, 5 ) AL PHMN, ALPHMX, BETAMN.BE TAMX ,XMAGMN,X MAGMX , OMEGMN, OMEGMX 
READ<5,5)  (ZALPHA( I ) ,I=1,NALPHA) 
REA0(5,5)  (ZBETAU  ).I  =  1,NBETA) 
REAO(5,5l  (ZGAMMA( I  ),  1  =  1 .NGAMMA I 
READ(5,5)  (ZMAGSQ< I) ,I=1,NMAG) 
REAO(5,5)  (ZOMEGAl I ), I =1 , NOMEGA ) 
LOE=NORN*l 
MOE=NORD«-l 

READ(5,5>  (HA  (I),  1  =  1, LOE) 
READ(5,5)  (HAT(I),  1=1, LOE) 
REA0I5.5)  (HBT(I),  1  =  1, LOE) 
REA0(5,5)  (HCT(I),  1=1, LOE) 
REA0(5,5)  (HO  ( I ),  1=1, MOE) 
READ(5,5)  (HOT(I),  1  =  1, MOE) 
READ(5,5)  (HET(I),  1  =  1, MOE) 
READ(5,5)  (HFT(I),  1=1. MOE) 

IF( (NGAMMA). GE.(  1) )  READ(5,5)  ( ZGAMMA( I ),  1= 1 , NGAMMA ) 
IF((NGAMMA).GE.(  1) )  READ(5,5)  (GAMHA  (I), 1=1, LOE) 
IFt (NGAMMA). GE. ( 1) )  REAO(5,5)  ( GAMHAT (  I  )  ,  1  =  1 ,LOE ) 
1F( (NGAMMA) .GE .( 1) )  REAO(5,5)  ( G AMHBT( I ) , 1= I , LOE ) 
IF( (NGAMMA). GE .( 1) )  READ(5,5)  ( GAMHC T( I ) , 1= 1 ,L OE ) 
IF((NGAMMA).GE.(  1) )  READ(5,5)  (GA"HD  (I), 1=1, MOE) 
1F( (NGAMMA). GE .  (  1)  )  READ<5,5)  (GAMHOTI I ) , 1  =  1, MOE) 
IF( (NGAMMA). GE.( 1) )  READ(5,5)  ( GAMHE T( I ) , 1= 1 , MOE ) 
IF( (NGAMMA). GE.( 1) )  REA0(5,5)  ( GAMHF T( I ) , I = 1 , MOE ) 
IF( (IWRITE).LT.(O)  )  GO  TO  2* 
i************************INPUT  DATA  PRINT  OUT  SECTION***************** 
WRITE (6. 7)  IRUN,NORN,NORD,KPARAM,( JOKER ( 2 ), JP ARAM ( l),I=l,KPARAM) 

7  FORMAT( 1H1.10X,' *****  THE  INPUT  DATA  FOP  RUN',!?,'  IS  *****•,//, 
C10X, 'ORDER  OF  NUMERATOR :'. I  2 ,//, 10X, • ORDER  OF  DENOMINATOR:',^,//, 
tlOX.'THE  FOLLOWING', 12, •  GRAPHICAL  CONFIGURATIONS  WILL  RF  COMPUTED 
6:  •  ,//,20X,8( A8,I 1.3X) ,/) 

IF( (NGAMMA). GE . ( 1) )  WRITE(6,8)  NGAMMA, ( ZGAMMA ( I ), I =1 , NGAMMA ) 

8  FORMATdOX,  12,'  GAMMA,  OR  THIRD  VAR.  PARAMETER  COEFFICIENTS  WILL  B 
CE  USED.' ,//,10X, 'THE  GAMMA  VALUES  ARE: • , I P8E1 2 . 3 , / , 10X, 1 PBE 1 2 . 3, / ) 

WRITEI6.10) 
10  FORMAT! 10X, 'THE  COEFFICIENTS  FOR  THE  NUMERATOR  ARE:',/) 

DO  11  1=1, LOE 

J=l-1 

IFINGAMMA  .LE.O)  WRITE(6,12)  HA (  1  )  ,HAT ( I )  .HRT <  I ) , HCT ( I ) , J 
12  F0RMAT(20X,'( • , 1 PE 10 . 3 , • ♦ • , 1 PE 10. 3 , '*ALPHA  ♦  •  ,  1PF  10  .  3, • *BET< 
)  *  S« 


E, 1PE10.3,'*ALPHA*BETA  )  *  S**( ',  12 ,')',/ ) 

IFINGAMMA  .GT.O)  WRITEI6,  9)  HA (  I ) , GAMHA (  I ) , HAT  (  I )  ,  GAMHAT  (  I )  ,  HBT  (  I 
O.GAMHBTI  I  ).HCT(  I  )  ,  GAMHC  T(  I  )  ,J 
9  FORMA T( 15X,M(',1PE10.3,'+',1PE10.3, «*GAMMA)*( • , 1PE10.3, •♦• , IPE10. 
C3,  «*GAMMA)*ALPHA  ....',// ,20X. '«•(',  I  PE  10.  3  ,  •  ♦  '  ,  1  PE10  .  3  ,  •  *  GAMMA  )  *BE 
£TA«-(  ',1PE10.3,'*',1PE10.3,'*GAMMA)*ALPHA*BETA)**(  •  ,12,  •  )  ',/) 
11  CONTINUE 

WRITE(6,13) 

13  FORMAT( /10X, • THE  COEFFICIENTS  FOR  THE  DENOMINATOR  ARE:',/) 
DO  14  1  =  1, MOE 

J=I-1 

IF(NGAMMA  .LE.O)  WRITE(6,12)  HD( I ) ,HDT ( I ) ,HET ( I ) , HFT ( I ) , J 
IF(NGAMMA  .GT.O)  WRITEI6,  9)  HD( I ) ,GAMHD( I ) , HOT ( I ) ,G AMHOT (  I  I , HET ( T 
E), GAMHE T( I),HFT(  I  )  ,GAMHFT( I ) ,J 

14  CONTINUE 
WRITE(6,15) 

15  F0RMATU5X,'   NORN       NORD      MODORD     NALPHA      NBFTA     N^A 
EG        NOMEGA   NGAMMA  ',/) 

WRITE(6,2  )  NORN, NORD,MODORD,NALPHA,NBETA,NMAG, NOMEGA, NGAMM<\ 
WRITE(6, 16) 

16  FORMAT! /,17X, • IMAG        IWRITE     IGRAPH      IPHASE       IQHAD      IPT 
ELAR      LXAXIS     LYAXIS   ',/) 

WRITE(6,2)   IMAG  ,  I  WR  I  TE  ,  IGR  APH  ,  I  PHASE  ,  I  QUAD,  I  POL  AR  ,  |_X  ,L  Y 
WPITE(6,17) 

17  F0RMAT(/,12X, •  CCNALPHA     CONBETA       CONMAG      CONOMEGA   ',/) 
WRITE (6, 6  )CONALP,CONBET,CONMAG,CONOME 

WRITE(6,18) 

18  F0RMAT(/,12X,'  ALPHMIN     ALPHMAX     BETAMIN       BETAMAX      XMAG 
EMIN     XMAGMAX      OMEGAMIN      OMEGAMAX   ',/) 

WRITE(6,6) AL PHMN, ALPHMX, BE TAMN, BET  AM x,XMAGMN,XMAGMX,OMEGMN,OMFGMX 
WRITE (6, 19)  JOKER (A) , ( ZAL PHA ( I ) , I =1 , NAL PHA ) 
WRITE(6,19)  J0KER(5) ,(ZBETA  ( I )  ,  I  =  1 , NBETA  ) 
WRITE  (6,  19)  JOKER  (7)  ,(  ZMAGSOU  I  ,I  =  1,NMAG) 
WRITE (6,  19)  J0KER(6) ,(ZOMEGA( I  ), 1  =  1 , NOMEGA) 

19  F0RMAT(/,10X, AH, 'CURVE  VALUE S: ' , 1P8E 12 . 3 , /3 IX , 1 P8 E 12 . 3 ) 
IF(IWRITE)  20,20,22 

20  WRITE(6,21) 

21  F0RMAT( /,10X, 'THERE  WILL  BE  NO  PRINTOUT  OF  COMPUTED  RESULTS'! 


GO  TO  24 
11    WRITEI6, 23) IWRI TE 

li    FU«M4T(/,10X, 'CJmpuTcT  RESULTS  WILL  HE  PRINTED  EVERY', 13,"  STFPS'l 
C*»****************iNITIALmNG  COMPUTATIONAL  SETT  I  NHS  ****************** 

24  IF(MOOORO.GP.2)  or  TO  25 
OUM«Yl=«ETAMN 
OU*MY2=BETA*X 
bETAMN=ALr>HMN 
BETAMX=ALPHMX 
ALPHM^=nUM«Yl 
ALPH«X=OU«"Y2 

25  DO     1000     IPARAM=1  ,  KPARAM 
NPARA«=JPAPAM<IPARAM) 
Mr,AMMa  =  NGAMMA+l 

00  1000     IGAMMA=1 ,MGA«MA 
GAMMA=ZGAMM4( [GiM^A) 
IKNGAMma.ge.  1»     GC    TO    517 

26  XS=  (  ALPF-MX-ALHHMN)  /LX 

IF(  (ABS( XS)  )  .LE.t  1.0E-10I )  XS=1.0E-10 

YS=(BETAMX-«fcTAMN) /LY 

IF(  (A8SI YS)  ) .Lf  .  (  l.OE-10)  )  YS=1.0F-10 

I X=-HETAMN/YS 

I Y=-ALPH«N/XS 

M00=1 

JCURVE=C 

IPAR7=1 

RCAOI 5,3)1 ITITLE (11,1=1,12) 

IF(  NPARAM.ME.7)G0  TO  34 

IF(  I»AG-1)27,2H,2<- 

27  0Amax=xmag"X 
CBMIN=XMAG*N 
GO  T~ 

29  CRMAX=20*ALflG10(  XMcGMX) 
GrtMIN=20*ALPG10( XMAG-N) 
GO  TO  30 

Z-i  CBMAX=10*ALDG10( XMAG«X) 
DBHIN=1C*AL0G1C( XMARMNI 

30  XS=(OM*AX-')BMIN)  /LX 
IY=0B"AX/XS 

OV     41      1=1.16 

S<  I ,  1)=0.0 
T(  I ,  1)=0.0 
U(  I  ,11=0.0 
V(  I ,  11=0.0 

31  W( I ,1)=0.0 

IF( IPHASf.EO.OI     GC     TO     32 
REAiM  5, 3)(  JTITLFI  1  ),  1  =  1,12) 
i2     IF( ( IPOl AS ) .LE.( 0) )    GJ    TO    33 
READ! 5,3) (KTITLFI  I  ),I  =  1,12) 

33  IF(MO0ORO-2)45,44,44 

34  GO    TO    (  35,40,  '.I,  <•*,  47,51,  59,59)  .NPAKAM 

35  NLI*IT=NO*tGA 

36  IF(  I-AG-U37,  3rt,  39 

37  OUMMY1=10.0**(CONMAG/10.0) 
IF(NPARAM-1  159,59,43 

3d  DUMMYl=CONHAG*CONHAG 

IF(  NPARAM-l)  59,59,43 
3=5  0UMMY1=CCN«AG 

IF(  NPARAM-D59,  59,^3 
40    NLIMIT=N«AG 
Hi     OMEGA=CrNOME 

IF(NPAkAM-2)59,5«5,a2 
^2  XS=( xmagmx-xmagmnS /LX 

|Y=-X«*AG«N/XS 
-3  IF(mgotbd-2)44,<,5,45 

44  NL1HIT=NBETA 
1F(NPARAM-6159,53,55 

45  NLI^IT=NALPHA 
IF(NPARAM-6) 59, 53,55 

46  NLI"IT=NM4G 
XS=(OMECMX-OMFGMM  /LX 
IY  =  -0»«EG"N/XS 

W     IF(MQ00R0-2)4R,49,49 

48  C0NST=CnN8£T 
IF(NPARAM-<,)52,  52,50 

49  CONST=CTNALP 

IF( NPAPAM-4152, 52, 50 

50  XS=(XMAGMX-XMAGMN) /LX 

1  Y  =  -X«AG*N/XS 


NLIMlT=NOMEGA 
GO    TO    59 

51  XS=!OMEC-MX-OMEGMN)/LX 
IY=-OMEGMN/XS 

GO  TO  36 

52  IF! IGRAPHI59, 59.54 

53  IF(  (  IGftAPH.LT.  0)  .OR.  (  I  GR  APH.GT  .  1 )  )  GO  TO  59 

54  IY=LX-IX 

55  IF(  IGRAPH.LE.-2I  GO  TO  59 

C**************«****LOG  GRID  GRAPHICAL  CONSTRUCTION********************* 
IFIOMEGMN.LT.IOOO. )YL0GMN=100. 
IF!0MEGMN.LT.100.  )YL0GMN=10. 
IFI0MEGMN.LT.10.   )YL0GMN=1. 
IF(OMEGMN.LT.l.     )YL0GMN=0.1 
IF!0MEGMX.GT.1.0  )YL0GMX=10.0 
1F( OMEGMX.GT. 10.0 )YL0GMX=1 00.0 
IF (OMEGMX.GT. 100. >YLOGMX= 1000.0 
IFIOMEGMX.GT.IOOG. ) YLOGMX= 10000. 0 
0UMMY3=  YLOGMX/YLCGMN 
NUMOEC=ALOG10(OUMMY3) 
IX=0 

DUMMY4=NUMDEC 
0UMMY5=LY 
YS=DUMMY4/DUMMY5 
IFINAL=LX+1 
JFINAL=LY*1 
DO  57   IA=1,IFINAL 
00  56   JA=1,JFINAL 
X! JA)=XS*( IA-l-IY) 

56  Y( JA)=YS*( JA-1) 

CALL  DRAW( J  A , X, Y , MOD ,0 ,L A < 50 ) ,ITITLE,XS,YS  ,0,IY,2,2,LX,LY,0,LAST) 

57  M0D  =  2 

DO  58  IB=1,NUMDEC 
DO  58  JB=1,10 

J=ALOG10(  YLOGMN)*1C.«-10.*<  IB-1)*10.+JB 
YLOG=JB 

DO  1058  KB=1,IFINAL 

X!KB)=XS*!KB-1-I Y) 

1058  Y<KB)=AL0G10!YL0G)+IB-1 

58  CALL  DRAWIKB,X,Y,MOD,0,LA( J) , I T I TLE , XS , YS  , 0 , I Y , 2 , 2 , LX , L Y , 0, L AST ) 
IF(  (  IPAR7) .NE.( 1 ) )  GO  TO  530 

59  IF!  IWRITE.LE.O)  GO  TO  200 
C»*****************GRAPHICAL   INFORMATION  PRINTOUT  S ECT I TN************** 

WRITE<6,60)  IRUN,  NPARAM 

60  FORMAT! 1H1,10X,' **********THE  PRELIMINARY  DATA  FOR  RUN  NUMBER'  1 12, 
£•,  PARAMETER  NUMBER', 12,'  IS    **********•.///) 

IF(NGAMMA.GE.l)  WRITE  (6,61)  I  GAMMA (  I  GAMMA) 

61  FORMAT! 10X, 'THE  FOLLOWING  CALCULATIONS  WILL  BE  3ASE0  ON  THF  THIRn 
^VARIABLE  PARAMETER  VALUE  ,  GAMMA  =  « , 1 PE 12  .  5 , // I 

WRITE  (6,62)  ( ITITLE( I ) ,1=1,12) 

62  FORMAT( 10X, 'THE  GRAPH  WILL  BE  T I TLED  :  •  , // , 10X ,6 48 , / / , 10X , 6  Aft , / ) 
IFINPARAM.NE.7)  GO  TO  80 

WRITE(6,83) 

WRITE(6,   63)nBMIN,DBMAX,XS 

63  FORMAT! 10X, 'MAGNITUDE,  THE  DEPENDENT  VARIABLE,  WILl  BE  PLOTTED  ON 
STHE  ORDINATE  FROM •  ,F 10 .4 , •  T0',F10.4,«  AT',F10.4,'  D8/D I  V  I S ION ' , / ) 

IF(IPHASE)  73,64,66 

64  WRITE(6,65) 

65  FORMAT* 10X, 'PHASE  CURVES  WILL  BE  PLOTTED  ON  MAGNITUDE  GRAPH',/) 
GO  TO  68 

66  WRITE(6,67) 

67  FORMAT! 10X, 'PHASE  CURVES  WILL  BE  PLOTTED  ON  SEPARATE  GRAPH',/) 

68  I F «  IQUAD.GT.O)  GO  TO  70 
WRITE(6,69) 

69  FORMAT! 10X, 'FROM  ZERO  DEGREES  THROUGH  -90  TO  -130  DEGRFFS',/) 
GO  TO  72 

70  WRITE(6,71> 

71  FORMAT! 10X, 'FROM  -90  DEGREES  THROUGH  -180  TO  -270  DEGRFFS',/) 

72  WRITE(6,62>  (  JT I TLE( I )  ,1 =1 , 1 2 ) 

73  IF( IPOLAR.LE.O)  GO  TO  75 
WRITE(6,74)  XS 

74  FORMATdOX,  'A  POLAR  PLOT  OF  MAGNITUDE  VS.  PHASE  WILL  BE  PtOTTFD,  W 
&ITH  MAGNITUDE  SCALED  4T',F10.5,'  DECIBELS   PER  HALF  INCH',/) 

WRITE(6,62)  (KTITLE!  I)  ,1=1,12) 

75  IFIMODORD-2)  127,  121,121 

76  FORMATdOX, A8, •  ,  THE    DEPENDENT  VARIABLE,  WILL  BE  PLOTTED  ON  THF 
(•.ORDINATE  FROM',  1PE10.  3,  •  TO' ,  1  PE  10  .  3  ,  •  WITH  A  SCALF  OF  •  ,  1PF  1  0.  *,  /  ) 

77  FORMAT! 10X.A8, • ,  THF  INDEPENDENT  VARIABLE,  WILL  BE  PLOTTFD  ON  THF 
&ABCISSA,  FROM' , 1PE10.3, •  TO' , 1 PE10 . 3 , •  WITK  - 


"H  A  SCALf  OF' , 1PE10.3,/ I 


76  FORMAT! 10X.AR, '  WILL  BE  HELD  CONSTANT  AT  THE  VALUF  OF •  ,  1  PE  1  2 . 5 , / ) 
79  FORMAT! 10X.A8, '  WILL  BE  GRAPHED  AS  ',13,'  CURVES,  CONSISTING  DF  TH 

&E  FOLLOWING  VALUES: '  ,16(  'CURVE', 13, •  =  • , 1 PF1 2 . 5  ,  /  )  ) 
8C  IFJ  IrtAG.E3.0I  WRITEI6.83) 

IF  (  IMAG.E3.  1  )  WR  ITK 6,8*) 

I^l  IMAG.E0.2)  WR  1TE( 6,85) 

81  IF(M001PD-2I  82,   9*, 1000 

82  WRITE(6,7o)  J.lKcRI  *)  ,  13ET  AMN.  BETAMX,  YS 

83  FORMAT! 1GX, 'ALL  MAGNITUDE  COMPUTATIONS  WILL  RE  IN  DECIBELS',/) 
8*  FORMAT! 10X, 'ALL  "AGNITUDE  COMPUTATIONS  ARE  ACTUAL  M AGN I TUDF • , / ) 
85  FORMAT! 10X, 'ALL  MAGNITUDE  COMPUTATIONS  ARE  MAGNITUDE  SQUARED',/) 

GO  TO  ! 86,86,121 ,130,130,121 ,200,200) ,NPARAM 
06  WRITE<6,77)  JOKER! 5) ,ALPHMN,ALPHMX,XS 

IF(NPARAM-l)   88,88,99 
88  wRITF(6,78)  JOK c K ( 7 ) , CONMAG 
9C  w^ITE(6,7?)J0KER!6),N0MFGA,( I ,ZOMFGA(I ) , I =1 , NOMEG A ) 

IFINPARAM-5120C, 125,200 
9*  WRITEI6.76)  JO* E R ( 5  I  ,  BET Amn , BE T AMX ,  YS 

96  GO  T0I97, 97,  127,  111,  111,127,200,200)  .NPARAM 

97  WRITE16.77)  JOK  E  R  (  * )  ,  ALPHMN,  AL  PHMX  ,  X  S 
IF (NPARAM-1J8H, 88,99 

99  WRITE(6,78)  JOK E R ( 6  )  .CGNOME 
101  WC  ITE!6,79)J0KER!7),NMAG,(  I .ZMAGSO!  I  )  ,I  =  l,NMAG) 

lF(NPARAM-*)200, 139,200 
111  WRITEI6.78)  J0KER!4) ,CCNALP 
113  WRITE!6,ll*) 
11*  FORMAT!  10X, • THE  COMBINED  COEFFICIENTS  OF  THE  NUMERATOR  APF',/) 

115  DO  116  I=1,L0E 

A!  I  )  =  HA(  I)  +HAT!  I  )*CC1NST 
8!  I  )=h3T(  I  )*HCT(  I  )*CONST 
KA=I-1 

116  WRITE(6,117)  AII),P(I),KA 

117  FORMAT!  ?0X,1H(  ,E  10.*,  3H  «•  ,  E10.*,16H  *  BETA  I  *  S**  (  ,  I  ?  ,  1  H  )  /  ) 
wRITE(6.136) 

OJ  119  1=1, MOE 

0! I )  =  HD< I)*HDT( I )»CONST 

fc(  I  )  =  HET|  I  )  «-HFT(  I  )*CONST 

KA=I-1 

119  wRITE!6, 117)  D(  I  ) ,E(  I  )  ,KA 

120  GJ  TO  (2CJ, 200, 2C0,U1,  90  ,  200  ,  I  39  ,200  )  ,  Nt>AR  AM 

121  wRITE(6,79)JOKER(5),NBETA,(I,ZBETA!I),t=l,NBfcTA) 
IF(NPARAM-b)  12*,  13P, 11  1 

12*  WRITE(6,78)  JOKEP ! 7) ,CONC*E 

125  WRITEI6.77)  JOKER! 7) , X"AGMN , XMAGMX , XS 

GO  TO  200 
12  7  WRI TEI6, 79) JOKE R ( * ) , NALPHA , < I , 7 ALPH A ( I  ) , I =1 , N AL °H A  ) 

IF(NPARAM-o) 12*, 138, 130 

130  WRITEI6.78)  JOK ER < 5 ) ,C CNBE T 

131  FORMAT!  lOX, 'BFTA   WILL  BE  HELD  CONSTANT  AT  THF   VALUE  OFt.FlO.*,/) 

132  WPITEI6.il*) 
DO  13*  1  =  1, LOE 

A!  I  )  =  HA(  I  I ♦HbT!  I  )*CCNST 
(i(  I  )  =  HAT(  I  )*HCT(  I  )*CONST 
KA=I-1 
13*  wPITc(6,135)A( I ) ,B( I ),KA 

135  FDRMAT(2CX, 1H(  ,E10.*,3H  ♦  ,E10.*,17H  *  ALPHA  )  *  S** (  ,  I  2 , 1 H ) / ) 
WRITEI6.136) 

136  FORMAT! 10X, • THF  CTKBINED  COEFFICIENTS  OF  THE  DENOMINATOR  ARF-',/) 
DO  137  I=1,MPE 

0!  I  )  =  H0(  I  H-HET!  I  )*CONST 
t!  I  )  =  HQT(  I  )+HFT(  I  )*CONST 
KA=I-1 

137  wRITE!6,135)  D( I ),E! I ) ,KA 
GJ  TO  120 

138  WRITEI6.78)  JUK ER ( 7 ) .CONMAG 

139  wRITE(6,l*0) 

1*0    FORMAT!  lOx,  "i«E',4,  THE     INDEPENDENT    VARIABLE, WHL     BE    PLOTTED    ON    THF 
F.ABCISSA'  ,/) 
IF (NP ARAM- 6) 14  1 , 1A2, 1** 

1*1     IF!  IGPAPHI  1*2,1*2, 1** 

1*2    WRITE(6,1*6)     OM£GMN,0MFGMx,XS 
GO    TO    200 

lO     IF!  (  IGRAPH.LT. 0)  .CR.  (  IGRAPH.GT.l)  )     GO    TO    1*2 

1**  WR ITE (6, 1*5)NU"DEC,YL0GMN 

1^5  FORMAT! 10X, 'WITH  A  SPAN  OF', 12,'  DEC AOFS , ST AR T I NG  AT  •, Fl 0 .*,//, 3CX 
I     • »**««M3TE***«*  '  ,// ,10X, 'FOR  THIS  LOG  PLOT,  IT  WILL  >M    NECESSARY 
CTO  SHIFT  THF  ORIGIN  TO  THE  LOWEP  RIGHT  HAND  CORNER ,',//, 10 X ,' W I TH 
CTHF  ABCISSA  RUNNING  UP  THE  RIGHT  HAND  FOGE,  AND  THF  ORDINATE  INCRF 


t,mr  ail.  liia  nur^^iiNo  ur  j  nc  niuni  numu     ruuc, 
CASING  TCWARDS  THE  LEFT  ALONG  THE  BOTTOM'//) 


K6  FORMAT! 10X, 'FROM  THE  INITIAL  VALUE  OF  •  ,  1PE10.3 ,  •  TO  THE  FINAL  VALUF 
G  0F'.1PE10.3,»  WITH  A  SCALE  OF  '.1PE10.3,'  RADIANS  PER  INCH*,/) 
C********************main  COMPUTATIONAL  SECTION  DO  LOOPS**************** 

200  IF! < IWRITE.LT.l) .OR.! IWRITE.GT.900) )  GO  TO  203 
WRITE(6,  201) IRUN,NPARAM,IWRITE 

201  FORMAT! 1H1.10X, '*********THE  OUTPUT  DATA  FOR  RUN  NUMBER1 , 12. • ,  PAR 
6AMETER  NUMBER', 12,'  I S**********  '  , //// ,  10X ,  I  2 ,  •  COMPUTATIONS  WILL 
GBE  PERFORMEO  BETWEEN  EACH  PRINTOUT',///) 

IF(NGAMMA.GT.O)  WRITE!6,202)  ZG AMMA( IGAMM A ) 

202  FORMAT! 10X, • THE  VALUE  OF  GAMMA,  THE  THIRD  VARIABLE  I S : • , 1 PEIO.  3, /  ) 

203  00  520  L0M=1,NLIMIT 
JCURVE=0 

DO  519  IC0EF=1,3,2 
ZSIGN=IC0EF-2 

IF! < IwRITE.LT.l) .DR.! IWRITE.GT.900) )  GO  TO  208 
WRITE<6,205  )LOM 

205  FORMAT!/, 10X, 'COMPUTATIONAL  RESULTS  FOR  CURVE  NUMBER', 13) 
IF  (NPARAM-7)2C6, 208,208 

206  WRITE(6,207  IZSIGN 

207  FORMAT! 1H+,55X, '  USING  A  QUADRATIC  FORMULA  COEFFICIENT  0F',F4.1,/) 

208  IP  =  0 
YMAGPK=0.0 
UMEGPK=0.0 
DO  51b  INCRMT=1,900 
GO  TO  ( 209,216,253,216,212,254,257,1000) ,NPAR4M 

209  OMEGA=ZnMEGA(LOM) 

211  CONST=ALPHMN*< 1NCRMT-1 ) * ( ALPHMX-ALPHMN) /899 .0 

212  IF(MODORD-2)260,213,213 

213  00    2H       1  =  1, LOE 
A(  I  )  =  HA(  I) +  HAT! I  )*CONST 

214  B<  I)  =  HBT( I  )*HCT<  I )*CONST 
DO  215  I  =1,M0E 
0(  I  )  =  HD<  II+HDT!  I  )*CONST 

215  E(  I)  =  HET( I )*HFT(  I)*CONST 
GO  TO  263 

216  IF( IMAG-1)217,218,219 

217  DUMMY 1=10.0**(ZMAGSQ(L CM)/ 10.0) 
GO  TO  220 

218  DUMMY l=ZMAGSO(LOM)*ZMAGSO(LnM) 
GO  TO  220 

219  DUMMY1=ZMAGS0IL0M» 

220  IP(NPARAM-2)211,211,212 

253  OMEGA=CONOME 

254  IF(MODORD-2)255. 256, 1000 

255  CONST=ZBETA!LOM) 
GO  TO  ?60 

256  CONST=ZALPHA(LOM) 
GO  TO  213 

257  IF(MOOOPD-2)25<5,258,258 

258  BETA=ZBETA(LOM) 
CONST  =  CONAI.P 
GO    TO    213 

259  8ETA=ZALPHA(L0M) 
C0NST=C0N9ET 

260  DO  261   1=1, LOE 
A(  I )  =  HA(  I1+H8T!  I  )*CONST 

261  B<  I  )  =  HAT(  I  )*-HCT<  I)*CONST 
00  262   1=1, MOE 
D(  I  )  =  HD(  D+HET!  I  )*CONST 

262  E(  I )  =  HDT( I  )*HFT(  I)*CONST 

263  GO  TO  (350,350,275,  290,264  ,29 1  ,294 , 1  000 ) , NP4R AM 

264  OMEGA=ZOMEGA(LOM) 
27  5  ZMAG=XMAGMN+< INCRMT-1)*( XMAGMX-XMAGMNI) /899.0 

IF( IMAG-1)276,277,278 

276  DUMMY1=10.0**(ZMAG/10.0) 
GO  TO  350 

277  DUMMY1=ZMAG*ZMAG 
GO  TO  350 

278  DUMMY1=7MAG 
GO  TO  350 

290  IF( IGRAPH)292,292,294 

291  IF(  (  IGRAPH.GE.C) .AND.(  IGRAPH.LE. 1) )  GH  TO  294 
29  2  OMEGA=OMEGMN+INCRMT*(nMEGMX-OMEGMN)/900.0 

GO  TO  350 
29  4  WINC=1.0+! INCRMT-l)*NUMDEC/899.0 

OMEGA=( YLOGMN* 10. 0**WINC 1/10.0 
C**************\JUMERATCR  COEFFICIENT  CALCULATIONS*********************** 
350  SNB2=0.0 

SN81=0.0 


SNBOO.O 

SUMN=0.0 

ND0UN=2*N0RN 

DO    373     IL03P=1,<. 

00    373    NP=2,L0E 

00     373    mo=1,NORN 

MP=MQ-1 

N0=NP-1 

IF< (NP+MJ)-Nn0UN)3t 1,36 1,373 
3ol     IF< <MP+N3)-NORN» 362,362,373 
362     IF(NP-MP)373, 373,363 
3b3    IF(MP) 364, 36^,366 
364    C=1.C 

GO    TO    367 
36o    C=2.0 

367    GO    TO     ( 36d, 369,370,371 ), ILOOP 
'-.»    CSNB2    =     <C*B(NP+MP  )*B(  NP-MP)  ) *0MEG A** ( 2*N0» *< -1 .0 ) ** ( MQ+ 1  ) 

SNR2  =  SN32  +  CSNP2 

GO  TO  373 

369  CSNB1  =  <C*A(NP  +  MP)*3(  NP-Mp)+c*4(NP-MP)*B(NP+MP)  )*PMFr,A**(  2*NQI*<- 
ll.O)**(MQ*U 

SN81  =  SNfll  ♦  CSNB1 
GO  TO  37  3 

370  CSNBO  =  (C*A(NP  +  MP)*A<  NP-MP)  )*C1MEGA**<  ?*Na)*<-1.0)**<MQ*l  I 
SNBO  =  SNBO  +  CSNRO 

GU  TO  37? 

371  COEFN  =  (C*A(NP  +  MP)*A<NP-MP)+C*A<NP+MP)*R(NP-MP)*BFTA«-C*A(NP-Mp)*e 
1<NP*mP)*BETA+C*H(NP+mP)*B<NP-MP) *BET A* BET  A  I *OMEGA**< 2*NQ ) * < -1 . 0 ) ** 
2<M0*1> 

SUMN  =  SUMN  ♦  COEFN 
373  CONTINUE 

5NB2  =  SNF32  ♦  B(  1)*R<1  ) 
SNB1  =  SN31  ♦  2.0*A(  1) *B< 1) 
SNBO=SNPO*A(  1  )*A(  1  ) 

SUMN  =  SU^N  +  A(  1)*A(  1  )+2.0*A<  1  )  *B(  1  >*8FTA*R(  1  >*R(  1  )*BETA*BETA 
C«*«««»»*»*«»pENriMINATOP  COEFFICIENT  CALCULATIONS*********************** 
STR2=0.0 
SOB1=0. 
S0B0=0. 
SUMD=0. 
N00UD=2*NJRD 
DO  385  IL03P=1,<. 
DO  385  NP=2,M0E 
DO  385  mc=1,N0R0 
MP=MQ-1 
N0=NP-1 
IF  (  (NP*MO)-NufJtn  13  75,  375,  385 

375  |F(  (MP*i*3)-N'0RD)  376,  5  7  6,385 

376  IF(NP-MP) 395,^85,377 

377  IF(«P) 378, 378,379 

378  C=1.0 

GO  TO  38C 
J79  C=2.0 

380  GO  TO  (  381,332,383,33^.1  ,  ILOOP 

381  CS0B2  =   C*E«NP  +  »'F)*Ef  NP-MP)*QMEGA**(2*NO)*(-1.0)**(MO+1) 
SDB2  =  SUB2  ♦  CS0B2 

GO  TO  38  5 
3d2  CSDB1  =  (C*J(NP  +  MP  )*E(NP-MP)+C*D(NP-mpi*e (MP  +  MP)  |*rMEG4**(  2*N0)*(- 
Ll.OI**(MQ+l| 
SOBi  =  SDil  ♦  CSDBl 
GO  TO  385 

383  CSDBO  =  (C*D(NP+MP)*D(NP-MP) )*OM£GA**( 2*N0)*(-I.n»**(MQ+l ) 
SDHO  =  S03D  +  CSDPC 

GO  TO  3Q5 

384  COEFD  =  (C*U(NP«-Mp  )*o<NP-MP)-»-C*Q<NP*MP)*F(MD-MP)  *RET4+C*0(NP-mp  )*F- 
1(NP*M?  )*BETA*C*E  (NP  +  »"P)*F(NP-MP)*BFTA*BETA)*nMFr,A«*(  ?*NQ ) * ( -1 . C ) ** 
2(«0+l ) 

SUMO  =  SU*D  ♦  COEFC 
395  CuNTINUF 

SUMD  =  SU'-«D  +  U(  1)*0(  1  )*2.0*n(  1  )*E(  1  )*BETA  +  F(1  )*F(  1  )*BFTA*BETA 

SDB2  =  SD82  ♦  t(  1)*E( 1 ) 

S0B1  =  SD81  ♦  2.0*D(  l)*E(  1) 

SDBO  =  SOBO  ♦  D( 1)*D( ! ) 

IF(NPARAM-7)',13,  3F6.1C00 
C****»*******************PARAM  7  PHASF  COMPUTATIONAL  SECT  I  "IN****  ******** 
386  T0E=L0E/2.0 

PNEVEN=0. 

PDEVE\=0. 

PN00D=0. 


P0000=0. 

J0E=L0E/2 

IF( JOE-TOE)387,388,388 

387  J0N=(L0E*l)/2 
MON=JON 

GO  TO  389 

388  JON=(LOE+2)/2 
MON=JON-l 
IF(NORN-l) 393,389,389 

389  DO  390  1=2, JON 

390  PNODD=PNCDD+< A( 2*1-2 )  +  8( 2*1-2 ) *BETA)*0«EGA**( 2*1-3 )*<-l. 01 **( I ) 
IF(NORN-l) 393,393,391 

391  00  392  1=2, HON 

39  2  PNEVEN=PNEVEN*(A<2*I-1 H-R < 2* I  -  1 ) *BET A) *OMEGA** < 2* I -2 >* < - 1 . )**< 1-1) 

393  PNEVEN=PNEVEN*A<  1)*B(1)*BETA 
T0E=M0E/2.0 

J0E=M0E/2 

IF< J0E-T0E)394,395,395 

394  JON=< MOE+1 )/2 
MON=JON 

GO  TO  396 

395  JON=(MOE+2)/2 
M0N=J0N-1 
IF(N0RD-1)400,396,396 

396  DO  397  1=2, JON 

39  7  PDO0D=PDO0O+(D< 2*1-2 )+E( 2* I -2 ) *BE TA ) *OMEGA**( 2* I -3 ) * ( -1 . 0 ) ** ( I ) 

IF (NORD- 1)400, 400, 39  8 

398  DO  399  1=2, MON 

399  PDEVEN=PDEVEN*(D< 2*1-1 )  +  E ( 2*1- I ) *BET A > *OMEGA** ( 2* I -2 )* ( -  1 .  )**<  1-1) 

400  PDEVEN=PD£VEN  +  D(  1)+E<1 )*BETA 

TANG=( PDEVEN*PNODD-PNEVEN*POODD) /  <  PNF  VFN*PDEVEN  +  PNODD*PDr)DD  ) 

RAD=ATAN{TANG) 

SINE=SIN(RAD) 

IF( TANG)402,4C2,401 

401  SINE=-SINE 

402  YMAGSU=ABS( SUMN/SUMD) 
YMAG=SQPT( YMAGSQ) 

YMAGDB=20.0*.4  34  2945*AL0G< YMAG) 
YMAGAB=ABS( YMAGOB) 

IF( YMAGPK-YMAGDB) 40  3,404,404 

403  YMAGPK=YMAGDB 
OMEGPK=OMEGA 

404  IF(DBMAX-YMAGDB)405,406,406 

405  YMAGOB=DBMAX 

406  IF( Y*AGDB-QBMIN)407,408,408 

407  YMAGDa=DBMIN 

40  8     IF(  (  XS*8.0)-YMAGAR)  409  ,410,4  10 

409  YMAGAB=8.0*XS 

410  GO    TO    426 

C*******»*******«0UAPRATIC     FORMULA     SOLUTION    OF     DEPENDENT    V  AP  I  4BL  E  ,  P  A<J  AM 
413    AC0=DUMMY1*SDB2-SNB2 
416    BC0=DUMMY1*SDB1-SNB1 

CC0=DUMMY1*SD80-SN80 

IF( ABS(ACO) .LT.l.OE-35)  GO  TO  425 

DPP=BC0*BC0-4.0*AC0*CC0 

IF(OPP.LT.O.O)  GO  TO  425 

DSQR=ZS1GN*SQRT(DPP) 

BETA=(DS0R-BCO)/(2.0*ACO) 

IF < BET A-BETAMN I  425,426,424 
424  IF(BETA-BETAMX)426,426,425 


425    IF(  IP-1>500,505,507 

42  6    IF<  (  IWRITE.LT.l)  .OP. ( IWRITE.GT.< 


900) )  GO  TO  450 
:*«************<<*COMPUTED  DATA  PRINTOUT  SECTION**' 
ISTEP=ISTEP+1 


IF(  IP)428,428,427 
427  IF(  ISTEP-IWRITF)4! 


i50,430,430 

428  WRITE(6,429)L0M 

429  FORMAT! /,10X, •     ALPHA        BETA      MAGNITUDE 
&  TANGENT   ;   CURVF  NUMBER • , I  2 , / ) 

430  IFIM0D0RD.E0.2)  GO  TO  431 

GO  TO  (433,435,436,437,439,441 ,443,1000) .NPARAM 

431  GO  TO  < 432,434,445,438,440,442,444, 1000) , NPARAM 

432  WRITE(6,6)    CONST  ,  BETA , CONMAG , OMFGA 
GO  TO  44  8 

433  WRITEI6.6)    BETA, CONS T , CONMAG , OMFGA 
GO  TO  44  8 

434  WRITE(6,6)    CONST, BETA, ZMAGSQ(LOM) , OMEGA 
GO  TO  448 

435  WRITE(6,6)    BE T A .CONST , ZMAGSQ < LOM) , OMEGA 


GO  TO  44  8 
43t.  WRITEI6.6)    BETA, CONST, ZMAG, OMEGA 
GO  TO  448 

437  WRITE16.6)    BETA, CONrtET.ZMAGSQILCM)  .OMEGA 
GO  TO  448 

438  WRITE<6,6)    CONAL P , BE T A , ZMAGSQ! LOM) , OMEG A 
GO  TO  448 

439  WRITt(6,6)    BE TA.CONBET , ZMAG, OMEGA 
GO  TO  448 

440  WRITEI6.6)    CONALP, BETA, ZMAG, OMEGA 
GJ  TO  44  8 

441  WRITE!6,6>    BE T A  ,  C ONS T , CCNMAG , OMEGA 
GO  TO  448 

442  WRITE<6,6)    CONST, BETA, CONMAG.OMFGA 
GO  TO  448 

-.43  wRITE(o,6)   ZALPHAILC") ,CONBET,YMAGDB, OMEGA, TANG 

GO  TO  44  8 
<".'<  WRITEI6.6)   CONALP,ZBETA(LOM)  .YMAGDB, OMEGA, TANG 

GO  TO  44  8 

443  wRITE!6,6)    CONST, BETA, ZMAG, OMEGA 
44e  ISTEP=0 

C**»»»*«**«*»***»«*»GRAPHICAL  ARRAY  STORAGE  SECTION********************* 
450  IP=IP+1 

Y(IP)=  BETA 

452  GO  TO  (453,453,460,470,460,470,475  ,  1000 )  , NP AR A M 

453  X( IP)=  CONST 
GO  TO  5C0 

450  X( IP)=ZMAG 

GO  TO  5C0 
470  X( IP)=JMPGA 

IF(NP<WAM- 6)471, 472, 475 
h71     IF( 1GRAPH) 500,500,474 
472     IF<  (  IGRAPH.LT.O)  .OP.  ( IGRAPH.GT. 1) )     GO    TO    500 

474  X( IP)=LX-BETA 
GO    TO    479 

475  X(  IP)=-YM4GPB 

UILOM, IP)=U(LOM, ip)+Y«AG 

VILOM, IP)=( ALOG10(CMEGA/YLOGMN)*LY/NUM0ECI*YS 

S(LCM,IP)=SINE*YMAGA8/(XS*TANG) 

T(LO",  IP)=SINF*Y*AGA8/XS 

IC(LOM)=  IP 

IF(  10UAD) 476,479,478 

476  b.(LOM,  IP)  =  -4.0*RAD/1.570ft 
IF( RA01477.479, 479 

477  W(LOM, IP )=8.0-U( LOM, IP) 
GO    TO    479 

47  8    rflLOM, IP)=XS*(LX/2.-lYI-RAD*XS*LX/3.l4l4 

IF(  (  IP)  .LE.I  1) )G0    TO    *79 

IM  <  YMAGD*).NE.(  CHMIM  )     GO    TO    479 

IF(  (WILOm,  IP  )  )  .GE.IWlLOM,  IP-1)  )  )G0    TO    47-3 

W(LOM, IP)=W(LUM, IP-1 ) 

V<  LOM,  ipj=v<LOM,  IP-1  ) 
4  79  Y( IPI=( ALOGlOITMFGt/Yl CGMN ) * L Y/ NUMOEC ) *Y S 

500  IF(  INCRMT.LT.900)  GO  TC  518 

r  »«•»»«*»*■************«****<-,  RAPHlCAL  CURVE  PL  HTTING  ***♦*****#*********** 
IF(  IP-1)  =01,505,507 

501  IF( JCURVF.GT.O)  GO  TO  503 
WRITE (6,  502)  LOM 

502  FORMAT!  10X, 'NO  GRAPH  POINTS  WERF  FOUND  FOR  CURVE  NUMF)FR  ',12,  /I 
GO  TO  518 

503  wRITE(6,  504)L(IM 

5„4  FORMAT! lCX, 'NO  FURTHER  POINTS  WERE  FOUND  FOR  CURVF  N).«,I3,/I 
GO  TO  518 

505  v.RITE(6,506) 

506  FOPMAT! 1CX, *ONL Y  OWE  POINT  WAS  FOUND,  AND  WILL  NOT  P^  "LOTTED',/) 
GO  TO  516 

507  J=L0"+50 
JCURVE=1 

IF! IGRAPH.LE.-2)  GC  TO  516 

IF!  (NPAPAM.EU.  7)  ,A*0.(  IPAR7.LT.  NPAR7)  )  GO    TO    513 
IFUZSIGN.GT. 0.0). AND. (JCURVF.GT.O))  J=50 
WPITEI6.51  DLPM 

511  FjRMAT! /,10X, 'THE    PREVIOUS    POINTS    WIIL     BE    PLOTTED    t-S    CUPVF     ',!?,/) 
IF(NP ARAM- 7)514, 512, 514 

512  WRITE(6,513)YMA;,FK  ,0«FGPK 

513  FORMAT! /,10X, -THE  MAXIMUM  VALUE  OF  MAGNITUDE  IS'.FIO.^,'  LOCATED  A 
61  OMEGA  EQUAL  TO',F10.4,/) 

J=50 

514  CALL  CRAW!  I  P  ,  X , Y  ,  MCD  ,0 ,L A ( J )  ,IT1TLE,XS,YS,IX,IY,2,2,LX,LY,1,LAST) 


516 

518 


519 
520 


M0D  =  2 
IP=0 

CONTINUE 

IFINPARAM.EQ.7)  GO  TO  520 
CONTINUE 
CONTINUE 

IF( ( IPAR7.GE.NPAR7).0R.(NPARAM.NE.7) >  GO  TO  529 
'♦♦GAMMA  &  PARAM-7  SUPERPOSITION  :  SEE  IN°UT  DATA  SYMBOL  NGAMMA**** 
IF( (NGAMMA.LE.O) .OR. ( I  GAMMA. GT.NGAMMA) )  GAMMA=0.0 
IF(  IPAR7.GT.DG0  TO  523 
00  521  1  =  1, LOE 
0L0HA(  I  )=HA(  I ) 
OLDHATI  I  )  =  HAT( I  ) 
OL0H8T(  I)  =  HBT( I ) 
OLDHCT<  I  )  =  HCT(  I ) 
DO  522  1=1, MOE 
OLDHO  (  I  )  =  HD  < I ) 
OLDHDTt I )=H0T( I ) 
OLDHET( I )=HET< I ) 
OLDHFTI I )=HFT( I ) 
IPAR7=IPAR7*1 
IF(NGAMMA.GT.O)  GO 
READ(5,5)  (HA(I).I^ 
REA0(5,5)  (HAT(I), 
PEA0(5,5)  (HBT(I), 
READ<5,5)  <HCT(I), 
KEAD(5,5)  (HD  ( I ), 
READ(5,5)  (HDT(I), 
READ(5,5)  (HET(I) 
1FT(  I  ) 


530 
531 
532 


539 
542 


0  526 

,LOE) 

=1,L0E) 

=1,L0E) 

=1,L0E) 

=1,M0E) 

=1,M0E) 

=1,M0E) 
READ(5,5>  (HFT(I),  1=1,  MOE) 
IF(  IPAR7.LT. NPAR7)  GO  TO  200 
DO  525  L0M=1,NLIMIT 
IP=IC(LOM) 
00  524  1=1, IP 
XI  I  »=-U(LOM,I ) 
Y(  I  )  =  V(LOM,I I 
J=L0M+1 

CALL  DRAW(  I P , X , Y , MCD ,Q ,L A < J ) ,ITITLE,XS,YS,IX,IY,2,?,LXfLY,l,LAST| 
00  527  1=1, LOE 

HA  (I)=OLOHA  ( I  )+GAMMA*GAMHA  (I) 
HAT( I )=CLOHAT< I ) +GAMMA*G AMHAT ( I ) 
HBT( I )=0L0H8T( I ) +G AMMA*G AMHB T ( I ) 
HCT( I )=OLDHCT< I )*GAMMA*GAMHCT< I ) 
DO  528  1  =  1, MOE 

HD  (I)=OLDHD  ( I )+GAMMA*GAMHD  (I) 
HDT( I )-OLOHDT( I )+GAMMA*GAMHDT< I ) 
HET(  I  )=OLOHET(  I  )  *GAMMA*G  AMHET(  I  ) 
HFT(  I  )=OL0HFT(  I  ) +G AMMA *G AMHF T<  I ) 
IF(  IGAMMA.LE.NGAMKA)  GO  TO  26 

IF< (NPAR AM.NE.7)  .OR. < IGRAPH.LE.-2) )  GO  TO  999 

*******#****#*******#p  AR  A"  7  PHASE  CURVES**********'***************' 
IPAR7=2 

IF(  IPHASE) 537,532,999 
DO  531  1=1,12 
ITITLEJ  I)  =  JTITLE( I  ) 
DO  536  L0*=1,NLIMIT 
I°=IC(LOM) 
00  533  1=1, IP 
X(  I  )=W(LOM, I ) 
Y(  I  )  =  V(LOM,I ) 

CALL    DRAM  IP  ,  X , Y  ,MCD ,0 , L A ( 50),JTITLE,XS,1.,0,IY,2,?,LX,LY,1,LAST) 
IPAR7=l 

IF(  I  IPOLAR).LE.( 0) )     GO    TC    999 

*********************P4RAM    7    POLAR    PI  OT    CURVES  ******************** 
M0D=1 

DO  542  1=1,8 
RADIUS=I/2.0 
DO  539  J=l,360 
RADIAN=6.28318/360*J 
X( J  )  =  RA0IUS*C0S( RADIAN) 
Y( J)=RADIUS*SIN(RADI AN) 

CALL  DRAW( J,X,Y,MOD,0,LA( 1 ) , KT I TLE , 1 . , I . ,4,4, 2 , 2 , 8 ,8  ,  0,  L AS T ) 
M00=2 

DO  543  J=10, 180,10 
X(  1)  =  X( J) 
X(2)=-X( 1) 
Y( 1)=Y( J) 
Y(2)=-Y(l) 
CALL  DRAW(2,X,Y,2,0,LA<l),KTITLF,l.,l.,4,4,?,2,R,q,0,LAST) 


00  545  L0M=1,NLIMIT 
[P=IC(LOM| 

DO  544  J  =  1,IP 
X( J  )  =  S(LOM,J) 

544  Y< J)=T(LOM, J) 

L  =  LOM«.l 

545  CALL    DRAM  I  P  ,  X  ,  Y  ,  2  ,0  ,L  A(  L)  .KTITLE.l.  ,1  .  ,4 ,4 ,2 , 2,8 , 8 , 0, L AST  ) 
<?^9     IF( <MOO).NE.(?> )     GC    TO    1000 

X(  I)     =     0.0 
Y(  1  )     =    CO 
XI2)     =     <LX-IY)*XS 
Y(  2)     =     CO 

CALL    OK  Aw    <2,X,Y,3,O,LA(50),ITITLE,1.0,l.O,O,d,2,2,LX,LY,l,LA?T| 
"00=1 

IF(  (  IPJR7I  .NE.(  I )  )  GC  TO  55 
1000  CONTINUE 


STOP 
tNO 

C*»»*«**»*SAMPLE    DATA    DE C* ****************** **** 

********** 

********* 

1 

4 

1 

2 

4 

6 

0 

3 

1 

8 

R 

8 

0 

10 

-1 

1 

I 

0.00047 

C036 

12.0 

10. 

-0.001 

0.007 

CO 

0.05 

-3.0 

n.o 

0.0 

-0.001 

0.0 

0.000*7 

0.002 

0.003 

0.004 

0.006 

0.0 

0.01 

0.02 

0.03 

0.036 

0.04 

0.045 

-3.0 

0.0 

3.C 

6.0 

CO 

12.0 

15.0 

1.0 

4.0 

ICO 

12.0 

15.0 

20.0 

30.0 

45.                  50. 

55. 

60. 

70. 

BO. 

100000.0 

(  BlAmk    cflto 

) 

fit.***       C««pN 

Cfii»«K    cfl^r' 

100000. 

5C0.C 

60. C 

1.0 

1000C0.0 
„  ,100000.0 

GLAVIS,     PAR4M     1,     3H0    OHLFR     T/A    F.R.,     OMEGA    CURVS 
0R0INATE  =  ALPHA,     Ar,C  I  S  S  A  =  8E  T  A  ,     C  ON*  AGN  I  TUDE  =  1  20« 
GLAVIS,     PARAM    2,     1KD    OPDER     T/A    F.R.,     MAG    CURVFS 
0RCINATE=ALPHC     4BCISSA=6ETA,    CONOMeG A=l OR AO /SEC 
GLAVIS,     PARA*    c,     -JBD    HRUER     T/A    F.C,     MAG    CURVFS 
Ofi^INATE=AL°HA,     ARC  I S S A=r*EG A ,     CONBET A  =  0 . 0 36 
GLAVIS,     PARAM    *,     IPO    ORDER     T/A    CB.,     BETA     CURVFS 
ORDlNATE  =  AL"HA,     ARC  I S S A=CM£GA ,     CONMAGN I TUO E= I  2  OH 


LOGIC  OUTLINE  OF  PROGRAM  PARAMS 

1.  Read  number  of  runs  submitted. 

2.  Set  run  number  equal  to  1. 

3.  Read  input  data  for  appropriate  run. 

4.  Print  input  data  for  appropriate  run,  if  requested. 

5.  Calculate  graph  scales  and  location  of  origin. 

6.  Set  NPARAM  equal  to  the  first  PARAM  value. 

7.  Select  initial  values  for  appropriate  PARAM. 

8.  Print  graphical  information  for  appropriate  PARAM,  if  requested. 

9.  Construct  logarithmic  grid  if  required. 

10.  Set  curve  number  equal  to  1. 

11.  Select  curve  value  for  appropriate  curve  number. 

12.  Set  abcissa  variable  to  minimum  allowable  value. 

13.  Set  number  of  graph  points  equal  to  0. 

14.  Compute  numerator  and  denominator  polynomial  coefficients. 

15.  Check  for  possible  division  by  0. 

16.  Calculate  ordinate  value,  using  the  quadratic  formula  for  PARAM-1 
through  PARAM- 6. 

17.  Check  if  calculated  ordinate  value  is  within  graph  limits.   If  yes, 
continue;  if  no,  go  to  statement  21. 

18.  Increase  number  of  points  by  1. 

19.  Print  computed  values. 

20.  Store  values  necessary  plotting. 

21.  Increase  abcissa  variable  by  pre -computed  increment. 

22.  Check  if  abcissa  variable  exceeds  maximum  graph  limit.   If  yes, 
continue;  if  no,  go  to  statement  14. 


176 


23.  Check  if  any  points  were  found  within  graph  limits.   If  none, 
so  state  and  continue  without  plotting.   If  one,  so  state 

and  continue  without  plotting.   If  more  than  one,  continue  and 
plot  the  graph  points. 

24.  Check  if  both  positive  and  negative  coefficients  were  used  in  the 
quadratic  formula  solution  where  applicable.   If  yes,  continue; 
if  no,  change  sign  and  return  to  statement  12. 

25.  Increase  curve  number  by  1. 

26.  Check  if  curve  number  exceeds  the  maximum  number  of  curves  stipulated 
in  input  data.   If  yes,  continue;  if  no,  return  to  statement  11. 

27.  Draw  graph  and  plot  computed  curve  data  on  same. 

28.  Increase  parameter  number  by  1. 

29.  Check  if  parameter  number  exceeds  maximum.   If  yes,  continue;  if  no, 
return  to  statement  7. 

30.  Increase  run  number  by  1. 

31.  Check  if  run  number  exceeds  maximum.   If  no,  return  to  statement  3; 
if  yes,  stop. 
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FLOW  CHART  FOR  PARAMETER  SIX 


(Print  Input  Data) 
(CalculateX+Y  Scales) 


:(Print  Graphical  Information) 


Magnitude 

-Squared < 

Actual 


Magnitude 
Actual 


(Magnitude  =  EXP  (Magnitude/8T5B~5g9)) 


(Magnitude  =  Magnitude  Q 


Digraphs  ^~o~~) 


(  jac  =  o  ) 
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(  DO   ICURVE  =   1,    NUMCURVES^) 


a/ 


(&MEGA  »  OMEGAMIN") 

^^ 

< =-< 


0URVE  =ZBETA  (ICURVE)) 


C 


CURVE  =  ZALPHA(ICURVE>) 


Calculate  Coefficient 
Variables 


A  A 


{ 


Calculate  Coefficient  ] 
Variabl  e  s J 


("IPOINT  =  ~0~") 


GRAPH  -1VL_^. 


OMEGA=OMEGA  +  SQUARE  GRID 
!     INCREMENT 


OMEGA^OMEGA+LOG  GRID 
INCREMENT 


A  A 


A  A 


"OMEGA- 


1 


■OMEGAMSX^ 


Compute  Numerator  coef ficientsj) 


Compute  Denominator  Coefficients' ^ 


0®(3)& 


I 


a  ^ 


<v  A 


Compute  A  for  S 


I2  -  4AC 


2A 
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A  A  A  A 


&=■&*■     -4Ad 


0(J>(3)© 


A  *  A 


(S)SQR  =  SIGN*SQRT(g) 


Cbeta  =  (-b  +  pswrgfo 


(Y( IPO INT)   »  BETA) 


<§> 


-£-- 
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0<£>(D 


AAA 


® 


J>RINT"NO  POINTS  WERE  FOUND")   *    (PRINT  "POINTS  WILL  BE  PLOTTED^ 


PRINT  "ONLY  1  POINT  WAS  FOUND", 

i        


SIGN  =    -SIGN 


^CONTINUE      ) 


(f~       END       ~~) 
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Algebraic  Methods  are  employed  for  frequency  response  studies  of  physical 
systems.   The  transfer  function  of  each  system's  mathematical  model  is  composed 
of  at  least  two  independent  variables  embedded  in  the  coefficients.   These, 
together  with  magnitude,  angular  frequency,  and  phase  angle  are  the  variable 
parameters  to  be  studied. 

Computer  programs  are  developed  for  graphically  representing  these 
variable  parameters.   Interpretation  of  these  results  is  made  to  evaluate  the 
utility  of  such  methods  in  design  and  analysis  studies. 
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